Historical Orthoimage - Historique des versions

Luc Girod : /* Relative orientation */

2022-09-14T12:41:23Z

‎Relative orientation

Luc Girod le 14 septembre 2022 à 12:37

2022-09-14T12:37:21Z

Luc Girod le 9 septembre 2022 à 06:06

2022-09-09T06:06:27Z

Luc Girod : /* DEM processing and orthorectification */

2019-11-18T14:51:03Z

‎DEM processing and orthorectification

Luc Girod : /* Relative orientation */

2018-08-13T07:46:09Z

‎Relative orientation

Luc Girod : /* Absolute orientation */ Add SH=HomolMasqFiltered to Campari cmd

2017-03-31T14:35:12Z

‎Absolute orientation : Add SH=HomolMasqFiltered to Campari cmd

Luc Girod : /* Relative orientation */

2017-03-28T07:43:04Z

‎Relative orientation

Luc Girod : /* Absolute orientation */

2017-02-16T17:51:59Z

‎Absolute orientation

Luc Girod : /* Description */

2017-02-16T17:31:20Z

‎Description

Luc Girod : /* Absolute orientation */

2017-02-12T19:44:08Z

‎Absolute orientation

@@ Ligne 159 : / Ligne 159 : @@
 To be able to ignore the fiducial marks and other inscriptions on the images that would yield nonsensical tie points, a mask need to be created.
-<pre>mm3d SaisieMasq "OIS-image.tif"</pre>
+<pre>mm3d SaisieMasq "OIS-AN_IMAGE_NAME.tif"</pre>
 Once created, the mask should be renamed <i>filtre.tif</i>.

@@ Ligne 113 : / Ligne 113 : @@
 </div>
-Then to input the image coordinate of the fiducial marks, you should use the SaisieAppuisInitQT command on each image like this (<i>id_fiducial.txt</i> is a text file with a point name on each line, see example bellow. BE CAREFULL TO NAME THE POINTS THE SAME WAY IN THE <i>MeasuresCamera.xml</i> and the <i>id_fiducial.txt</i>).
+Then to input the image coordinate of the fiducial marks, you should use the SaisieAppuisInit command on each image like this (<i>id_fiducial.txt</i> is a text file with a point name on each line, see example bellow. BE CAREFULL TO NAME THE POINTS THE SAME WAY IN THE <i>MeasuresCamera.xml</i> and the <i>id_fiducial.txt</i>).
 <div class="toccolours mw-collapsible mw-collapsed" style="background-color: Lavender">
 <h6 style="font-family: Helvetica:font-size: 40px">Example of <i>id_fiducial.txt</i> -->Click Expand to reveal--></h6>
@@ Ligne 130 : / Ligne 130 : @@
 </div>
-<pre>mm3d SaisieAppuisInitQT "image.tif" NONE id_fiducial.txt MeasuresIm-image.tif.xml </pre>
+<pre>mm3d SaisieAppuisInit "AN_IMAGE_NAME.tif" NONE id_fiducial.txt MeasuresIm-AN_IMAGE_NAME.tif.xml </pre>
-The resulting <i>MeasuresIm-image.tif-S2D.xml</i> file (created in the image folder) should be moved in the <i>Ori-InterneScan</i> directory and renamed <i>MeasuresIm-image.tif.xml</i> (note the removal of "-S2D").
+The resulting <i>MeasuresIm-AN_IMAGE_NAME.tif-S2D.xml</i> file (created in the image folder) should be moved in the <i>Ori-InterneScan</i> directory and renamed <i>MeasuresIm-AN_IMAGE_NAME.tif.xml</i> (note the removal of "-S2D").
-If you have images where the fiducial marks are easily recognizable (they look like targets, not just a dot), and if the images are already close to be aligned, you can use Kugelhupf to compute the position of the points starting with the second image (appearance and position of the points are dictated by the first image that you processed manually).
+If you have images where the fiducial marks are easily recognizable (they look like targets, not just a dot), and if the images are already close to be aligned (close to the value of SearchIncertitude, in pixels), you can use Kugelhupf to compute the position of the points starting with the second image (appearance and position of the points are dictated by the first image that you processed manually).
-<pre>mm3d Kugelhupf .*tif Ori-InterneScan/MeasuresIm-image.tif.xml SearchIncertitude= ??</pre>
+<pre>mm3d Kugelhupf .*tif Ori-InterneScan/MeasuresIm-AN_IMAGE_NAME.tif.xml SearchIncertitude= ??</pre>
 Then all the images can be re-sampled to fit in the same geometry and can therefore be processed like digital images. The user need to input the scan resolution (in the example line, 0.025 -> 0.025mm=25microns). This process is slow (ca. a minute per image), but is parallelized. Note the use of the '''BoxCh''' option, which is critical to ensure that the final image is indeed the size described in the ''<SzCaptMm>'' tag of <i>MicMac-LocalChantierDescripteur.xml</i>.
