Historical Orthoimage : Différence entre versions

De MicMac
Aller à : navigation, rechercher
(Description)
Ligne 12 : Ligne 12 :
  
 
=Tutorial=
 
=Tutorial=
==Preamble==
 
MicMac use EXIF metadat in order to determine image format and focal length. However, historical images doesn't have such metadata, so we have first to create a file <i>MicMac-LocalChantierDescripteur.xml</i> like this :
 
  
 
==Internal Orientation==
 
==Internal Orientation==
 +
MicMac use EXIF metadat in order to determine image format and focal length. However, historical images often don't have such metadata, so we have first to create a file <i>MicMac-LocalChantierDescripteur.xml</i> like this :
 +
 +
<pre>
 +
 +
<Global>
 +
  <ChantierDescripteur >
 +
 +
<!-- Define a camera model (name and sensor/film size) -->
 +
    <LocCamDataBase>
 +
        <CameraEntry>
 +
              <Name> ZeissRMKATOP15  </Name>
 +
              <SzCaptMm>  226.004 226.008  </SzCaptMm> <!-- MidSideFiducials or "MaxFidX-MinFidX MaxFidY-MinFidY"-->
 +
              <ShortName> Zeiss RMK A Top15* and Zeiss Pleogon A3/4 </ShortName>
 +
        </CameraEntry>
 +
    </LocCamDataBase>
 +
 +
<!-- Associate images to a camera model -->
 +
    <KeyedNamesAssociations>
 +
            <Calcs>
 +
                <Arrite>  1 1 </Arrite>
 +
                <Direct>
 +
                      <PatternTransform> .*  </PatternTransform> <!-- Regular expressions of the group of images with the following camera model -->
 +
                      <CalcName> ZeissRMKATOP15 </CalcName> <!-- Name of the camera for these images -->
 +
                </Direct>
 +
            </Calcs>
 +
            <Key>  NKS-Assoc-STD-CAM </Key>
 +
    </KeyedNamesAssociations>
 +
 +
<!-- Associate images to a focal length -->
 +
    <KeyedNamesAssociations>
 +
            <Calcs>
 +
                <Arrite>  1 1 </Arrite>
 +
                <Direct>
 +
                      <PatternTransform> .*  </PatternTransform> <!-- Regular expressions of the group of images with the following focal length -->
 +
                      <CalcName> 153.664 </CalcName> <!-- See calibration report -->
 +
                </Direct>
 +
            </Calcs>
 +
            <Key>  NKS-Assoc-STD-FOC  </Key>
 +
    </KeyedNamesAssociations>
 +
 +
  </ChantierDescripteur>
 +
</Global>
 +
<pre>
 +
  
 
<pre>mm3d SaisieAppuisInitQT "image.tif" NONE id_appuis.txt image.tif.xml </pre>
 
<pre>mm3d SaisieAppuisInitQT "image.tif" NONE id_appuis.txt image.tif.xml </pre>
 
<pre>mm3d Kugelhupf 1987.*tif Ori-InterneScan/image.tif.xml SearchIncertitude= ??</pre>
 
<pre>mm3d Kugelhupf 1987.*tif Ori-InterneScan/image.tif.xml SearchIncertitude= ??</pre>
 +
 +
== External orientation ==
 +
 
<pre>mm3d OriConvert OriTxtInFile GPS_sommets.txt Sommets NameCple=Couples.xml</pre>
 
<pre>mm3d OriConvert OriTxtInFile GPS_sommets.txt Sommets NameCple=Couples.xml</pre>
 
<pre>mm3d Tapioca File Couples.xml 2000</pre>
 
<pre>mm3d Tapioca File Couples.xml 2000</pre>
Ligne 35 : Ligne 80 :
 
<pre>mm3d AperiCloud 1987.*tif TerrainFinal SH=HomolMasqFiltered</pre>
 
<pre>mm3d AperiCloud 1987.*tif TerrainFinal SH=HomolMasqFiltered</pre>
 
<pre>meshlab AperiCloud_TerrainFinal.ply</pre>
 
<pre>meshlab AperiCloud_TerrainFinal.ply</pre>
 +
 +
== DEM processing and orthorectification==
 +
 
<pre>mm3d Tarama 1987.*tif TerrainFinal
 
<pre>mm3d Tarama 1987.*tif TerrainFinal
 
mm3d SaisieMasqQT TA/TA_LeChantirer.tif</pre>
 
mm3d SaisieMasqQT TA/TA_LeChantirer.tif</pre>

Version du 2 janvier 2017 à 15:10

Picto-liste.png Tutorials index

Description

This tutorial will present the method to process historical aerial images into DEM and Orthoimages. With this kind of products, you can monitor changes in an arean (urbanization, landscape changes, etc...).

The USGS NAPP program offers a large amount of free scanned images over the continental US (mostly), often with calibration data, though the Earth Explorer.

