[[Image:picto-liste.png|25px|link=Tutorials]] [[Tutorials|Tutorials index]]
==Description==
The post by user mathias.f on the MicMac forum will be translated into forum format soon : [http://forum-micmac.forumprod.com/vault-or-dome-develop-developpe-d-une-vo-te-ou-d-un-dome-t1597.html Forum topic]

Unwrap Vault

2018-08-02T12:47:56Z

Luc Girod :

The post by user mathias.f on the MicMac forum will be translated into forum format soon : [http://forum-micmac.forumprod.com/vault-or-dome-develop-developpe-d-une-vo-te-ou-d-un-dome-t1597.html Forum topic]

[[Tutorials|Back to tutorials]]

2017-07-04T09:34:05Z

Luc Girod : /* With Setup (recommended) */

French version below.
==With Setup (recommended)==

*NEW! Download the binaries from the GitHub:<br>
<code>https://github.com/micmacIGN/micmac/releases</code><br>

*Download the setup from :<br>
<code>http://logiciels.ign.fr/?Telechargement,20</code><br>
Choose the 32 or 64bits version according to your Windows system ("Configuration pannel/System"). Then, follow the differents step of the setup.
*We recommand you to install MicMac in a directory where the path doesn't include any space or special character (i.e.: " ","&","@" etc...). The common error is to install it in "C:\Program Files (x86)" which include space.

==With binaries==
===Download & Install===
*Download the binaries from :<br>
<code>http://logiciels.ign.fr/?Telechargement,20</code><br>
Choose the 32 or 64bits version according your Windows system ("Configuration pannel/System").Then extract the ".zip" file where you want to install MicMac.
*We recommand you to install MicMac in a directory where the path doesn't include any space or special character (i.e.: " ","&","@" etc...). The common error is to install it in "C:\Program Files (x86)" which include space.

===Set environnemental path===
It is possible that you should have to indicate to your computer where the MicMac binaries are. You need to open the configuration pannel -> System -> Set parameters (Screenshot 1)
[[Image:Capture ecran1.png|thumb|250px||alt=Exemple d'image en hauteur|Screenshot 1]]

Then Advanced parameters -> Environnemental variable (Screenshot 2).
[[Image:Capture ecran2.png|thumb|250px||alt=Exemple d'image en hauteur|Screenshot 2]]

Double-click in the path variable and add at the end (Screenshot 3):<br>
<code>;installation_directory\micmac\bin (for example ";C:\micmac\bin")</code>
[[Image:Capture ecran3.jpg|thumb|250px||alt=Exemple d'image en hauteur|Screenshot 3]]

==Check Installation==
===List of micmac tools===
To check if MicMac is well installed, you can type in a terminal :
<pre>mm3d</pre>
This will return the list of all tools available for micmac library.<br>

===Check dependencies===
Sometimes, even when the MicMac installation is good, you can encounter messages in the terminal such as :
*"Warning Exiftool not installed"
This warnings means that MicMac can't find an external tools.
MicMac integrates a module which can check if all dependencies are correctly installed. It can be called by typing :<br>
<pre>
mm3d CheckDependencies
</pre>

[[Image:checkdependencies.png|thumb|250px||alt=Exemple d'image en hauteur|Screenshot 4]]

You will get an output like Screenshot 4, with different informations :
*MicMac revision : micmac source code version number
*Qt : "enabled" permit you to use graphical interfaces such SaisieAppuisQT
*Kakadu : Image library for JPEG2000 management
*MicMac directory : Directory where MicMac is installed
*make :
*exiftool : library of image metadata management
*exiv2 : library of image metadata management
*convert :
*proj : library for cartographic projection
*cs2cs :

==Option : Install Exiv2==
*You can download the exe file from :<br>
<code>http://www.exiv2.org/download.html</code><br>
click on "Download exiv2.exe"

*Move the "exiv2.exe" file in the MicMac binaries directory ("C:\micmac\bin")

==Known problems==

==How to use MicMac with Windows==

This video shows some basic steps on how to use MicMac on Windows.

{{#ev:youtube|https://www.youtube.com/watch?v=tezTFRanTYA|500x300px}}

=French Version : Installation MicMac pour windows=

==Avec le setup==

*Cette solution consiste à installer MicMac en téléchargeant le setup proposé sur le site de l’IGN. Attention : L’installation par setup est uniquement recommandée si la version d’installation est la plus récente disponible sur le site.

*Télécharger le setup depuis http://logiciels.ign.fr/?Micmac. Choisir la version 32 ou 64bits en fonction de son type de système. Suivre les étapes d’installation semblable aux setups des logiciels Windows.

==Avec les binaires==

*Télécharger les binaires précompilés depuis http://logiciels.ign.fr/?Micmac. Choisir la bonne version 32 ou 64 bit selon votre ordinateur. (Pour connaître la version du système de son ordinateur : Démarrer/Panneau de configuration/Système).

*Décompresser le fichier micmac_bin_windows_x86_64_rev3982.zip dans le répertoire voulu. Attention : proscrire les arborescences contenant des espaces ou/et des accents) !

*Aller dans « Panneau de configuration Système », cliquer sur « Modifier les paramètres »

*Cliquer ensuite sur « Paramètres Système Avancés » et enfin sur « Variables d’environnement. »

*Dans « Variables Système » choisir la variable « Path » et cliquer sur « Modifier ».

*Dans « Valeur de la variable », rajouter à la fin un point-virgule suivi du nom du répertoire contenant les binaires ;répertoire_des_binaires. Attention : proscrire les arborescences contenant des espaces ou/et des accents) !

*Quitter le panneau de configuration en cliquant sur OK à chaque fois.

*Lancer une invite de commande (Installation_MicMac_Windows) et taper mm3d pour vérifier que tout est bien installé. La liste des commandes disponibles doit s’afficher.

*Taper ensuite : mm3d CheckDependencies

*Vérifier que tout est installé en vérifiant que chaque paquet est « found ». Si un paquet est manquant, se référer au fichier LisezMoi téléchargé avec les binaires.

==Optionnel : Installation de Exiv2==
*Se rendre sur http://www.exiv2.org/download.html puis télécharger l’éxecutable Windows en cliquant sur « Download exiv2.exe ».

*Il faut ensuite placer le fichier exiv2.exe à côté des binaires de MicMac (par défaut C:\MicMac64bits\bin)

[[Category:Installation]]

AperiCloud

2017-06-29T14:05:15Z

Luc Girod : /* AperiCloud doesn't support 64 bits ply */

[[Image:picto-liste.png|25px]] [[Command|List of commands]]
==Description==
AperiCloud is used for generating a visualization of and sparse 3D model and cameras position, computed by [[Tapas]] for example.

===Syntax===
The global syntax for AperiCloud is :
<pre>mm3d AperiCloud FullName Orientation NamedArgs</pre>

===Results===
AperiCloud, generate a sparse 3D cloud (".ply") with the position and orientation of each camera. If you don't use "Out" option, the file will be stored in your working directory at the name : <i>AperiCloud-Ori_Name.ply</i>. You can use MeshLab to visualize the file.<br>
[[Image:01_Gravillonn_AO.jpg|x200px]]
[[Image:Pierrerue2.png|x200px]]
[[Image:Exercice_Fontaine1.png|x200px]]

===Workflow===
AperiCloud is a tool dedicated to visualization of orientations. So you can use it every times you process a orientation. When you have processed the final orientation, you can define a 3D mask with the ".ply" file generated by AperiCloud (SaisieMasqQT).<br>
[[Image:Picto-previous.png|20px]] Previous Command : [[Tapas]],[[GCPBascule]],[[Campari]]<br>
[[Image:Picto-next.png|20px]] Next Command : [[SaisieMasqQT]].<br>


===Help===
You can access to the help by typing :
<pre>mm3d AperiCloud -help</pre>

Mandatory unnamed args :
*string :: {Full name (Dir+Pattern)}
*string :: {Orientation directory}
Named args :
*[Name=ExpTxt] INT :: {Point in txt format ? (Def=false)}
*[Name=Out] string :: {Result (Def=AperiCloud.ply)}
*[Name=Bin] INT :: {Ply in binary mode (Def=true)}
*[Name=RGB] INT :: {Use RGB image to texturate points (Def=true)}
*[Name=SeuilEc] REAL :: {Max residual (Def=10)}
*[Name=LimBsH] REAL :: {Limit ratio base to height (Def=1e-2)}
*[Name=WithPoints] bool :: {Do we add point cloud? (Def=true) }
*[Name=CalPerIm] bool :: {If a calibration per image was used (Def=false)}
*[Name=Focs] Pt2dr :: {Interval of Focal}
*[Name=WithCam] bool :: {With Camera (Def=true)}
*[Name=ColCadre] Pt3di :: {Col of camera rect Def= 255 0 0 (Red)}
*[Name=ColRay] Pt3di :: {Col of camera rect Def= 0 255 0 (Green)}
*[Name=SH] string :: {Set of Hom, Def="", give MasqFiltered for result of HomolFilterMasq}

===Example===

In example, you can launch this command in the Mur Saint Martin dataset :
<pre>mm3d AperiCloud «.*JPG » Out=MEP </pre>

== Known limitations ==
=== AperiCloud doesn't support 64 bits ply ===
Use of AperiCloud with coordinates with many digits generates models which seeems sliced, or very sparse, because this tool doesn't support 64 bits ply models.

Due to this limitations, it is recommended to use local or truncated coordinate. Or not to use AperiCloud after georeferencement in a projected coordinates system.
Look to [http://forum-micmac.forumprod.com/point-cloud-seems-sliced-and-very-sparse-after-coord-chgt-t1270.html| this thread] for an example and further information.

For precision and accuracy reasons, it is anyways always recommended to work in a local system at least until just before the correlation step (after Campari), and use [[ChgSysCo]] to then get to the projected coordinate system.

2017-05-31T13:51:32Z

Luc Girod : /* Preparing the data */

==Description==
This tutorial is there to give the workflow necessary to compute MMASTER DEMs from ASTER L1A data. The method is described in details in an upcoming paper (will be linked here as soon as it is published).

==Requirement==
Part of the algorithms used in MMASTER rely on the ALGLIB library. Since it is licensed under the GNU GPL license, it cannot be part of the standard MicMac release (under CeCILL-B license). For that reason, it is necessary to compile MicMac from source, getting the IncludeALGLIB branch of the Git repository.

You also need to have Gdal installed on your system.

==Download==
Any AST_L1A scene can work here, you can download them from [https://earthdata.nasa.gov/]. When chosing the file format, choose GeoTiff (NOT hdf).

==Preparing the data==
Ecaples are given using a scenenamed ''AST_L1A_00302212004225011_20170204145459_24750'' and the UTM4N cartographic projection.

First extract the GeoTiff L1A to a folder named SceneName/RawData (ex AST_L1A_00302212004225011_20170204145459_24750/RawData).

Then from the folder where **SCENE_FOLDER** (ex AST_L1A_00302212004225011_20170204145459_24750/) is situated run :
*<pre>WorkFlowMMASTER.sh -s AST_L1A_00302212004225011_20170204145459_24750 -z "4 +north"</pre>
*extra options (given here with default values) : <pre>-t 30 -f 1 -n false -c 0.7</pre>

==Processing==
The following code can be in a bash (.sh) file and directly ran.
<pre>

#Fixed symboles
Nx="_3N.xml"
Bx="_3B.xml"
Nt="_3N.tif"
Bt="_3B.tif"
Bcor="_3B.tif_corrected.tif"
RPC="RPC_"
scene_set=0
utm_set=0
# add default values for ZoomF, RESTERR, CorThr and NoCorDEM
ZoomF=1
RESTERR=30
CorThr=0.7
NoCorDEM=false

while getopts "s:z:c:nf:t:h" opt; do
case $opt in
h)
echo "Run the second step in the MMASTER processing chain."
echo "usage: WorkflowMMASTER.sh -s SCENENAME -z 'UTMZONE' -f ZOOMF -t RASTERRES -h"
echo " -s SCENENAME: Aster scenename/folder where data is located."
echo " -z UTMZONE : UTM Zone of area of interest. Takes form 'NN +north(south)'"
echo " -c CorThr : Correlation Threshold for estimates of Z min and max (optional, default : 0.7)"
echo " -n NoCorDEM : Compute DEM with the uncorrected 3B image (computing with correction as well)"
echo " -f ZOOMF : Run with different final resolution (optional; default: 1)"
echo " -t RESTERR : Run with different terrain resolution (optional; default: 30)"
echo " -h : displays this message and exits."
echo " "
exit 0
;;
n)
NoCorDEM=$OPTARG
;;
s)
name=$OPTARG
scene_set=1
;;
z)
UTM=$OPTARG
utm_set=1
;;
c)
CorThr=$OPTARG
echo "CorThr set to $CorThr"
;;
f)
ZoomF=$OPTARG
;;
t)
RESTERR=$OPTARG
;;
\?)
echo "RunMicMacAster.sh: Invalid option: -$OPTARG" >&2
exit 1
;;
:)
echo " WorkflowMMASTER.sh: Option -$OPTARG requires an argument." >&2
exit 1
;;
esac
done

#Variable symboles
echo $name
echo $UTM
cd $name
pwd
cd RawData
pwd
mm3d SateLib ASTERGT2MM $name
cd ..

mm3d SateLib Aster2Grid $name$Bx 20 "+proj=utm +zone=$UTM +datum=WGS84 +units=m +no_defs" HMin=-500 HMax=9000 expDIMAP=1 expGrid=1
mm3d SateLib Aster2Grid $name$Nx 20 "+proj=utm +zone=$UTM +datum=WGS84 +units=m +no_defs" HMin=-500 HMax=9000 expDIMAP=1 expGrid=1
mm3d SateLib Aster2Grid "FalseColor_$name.xml" 20 "+proj=utm +zone=$UTM +datum=WGS84 +units=m +no_defs" HMin=-500 HMax=9000 expDIMAP=1 expGrid=1

mm3d Malt Ortho ".*$name(|_3N|_3B).tif" GRIBin ImMNT="$name(_3N|_3B).tif" MOri=GRID ZMoy=2500 ZInc=2500 ZoomF=8 ZoomI=32 ResolTerrain=30 NbVI=2 EZA=1 Regul=0.1 DefCor=$CorThr DoOrtho=0 DirMEC=MEC-Mini

gdalinfo -nomd -norat -noct -nofl -stats MEC-Mini/Z_Num6_DeZoom8_STD-MALT.tif > gdalinfo.txt
deminfo=$(grep -P 'Minimum+' gdalinfo.txt)
Min=$(echo $deminfo | cut -d, -f1 | tr -d ' ' | tr -d 'Minimum=' | xargs printf "%.0f")
Max=$(echo $deminfo | cut -d, -f2 | tr -d ' ' | tr -d 'Maximum=' | xargs printf "%.0f")

echo Min=$Min
echo Max=$Max

#Filter obvious error in min/max (limit to earth min/max)
Min=$((($Min)<-420 ? -420 : $Min))
Max=$((($Max)>8850 ? 8850 : $Max))
#next 2 lines is basically if the auto min/max function failed / DEM is really bad, happen if a lot of sea or a lot of clouds
Min=$((($Min)>8850 ? -420 : $Min))
Max=$((($Max)<-420 ? 8850 : $Max))
#From min/max, compute the nb of grids needed in Z and the values for ZMoy and Zinc
DE=$(echo $Max - $Min| bc )
NbLvl=$(echo $DE/200| bc )
NbLvl=$((($NbLvl)<10 ? 10 : $NbLvl))
Mean=$(echo $Max + $Min| bc )
Mean=$(echo $Mean/2| bc )
Inc=$(echo $Max - $Mean| bc | xargs printf "%.0f")
echo Min=$Min
echo Max=$Max
echo NbLvl=$NbLvl
echo Mean=$Mean
echo Inc=$Inc
echo Min Max NbLvl Mean Inc >> Stats.txt
echo $Min $Max $NbLvl $Mean $Inc >> Stats.txt

#Re compute RPCs with updated min/max
mm3d SateLib Aster2Grid $name$Bx $NbLvl "+proj=utm +zone=$UTM +datum=WGS84 +units=m +no_defs" HMin=$Min HMax=$Max expDIMAP=1 expGrid=1
mm3d SateLib Aster2Grid $name$Nx $NbLvl "+proj=utm +zone=$UTM +datum=WGS84 +units=m +no_defs" HMin=$Min HMax=$Max expDIMAP=1 expGrid=1
mm3d SateLib Aster2Grid "FalseColor_$name.xml" $NbLvl "+proj=utm +zone=$UTM +datum=WGS84 +units=m +no_defs" HMin=$Min HMax=$Max expDIMAP=1 expGrid=1

mm3d MMTestOrient $name$Bt $name$Nt GRIBin PB=1 MOri=GRID ZoomF=1 ZInc=$Inc ZMoy=$Mean

# if we want to compute the uncorrected DEM
if [ "$NoCorDEM" = true ]; then #check variable name!
mm3d Malt Ortho ".*$name(|_3N|_3B).tif" GRIBin ImMNT="$name(_3N|_3B).tif" ImOrtho="FalseColor_$name.tif" MOri=GRID ZInc=$Inc ZMoy=$Mean ZoomF=1 ZoomI=32 ResolTerrain=30 NbVI=2 EZA=1 DefCor=0 Regul=0.1 ResolOrtho=2 DirMEC=MEC-NoCor
fi

#Applying correction to the 3B image
mm3d SateLib ApplyParralaxCor $name$Bt GeoI-Px/Px2_Num16_DeZoom1_Geom-Im.tif FitASTER=1 ExportFitASTER=1
mkdir ImOrig
mv $name$Bt ImOrig/$name$Bt
mv $name$Bcor $name$Bt

# Correlation with corrected image
mm3d Malt Ortho ".*$name(|_3N|_3B).tif" GRIBin ImMNT="$name(_3N|_3B).tif" ImOrtho="FalseColor_$name.tif" MOri=GRID ZInc=$Inc ZMoy=$Mean ZoomF=$ZoomF ZoomI=32 ResolTerrain=$RESTERR NbVI=2 EZA=1 DefCor=0 Regul=0.1 ResolOrtho=2

mm3d GrShade MEC-Malt/Z_Num9_DeZoom1_STD-MALT.tif ModeOmbre=IgnE
mm3d Tawny Ortho-MEC-Malt/ RadiomEgal=0
mm3d Nuage2Ply MEC-Malt/NuageImProf_STD-MALT_Etape_9.xml Out=$name.ply Attr=Ortho-MEC-Malt/Orthophotomosaic.tif

#Postprocessing to have easy to use geotifs
mkdir -p OUTPUT
outdir=$(pwd)/OUTPUT
cd MEC-Malt
finalimgs=($(ls Z_Num*_DeZoom1_STD-MALT.tif))
finalmsks=($(ls AutoMask_STD-MALT_Num*.tif))
finalcors=($(ls Correl_STD-MALT_Num*.tif))
lastimg=${finalimgs[-1]}
lastmsk=${finalmsks[-1]}
lastcor=${finalcors[-1]}
# strip the extension
laststr="${lastimg%.*}"
maskstr="${lastmsk%.*}"
corrstr="${lastcor%.*}"

cp -v $laststr.tfw $maskstr.tfw
cp -v $laststr.tfw $corrstr.tfw
# now, assign the CRS we got to the mask, dem, and apply.
echo "Georeferencing correlation mask"
gdal_translate -a_nodata 0 -a_srs "+proj=utm +zone=$UTM +ellps=WGS84 +datum=WGS84 +units=m +no_defs" $lastcor $dir\_CORR.tif
echo "Creating temporary georeferenced DEM"
gdal_translate -a_srs "+proj=utm +zone=$UTM +ellps=WGS84 +datum=WGS84 +units=m +no_defs" $lastimg tmp_geo.tif
echo "Creating temporary georeferenced Mask"
gdal_translate -a_srs "+proj=utm +zone=$UTM +ellps=WGS84 +datum=WGS84 +units=m +no_defs" -a_nodata 0 $lastmsk tmp_msk.tif
cd ../
if [ -d "Ortho-MEC-Malt" ]; then
cd Ortho-MEC-Malt
echo "Creating double size correlation mask for ortho"
gdal_translate -tr 15 15 -a_srs "+proj=utm +zone=$UTM +ellps=WGS84 +datum=WGS84 +units=m +no_defs" -a_nodata 0 ../MEC-Malt/$lastmsk tmp_mskDouble.tif
echo "Creating temporary georeferenced ortho"
gdal_translate -a_srs "+proj=utm +zone=$UTM +ellps=WGS84 +datum=WGS84 +units=m +no_defs" Orthophotomosaic.tif tmp_V123.tif
cd ../
fi
cd MEC-Malt
# apply the mask
echo "Applying mask to georeferenced DEM"
gdal_calc.py -A tmp_msk.tif -B tmp_geo.tif --outfile=$dir\_Z.tif --calc="B*(A>0)" --NoDataValue=-9999
cp -v $dir\_Z.tif $outdir/$datestr #might be good to code orig. wd here.
gdaldem hillshade $dir\_Z.tif $outdir/$datestr/$dir\_HS.tif
gdal_calc.py -A $dir\_CORR.tif --outfile=$outdir/$datestr/$dir\_CORR.tif --calc="((A.astype(float)-127)/128)*100" --NoDataValue=-9999
#cp -v $dir\_CORR.tif $outdir/$datestr #might be good to code orig. wd here.
rm -v tmp_msk.tif tmp_geo.tif
rm -v $dir\_Z.tif $dir\_CORR.tif
cd ../
if [ -d "Ortho-MEC-Malt" ]; then
cd Ortho-MEC-Malt
gdal_calc.py -B tmp_mskDouble.tif -A tmp_V123.tif --outfile=$dir\_V123.tif --calc="((A!=255)*(A+1)+(A==255)*A)*(B>0)" --NoDataValue=0 --allBands=A
#Expression complicated to solve real 0 values not being NoData and 255 no being +1-ed to 0
rm -v tmp_V123.tif tmp_mskDouble.tif
cp -v $dir\_V123.tif $outdir/$datestr
rm -v $dir\_V123.tif
cd ../
fi
echo "MMASTER processing is complete."

</pre>