-<pre>mm3d ReSampFid ".*.tif" 0.025 BoxCh=[0,0,223,223]</pre>
+<pre>mm3d ReSampFid ".*.tif" 0.025 BoxCh=[0,0,226.004,226.008]</pre>
 The user should then move the original images to a sub-folder, or state OIS.*.tif as the regular expression in futur steps.
@@ Ligne 159 : / Ligne 159 : @@
 To be able to ignore the fiducial marks and other inscriptions on the images that would yield nonsensical tie points, a mask need to be created.
-<pre>mm3d SaisieMasqQT "OIS-image.tif"</pre>
+<pre>mm3d SaisieMasq "OIS-image.tif"</pre>
 Once created, the mask should be renamed <i>filtre.tif</i>.
@@ Ligne 209 : / Ligne 209 : @@
 </div>
-The points can be given in any type of coordinates and be transposed, but MicMac requires a somehow euclidean system to work with (not Lat Long, but Easting Northing).
+The points can be given in any type of coordinates and be transposed, but MicMac requires a somehow Euclidean system to work with (not Lat Long, but Easting Northing).
-To convert the poinst in MicMac xml format:
+To convert the points in MicMac xml format:
 <pre>mm3d GCPConvert AppInFile GCPs.txt</pre>
@@ Ligne 263 : / Ligne 263 : @@
 </div>
-<pre>mm3d SaisieAppuisInitQT "OIS-image1.tif" Relative id_GCPs.txt MeasuresInit.xml
+<pre>mm3d SaisieAppuisInit "OIS-image1.tif" Relative id_GCPs.txt MeasuresInit.xml
-mm3d SaisieAppuisInitQT "OIS-image2.tif" Relative id_GCPs.txt MeasuresInit.xml
+mm3d SaisieAppuisInit "OIS-image2.tif" Relative id_GCPs.txt MeasuresInit.xml
-mm3d SaisieAppuisInitQT "OIS-image3.tif" Relative id_GCPs.txt MeasuresInit.xml</pre>
+mm3d SaisieAppuisInit "OIS-image3.tif" Relative id_GCPs.txt MeasuresInit.xml</pre>
 Use these points to get into the cartographic coordinate of the GCPs.
@@ Ligne 271 : / Ligne 271 : @@
 Then input all the GCPs on all the images (or at least quite a few), using pre-pointed approximate GCPs:
-<pre>mm3d SaisieAppuisPredicQT "OIS-imageN.tif" TerrainInit GCPs.xml MeasuresFinales.xml</pre>
+<pre>mm3d SaisieAppuisPredic "OIS-imageN.tif" TerrainInit GCPs.xml MeasuresFinales.xml</pre>
 Again, use these points (now more numerous) to get into the cartographic coordinate of the GCPs.
@@ Ligne 287 : / Ligne 287 : @@
 Create a pseudo orthoimage (with a "flat" terrain as target) to be able to draw a mask on the approximate area of interest.
 <pre>mm3d Tarama "OIS-.*tif" TerrainFinal
-mm3d SaisieMasqQT TA/TA_LeChantier.tif</pre>
+mm3d SaisieMasq TA/TA_LeChantier.tif</pre>
 Compute the DEM (DEM is the file called ''MEC-Malt/Z_Num7_DeZoom2_STD-MALT.tif'' , or similar).
@@ Ligne 293 : / Ligne 293 : @@
 Once the DEM is generated, other products can be generated :
 *An Orthoimage:

@@ Ligne 138 : / Ligne 138 : @@
 <pre>mm3d Kugelhupf .*tif Ori-InterneScan/MeasuresIm-image.tif.xml SearchIncertitude= ??</pre>