If you are looking for a special area in France, you can use the Geoportail (IGN) to download your own images and process it.

Download

Presentation

Tutorial

Internal Orientation

MicMac use EXIF metadat in order to determine image format and focal length. However, historical images often don't have such metadata, so we have first to create a file MicMac-LocalChantierDescripteur.xml like this :


<Global>
  <ChantierDescripteur >

	<!-- Define a camera model (name and sensor/film size) -->
    <LocCamDataBase>
        <CameraEntry>
              <Name> ZeissRMKATOP15  </Name>
              <SzCaptMm>  226.004 226.008  </SzCaptMm> <!-- MidSideFiducials or "MaxFidX-MinFidX MaxFidY-MinFidY"-->
              <ShortName> Zeiss RMK A Top15* and Zeiss Pleogon A3/4 </ShortName>
         </CameraEntry>
    </LocCamDataBase>

	<!-- Associate images to a camera model -->
    <KeyedNamesAssociations>
            <Calcs>
                 <Arrite>  1 1 </Arrite>
                 <Direct>
                       <PatternTransform> .*  </PatternTransform> <!-- Regular expressions of the group of images with the following camera model -->
                       <CalcName> ZeissRMKATOP15 </CalcName> <!-- Name of the camera for these images -->
                 </Direct>
             </Calcs>
             <Key>   NKS-Assoc-STD-CAM </Key>
    </KeyedNamesAssociations>
	
	<!-- Associate images to a focal length -->
    <KeyedNamesAssociations>
            <Calcs>
                 <Arrite>  1 1 </Arrite>
                 <Direct>
                       <PatternTransform> .*  </PatternTransform> <!-- Regular expressions of the group of images with the following focal length -->
                       <CalcName> 153.664 </CalcName>	<!-- See calibration report -->
                 </Direct>
             </Calcs>
             <Key>   NKS-Assoc-STD-FOC  </Key>
    </KeyedNamesAssociations>
	
  </ChantierDescripteur>
</Global>
<pre>


<pre>mm3d SaisieAppuisInitQT "image.tif" NONE id_appuis.txt image.tif.xml 
mm3d Kugelhupf 1987.*tif Ori-InterneScan/image.tif.xml SearchIncertitude= ??

External orientation

mm3d OriConvert OriTxtInFile GPS_sommets.txt Sommets NameCple=Couples.xml
mm3d Tapioca File Couples.xml 2000
mm3d SaisieMasqQT "image.tif"
mm3d HomolFilterMasq 1987.*tif GlobalMasq=filtre.tif
mm3d Tapas RadialBasic 1987.*tif InCal=CalibInit Out=Relative SH=HomolMasqFiltered LibFoc=0
mm3d AperiCloud 1987.*tif Relative SH=HomolMasqFiltered
meshlab AperiCloud_Relative.ply
mm3d SaisieAppuisInitQT 1987_FR4053_07.*tif Relative id_appuis.txt MesuresInit.xml
mm3d SaisieAppuisInitQT 1987_FR4053_08.*tif Relative id_appuis.txt MesuresInit.xml
mm3d SaisieAppuisInitQT 1987_FR4053_09.*tif Relative id_appuis.txt MesuresInit.xml
mm3d GCPBascule 1987.*tif Relative TerrainInit Mesures_appuis.xml  MesuresInit-S2D.xml
mm3d SaisieAppuisPredicQT 1987.*tif TerrainInit Mesures_appuis.xml MesuresFinales.xml
mm3d GCPBascule 1987.*tif TerrainInit TerrainBrut Mesures_appuis.xml  MesuresFinales-S2D.xml
mm3d Campari 1987.*tif TerrainBrut TerrainFinal GCP=[Mesures_appuis.xml,0.1,MesuresFinales-S2D.xml,0.5]<pre>
<pre>mm3d AperiCloud 1987.*tif TerrainFinal SH=HomolMasqFiltered
meshlab AperiCloud_TerrainFinal.ply

DEM processing and orthorectification

mm3d Tarama 1987.*tif TerrainFinal
mm3d SaisieMasqQT TA/TA_LeChantirer.tif
mm3d Malt Ortho 1987.*tif TerrainFinal MasqImGlob=filtre.tif NbVI=2 ZoomF=1 ResolTerrain=0.5 DefCor=0 CostTrans=4
mm3d GrShade  Z_Num8_DeZoom2_STD-MALT.tif ModeOmbre=IgnE Out=ModeleOmbre.tif Mask=MEC-Malt/AutoMask_STD-MALT_Num_7.tif
mm3d to8Bits MEC-Malt/Z_Num8_DeZoom2_STD-MALT.tif Out=hypso.tif Coul=1 Dyn=3 Mask=MEC-Malt/AutoMask_STD-MALT_Num_7.tif
mm3d Tawny Ortho-MEC-Malt Out=Mosaique.tif