-Then all the images can be re-sampled to fit in the same geometry and can therefore be processed like digital images. The user need to input the scan resolution (in the example line, 0.025 -> 0.025mm=25microns). This process is slow (ca. a minute per image), but is parallelised.
+Then all the images can be re-sampled to fit in the same geometry and can therefore be processed like digital images. The user need to input the scan resolution (in the example line, 0.025 -> 0.025mm=25microns). This process is slow (ca. a minute per image), but is parallelized. Note the use of the '''BoxCh''' option, which is critical to ensure that the final image is indeed the size described in the ''<SzCaptMm>'' tag of <i>MicMac-LocalChantierDescripteur.xml</i>.
-<pre>mm3d ReSampFid ".*.tif" 0.025</pre>
+<pre>mm3d ReSampFid ".*.tif" 0.025 BoxCh=[0,0,223,223]</pre>
 The user should then move the original images to a sub-folder, or state OIS.*.tif as the regular expression in futur steps.

@@ Ligne 290 : / Ligne 290 : @@
 Compute the DEM (DEM is the file called ''MEC-Malt/Z_Num7_DeZoom2_STD-MALT.tif'' , or similar).
-<pre>mm3d Malt Ortho "OIS-.*tif" TerrainFinal MasqImGlob=filtre.tif NbVI=2 ZoomF=2 ResolTerrain=0.5 DefCor=0 CostTrans=4</pre>
+<pre>mm3d Malt Ortho "OIS-.*tif" TerrainFinal MasqImGlob=filtre.tif NbVI=2 ZoomF=2 ResolTerrain=0.5 DefCor=0 CostTrans=4 EZA=1</pre>
 Once the DEM is generated, other products can be generated :

@@ Ligne 159 : / Ligne 159 : @@
 To be able to ignore the fiducial marks and other inscriptions on the images that would yield nonsensical tie points, a mask need to be created.
-<pre>mm3d SaisieMasqQT ""OIS-image.tif"</pre>
+<pre>mm3d SaisieMasqQT "OIS-image.tif"</pre>
 Once created, the mask should be renamed <i>filtre.tif</i>.

@@ Ligne 277 : / Ligne 277 : @@
 Perform a bundle adjustment and a refinement of the camera calibration using the GCPs. The numerical values in the GCP option are the estimate of the quality of your GCPs (the first in meters in the world coordinate and the second in pixels in your input).
-<pre>mm3d Campari "OIS-.*tif" TerrainBrut TerrainFinal GCP=[GCPs.xml,5,MeasuresFinales-S2D.xml,2] AllFree=1</pre>
+<pre>mm3d Campari "OIS-.*tif" TerrainBrut TerrainFinal GCP=[GCPs.xml,5,MeasuresFinales-S2D.xml,2] SH=HomolMasqFiltered AllFree=1</pre>
 Visualize if wanted:

@@ Ligne 171 : / Ligne 171 : @@
 Then run the whole set with the calibration as input:
-<pre>mm3d Tapas RadialBasic "OIS.*tif" InOri=CalibInit Out=Relative SH=HomolMasqFiltered LibFoc=0</pre>
+<pre>mm3d Tapas RadialBasic "OIS.*tif" InCal=CalibInit Out=Relative SH=HomolMasqFiltered LibFoc=0</pre>
 <div class="toccolours mw-collapsible mw-collapsed" style="background-color: Lavender">

← Version précédente		Version du 16 février 2017 à 17:31
Ligne 6 :		Ligne 6 :

	If you are looking for a special area in France, you can use the Geoportail ([[IGN]]) to download your own images and process it.		If you are looking for a special area in France, you can use the Geoportail ([[IGN]]) to download your own images and process it.
		+
		+	This tutorial is mostly designed and maintained by [[Luc Girod]], if you have questions about it, please get in touch with him.

	=Download=		=Download=

@@ Ligne 235 : / Ligne 235 : @@
 </pre>
-<pre>mm3d GCPConvert AppInFile GCPs.txt ChgSys=wgs84@MySystem.xml</pre>
+<pre>mm3d GCPConvert AppInFile GCPs.txt ChgSys=DegreeWGS84@MySystem.xml</pre>
 You can also use "ChgSys=MySystem1.xml@MySystem2.xml" if your input is in a different specified system.