Add is_lmi/ for Iceland

12 files changed, 581 insertions, 1 deletions

diff --git a/CMakeLists.txt b/CMakeLists.txt
index 3b5768b..415350e 100644
--- a/CMakeLists.txt
+++ b/CMakeLists.txt
@@ -32,6 +32,7 @@ set(CPACK_SOURCE_INSTALLED_DIRECTORIES
     "${CMAKE_CURRENT_SOURCE_DIR}/dk_sdfe/" "."
     "${CMAKE_CURRENT_SOURCE_DIR}/es_cat_icgc/" "."
     "${CMAKE_CURRENT_SOURCE_DIR}/fr_ign/" "."
+    "${CMAKE_CURRENT_SOURCE_DIR}/is_lmi/" "."
     "${CMAKE_CURRENT_SOURCE_DIR}/nz_linz/" "."
     "${CMAKE_CURRENT_SOURCE_DIR}/pt_dgt/" "."
     "${CMAKE_CURRENT_SOURCE_DIR}/se_lantmateriet/" "."
diff --git a/README.DATUMGRID b/README.DATUMGRID
index 05a72c0..9af9569 100644
--- a/README.DATUMGRID
+++ b/README.DATUMGRID
@@ -55,6 +55,7 @@ For example us_noaa for files originating from US NOAA.
 * [dk_sdfe: Denmark SDFE](https://github.com/OSGeo/proj-datumgrid-geotiff/blob/master/dk_sdfe/)
 * [es_cat_icgc: Spain Catalonia ICGC](https://github.com/OSGeo/proj-datumgrid-geotiff/blob/master/es_cat_icgc/)
 * [fr_ign: France IGN](https://github.com/OSGeo/proj-datumgrid-geotiff/blob/master/fr_ign/)
+* [is_lmi: Iceland LMI](https://github.com/OSGeo/proj-datumgrid-geotiff/blob/master/is_lmi/)
 * [nz_linz: New Zealand LINZ](https://github.com/OSGeo/proj-datumgrid-geotiff/blob/master/nz_linz/)
 * [pt_dgt: Portugal DGT](https://github.com/OSGeo/proj-datumgrid-geotiff/blob/master/pt_dgt/)
 * [se_lantmateriet: Sweden Lantmateriet](https://github.com/OSGeo/proj-datumgrid-geotiff/blob/master/se_lantmateriet/)
diff --git a/index.html b/index.html
index f50dbd8..6386390 100644
--- a/index.html
+++ b/index.html
@@ -154,6 +154,14 @@ Those files are designed to be mainly used by <a href="https://proj.org">PROJ</a
 <li><a href="gr3df97a.tif">gr3df97a.tif</a>: Geocentric translation from NTF (IGNF:NTF) to RGF93 (EPSG:4964). Last modified: 2019/12/17</li>
 <li><a href="ntf_r93.tif">ntf_r93.tif</a>: NTF (EPSG:4275) to RGF93 (EPSG:4171). Last modified: 2019/12/16</li>
 </ul><ul>
+<li><a href="README.is_lmi">README.is_lmi</a></li>
+<li><a href="ISL">ISL</a>. Last modified: 2019/12/18</li>
+<li><a href="ISN2004_ISN2016.tif">ISN2004_ISN2016.tif</a>: ISN2004 (EPSG:5324) to ISN2016 (EPSG:8086). Last modified: 2019/12/18</li>
+<li><a href="ISN93_ISN2016.tif">ISN93_ISN2016.tif</a>: ISN93 (EPSG:4659) to ISN2016 (EPSG:8086). Last modified: 2019/12/18</li>
+<li><a href="Icegeoid_ISN2004.tif">Icegeoid_ISN2004.tif</a>: ISN2004 (EPSG:5323) to ISH2004 height (EPSG:8089). Last modified: 2019/12/18</li>
+<li><a href="Icegeoid_ISN2016.tif">Icegeoid_ISN2016.tif</a>: ISN2016 (EPSG:8085) to ISH2004 height (EPSG:8089). Last modified: 2019/12/18</li>
+<li><a href="Icegeoid_ISN93.tif">Icegeoid_ISN93.tif</a>: ISN93 (EPSG:4945) to ISH2004 height (EPSG:8089). Last modified: 2019/12/18</li>
+</ul><ul>
 <li><a href="README.nz_linz">README.nz_linz</a></li>
 <li><a href="nzgd2kgrid0005.tif">nzgd2kgrid0005.tif</a>: NZGD49 (EPSG:4272) to NZGD2000 (EPSG:4167). Last modified: 2019/12/16</li>
 </ul><ul>
@@ -280,7 +288,7 @@ Those files are designed to be mainly used by <a href="https://proj.org">PROJ</a
 <li><a href="stpaul.tif">stpaul.tif</a>: NAD27 (EPSG:4267) to NAD83 (EPSG:4269). Last modified: 2019/12/16</li>
 </ul>
 <p>
-Total size of content: 463 MB
+Total size of content: 466 MB
 </p>
 <div style="text-align: center;">
 <a href="https://aws.amazon.com/what-is-cloud-computing">
diff --git a/is_lmi/.github/README.md b/is_lmi/.github/README.md
new file mode 120000
index 0000000..f4452fd
--- /dev/null
+++ b/is_lmi/.github/README.md
@@ -0,0 +1 @@
+../README.is_lmi
\ No newline at end of file
diff --git a/is_lmi/ISL b/is_lmi/ISL
new file mode 100644
index 0000000..abcefd6
--- /dev/null
+++ b/is_lmi/ISL
@@ -0,0 +1,470 @@
+#------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+#
+#    PROJ implementations of current and legacy Icelandic coordinate systems
+#
+#------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+#
+#  It is the official definition and implementation of transformations between different Iceland's systems.
+#  
+#  
+#
+#  Revision history:
+#  -----------------
+#
+#  2019-11-30  Initial setup.
+#              Dalia Prizginiene <dalia@lmi.is>
+#
+#  
+#-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+#
+<metadata> +version=1.0.0 +origin=LMI +lastupdate=2019-11-30
+#
+#--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+#  Vertical datum 
+#--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+
+<ISH2004> proj=vgridshift grids=Icegeoid_ISN2004.gtx  # current, geoid height calculated using ISN2004 datum
+#--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+<ISN93_ISH2004> proj=vgridshift grids=Icegeoid_ISN93.gtx  # Icegeoid_ISN2004.gtx plus corrections for vertical difference between ISN93 and ISN2004.
+<ISN2016_ISH2004> proj=vgridshift grids=Icegeoid_ISN2016.gtx  # Icegeoid_ISN2004.gtx plus corrections for vertical difference between ISN2016 and ISN2004.
+
+#--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+#  Deformation grids 
+#--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+<ISN93_16h> proj=hgridshift grids=ISN93_ISN2016.gsb
+<ISN93_16v> proj=vgridshift grids=ISN93_ISN2016.gtx 
+
+<ISN04_16h> proj=hgridshift grids=ISN2004_ISN2016.gsb
+<ISN04_16v> proj=vgridshift grids=ISN2004_ISN2016.gtx
+
+<D_93_16> 
+    proj=pipeline 
+    step init=ISL:ISN93_16h   
+    step inv init=ISL:ISN93_16v 
+    
+<D_04_16> 
+    proj=pipeline 
+    step init=ISL:ISN04_16h   
+    step inv init=ISL:ISN04_16v
+
+#----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+#  Velocity model
+#----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+#
+# xy_grids=ISN_vel_beta.ct2 
+# z_grids=ISN_vel_beta.gtx
+#
+#----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+
+
+#----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+# LCC ISN2016 (epoch 2016.5, IGS2014, GRS80):  The primary system  
+#----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+<ISN2016_LAM>  proj=lcc  lat_1=64.25 lat_2=65.75 lat_0=65 lon_0=-19 x_0=2700000 y_0=300000  ellps=GRS80
+
+#----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+<gie>
+operation   init=ISL:ISN2016_LAM
+tolerance   0.1 mm
+accept      -17 66
+expect      2790800.4329    412932.6164 
+</gie>
+
+
+#-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+# LCC ISN2004 (epoch 2004.6, IGb00, GRS80):  The other systems 
+#-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+<ISN2004_LAM>   proj=lcc  lat_1=64.25 lat_2=65.75 lat_0=65 lon_0=-19 x_0=1700000 y_0=300000  ellps=GRS80
+
+#-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+<gie>
+operation   init=ISL:ISN2004_LAM
+tolerance   0.1 mm
+accept      -17 66
+expect      1790800.4329    412932.6164
+</gie>
+
+
+#------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+# LCC ISN93 (epoch 1993.6, ITRF93, GRS80):  The other systems 
+#------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+<ISN93_LAM>     proj=lcc  lat_1=64.25 lat_2=65.75 lat_0=65 lon_0=-19 x_0=500000  y_0=500000  ellps=GRS80
+
+#------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+<gie>
+operation   init=ISL:ISN93_LAM
+tolerance   0.1 mm
+accept      -17 66
+expect      590800.4329    612932.6164
+</gie>
+
+
+#-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+# LCC Hjorsey1955 (ellipsoid: International 1924):  The historical systems 2D
+#-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+<HJORSEY_LAM> 
+    proj=pipeline
+    step proj=lcc   lat_1=65  lat_0=65 lon_0=-18 x_0=-500000 y_0=500000 ellps=intl
+    step proj=axisswap  order=-1,2,3,4 
+
+#-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+<gie>
+operation   init=ISL:HJORSEY_LAM
+tolerance   0.1 mm
+accept      -17 66
+expect      454587.8410    611870.0974
+</gie>
+
+#--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+# LCC Reykjavik1900 (ellipsoid: Danish):  The historical systems 2D
+#--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+<RJ_LAM> proj=lcc lat_1=65  lat_0=65 lon_0=-19.022125 x_0=0 y_0=0 ellps=danish
+
+#--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+<gie>
+operation   init=ISL:RJ_LAM
+tolerance   0.1 mm
+accept      -17 66
+expect      91794.9450    112953.3882
+</gie>
+
+#---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+#  LCC for Kinematic datum shifting/ dynamic transformation
+#---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+<KIN_LAM>  proj=lcc  lat_1=64.25 lat_2=65.75 lat_0=65 lon_0=-19 x_0=4700000 y_0=300000  ellps=GRS80
+#---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+<gie>
+operation   init=ISL:KIN_LAM
+tolerance   3 cm
+accept      -17 66 0 2017.5
+expect      4790800.4329    412932.6164
+</gie>
+
+#----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+#                                      Helmert transformations
+#----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+# Helmert transformations are used to transform between Reykjavik1900 2D and ISN93 2D also between Hjorsey1955 2D and ISN93 2D
+#----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+<REYK19_ISN93>    proj=helmert  convention=coordinate_frame x=556.620 y=-167.901 z=945.664 rx=4.154 ry=-0.269 rz=-2.279 s=-3.729 
+<HJ_ISN93>        proj=helmert  x=-72.4 y=46.8 z=-82.7
+<REYK19_HJ>       proj=helmert  convention=coordinate_frame x=629.020 y=-214.701 z=1028.364 rx=4.154 ry=-0.269 rz=-2.279 s=-3.729 
+#
+#
+#-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+#                                      Actual transformations
+#-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+#
+#-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+# Reykjavik1900 (Geodetic) 2D -> ISN93 (LCC) 2D
+#-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+# from west of Greenwich
+<RJ_G> 
+    proj=pipeline
+    step proj=push +v_3
+    step proj=cart   ellps=danish  
+    step init=ISL:REYK19_ISN93   
+    step proj=cart  inv ellps=GRS80 
+    step proj=pop +v_3
+    step init=ISL:ISN93_LAM
+
+# Example:
+# echo  -21.930875   64.1421888889 0 0  | cct  init=ISL:RJ_G   
+
+# from west of Copenhagen
+<RJ_W> 
+    proj=pipeline 
+    step inv proj=longlat  pm=copenhagen
+    step proj=push +v_3
+    step proj=cart   ellps=danish  
+    step init=ISL:REYK19_ISN93 
+    step proj=cart  inv ellps=GRS80 
+    step proj=pop +v_3
+    step init=ISL:ISN93_LAM
+
+# Example:
+# echo  -34.50875   64.1421888889 0 0  | cct  init=ISL:RJ_W   
+#------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+<gie>
+operation   init=ISL:RJ_G
+tolerance   50 m
+accept      -21.930875   64.1421888889 0
+expect      357557.5871   407700.8896  0
+</gie>
+#------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+<gie>
+operation   init=ISL:RJ_W
+tolerance   50 m
+accept      -34.50875   64.1421888889  0
+expect      357557.5871   407700.8896  0
+</gie>
+#-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+# Hjorsey1955 (Geodetic) 2D -> ISN93 (LCC) 2D
+#-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+                
+<HJORSEY_ISN93> 
+    proj=pipeline  
+    step proj=push +v_3
+    step proj=cart  ellps=intl 
+    step init=ISL:HJ_ISN93 
+    step proj=cart inv ellps=GRS80 
+    step proj=pop +v_3
+    step init=ISL:ISN93_LAM  
+                  
+#-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+<gie>
+operation   init=ISL:HJORSEY_ISN93
+tolerance   15 m
+accept      -21.9933555556   64.1533913888889 0
+expect      354399.1771   409035.1420  0
+</gie>
+#-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+# ISN93 (Geodetic) 3D -> ISN2004 (LCC) 3D
+#-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+
+<ISN93_ISN04> 
+    proj=pipeline  
+    step init=ISL:D_93_16   
+    step inv init=ISL:D_04_16 
+    step init=ISL:ISN2004_LAM  
+
+# Example:
+# echo -17.981614272222  66.5274148361111 75.737 0 | cct init=ISL:ISN93_ISN04
+#--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+<gie>
+operation   init=ISL:ISN93_ISN04
+tolerance   5 cm
+accept      -17.981614272222  66.5274148361111 75.737 0
+expect      1745293.041  470683.539 75.653
+</gie>
+
+#-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+# ISN93 (LCC) 3D -> ISN2004 (LCC) 3D
+#-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+
+<ISN93_L_ISN04> 
+    proj=pipeline 
+    step inv init=ISL:ISN93_LAM 
+    step init=ISL:D_93_16   
+    step inv init=ISL:D_04_16 
+    step init=ISL:ISN2004_LAM  
+
+# Example:
+# echo 545293.098 670683.324   75.737 0 | cct init=ISL:ISN93_L_ISN04
+#--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+<gie>
+operation   init=ISL:ISN93_L_ISN04
+tolerance   5 cm
+accept      545293.098 670683.324   75.737 0
+expect      1745293.0413 470683.5387 75.6532
+</gie>
+
+#-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+# ISN93 (Geodetic) 3D -> ISN2004 (Geodetic + Vertical datum)
+#-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+
+<ISN93_ISN04_v> 
+    proj=pipeline  
+    step init=ISL:D_93_16   
+    step inv init=ISL:D_04_16 
+    step init=ISL:ISH2004 
+
+# Example:
+# echo -17.981614272222  66.5274148361111 75.737 0   | cct  init=ISL:ISN93_ISN04_v
+#--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+<gie>
+operation   init=ISL:ISN93_ISN04_v
+tolerance   5 cm
+accept      -17.981614272222  66.5274148361111 75.737 0
+expect      -17.9816154690  66.5274167690  11.6147
+</gie>
+
+#--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+# ISN2004 (Geodetic) 3D -> ISN2016 (LCC) 3D
+#--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+
+<ISN04_ISN16> 
+    proj=pipeline   
+    step init=ISL:D_04_16  
+    step init=ISL:ISN2016_LAM
+ 
+#--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+<gie>
+operation   init=ISL:ISN04_ISN16
+tolerance   5 cm
+accept      -17.98161545  66.5274167694444 75.653 0 
+expect      2745292.997 470683.7597 75.618
+</gie>
+
+#--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+# ISN2004 (LCC) 3D -> ISN2016 (LCC) 3D
+#--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+
+<ISN04_L_ISN16>
+    proj=pipeline 
+    step inv  init=ISL:ISN2004_LAM  
+    step init=ISL:D_04_16  
+    step init=ISL:ISN2016_LAM
+ 
+#--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+<gie>
+operation   init=ISL:ISN04_L_ISN16
+tolerance   5 cm
+accept      1745293.042 470683.539 75.653 0 
+expect      2745292.9970 470683.7595 75.6182
+</gie>
+
+#--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+# ISN2004 (Geodetic) 3D -> ISN2016 (Geodetic + Vertical datum)
+#--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+
+<ISN04_ISN16_v>
+    proj=pipeline   
+    step init=ISL:ISN04_16h   
+    step init=ISL:ISH2004 
+
+# Example:
+# echo -17.98161545  66.5274167694444 75.653 0   | cct  init=ISL:ISN04_ISN16_v 
+#--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+<gie>
+operation   init=ISL:ISN04_ISN16_v
+tolerance   5 cm
+accept      -17.98161545  66.5274167694444 75.653 0 
+expect      -17.9816163818 66.5274187527  11.6145
+</gie>
+
+#-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+# ISN93 (Geodetic) 3D -> ISN2016 (LCC) 3D
+#-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+
+<ISN93_ISN16> 
+    proj=pipeline   
+    step init=ISL:D_93_16  
+    step init=ISL:ISN2016_LAM  
+ 
+#--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+<gie>
+operation   init=ISL:ISN93_ISN16
+tolerance   5 cm
+accept      -17.981614272222  66.5274148361111 75.737 0
+expect      2745292.996  470683.760 75.618
+</gie>
+
+#-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+# ISN93 (LCC) 3D -> ISN2016 (LCC) 3D
+#-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+
+<ISN93_L_ISN16> 
+    proj=pipeline 
+    step inv   init=ISL:ISN93_LAM  
+    step init=ISL:D_93_16  
+    step init=ISL:ISN2016_LAM  
+ 
+#--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+<gie>
+operation   init=ISL:ISN93_L_ISN16
+tolerance   5 cm
+accept      545293.098 670683.324   75.737 0
+expect      2745292.996 470683.759 75.618 
+</gie>
+
+#-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+# ISN93 (Geodetic) 3D -> ISN2016 (Geodetic + Vertical datum)
+#-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+
+<ISN93_ISN16_v> 
+    proj=pipeline   
+    step init=ISL:D_93_16  
+    step init=ISL:ISN04_16v 
+    step init=ISL:ISH2004 
+
+# Example:
+# echo -17.981614272222  66.5274148361111 75.737 0   | cct  init=ISL:ISN93_ISN16_v 
+#--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+<gie>
+operation   init=ISL:ISN93_ISN16_v
+tolerance   5 cm
+accept      -17.981614272222  66.5274148361111 75.737 0
+expect      -17.9816164008 66.5274187522  11.6147
+</gie>
+
+#-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+# LCC ISN93 3D -> ISN93 (Geodetic + Vertical datum)
+#-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+
+<ISN93_ISH2004_v> 
+    proj=pipeline   
+    step inv init=ISL:ISN93_LAM   
+    step init=ISL:ISN93_ISH2004 
+
+# Example:
+# echo 545293.098 670683.324   75.737  0   | cct  init=ISL:ISN93_ISH2004_v 
+#--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+<gie>
+operation   init=ISL:ISN93_ISH2004_v
+tolerance   5 cm
+accept      545293.098 670683.324   75.737  0
+expect      -17.9816142725   66.5274148323   11.6147
+</gie>
+
+#-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+# LCC ISN2016 3D -> ISN2016 (Geodetic + Vertical datum)
+#-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+
+<ISN2016_ISH2004_v> 
+    proj=pipeline    
+    step inv init=ISL:ISN2016_LAM    
+    step init=ISL:ISN2016_ISH2004 
+
+# Example:
+# echo 2745292.998 470683.759  75.618   0   | cct  init=ISL:ISN2016_ISH2004_v 
+#--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+<gie>
+operation   init=ISL:ISN2016_ISH2004_v
+tolerance   5 cm
+accept      2745292.998 470683.759  75.618   0
+expect      -17.9816163630   66.5274187459   11.6143
+</gie>
+
+#---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+# ISN_DRF (Geodetic): epoch 2018.679 -> ISN_DRF (LCC): epoch  2017.5  3D
+#---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+
+<DRF_E>  
+    proj=pipeline  
+    step proj=cart  ellps=GRS80 
+    step proj=deformation  dt=-1.179 xy_grids=ISN_vel_beta.ct2 z_grids=ISN_vel_beta.gtx  
+    step proj=cart inv ellps=GRS80 
+    step init=ISL:KIN_LAM 
+
+# Example:
+# echo -18.1224832715  65.6854274786 134.2043   | cct init=ISL:DRF_E
+#---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+
+<gie>
+operation   init=ISL:DRF_E
+tolerance   3 cm
+accept      -18.1224832715  65.6854274786 134.2043
+expect      4740331.7212 376698.7742 134.2005
+</gie>            
+
+#----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+# ISN_DRF (Geodetic): epoch 2018.679 -> ISN2016 (LCC)    3D
+#----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+
+<DRF_2016>  
+    proj=pipeline 
+    step proj=cart ellps=GRS80 
+    step proj=deformation dt=-2.179 xy_grids=ISN_vel_beta.ct2 z_grids=ISN_vel_beta.gtx 
+    step proj=cart inv ellps=GRS80 
+    step init=ISL:ISN2016_LAM
+
+# Example:
+# echo -18.1224832715  65.6854274786 134.2043    | cct init=ISL:DRF_2016
+#----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+
+<gie>
+operation   init=ISL:DRF_2016
+tolerance   3 cm
+accept      -18.1224832715  65.6854274786 134.2043 
+expect      2740331.7315 376698.7505 134.1972
+</gie>   
diff --git a/is_lmi/ISN2004_ISN2016.tif b/is_lmi/ISN2004_ISN2016.tif
new file mode 100644
index 0000000..3ff07f6
--- /dev/null
+++ b/is_lmi/ISN2004_ISN2016.tif
diff --git a/is_lmi/ISN93_ISN2016.tif b/is_lmi/ISN93_ISN2016.tif
new file mode 100644
index 0000000..48afcd1
--- /dev/null
+++ b/is_lmi/ISN93_ISN2016.tif
diff --git a/is_lmi/Icegeoid_ISN2004.tif b/is_lmi/Icegeoid_ISN2004.tif
new file mode 100644
index 0000000..fad3c6f
--- /dev/null
+++ b/is_lmi/Icegeoid_ISN2004.tif
diff --git a/is_lmi/Icegeoid_ISN2016.tif b/is_lmi/Icegeoid_ISN2016.tif
new file mode 100644
index 0000000..c24075e
--- /dev/null
+++ b/is_lmi/Icegeoid_ISN2016.tif
diff --git a/is_lmi/Icegeoid_ISN93.tif b/is_lmi/Icegeoid_ISN93.tif
new file mode 100644
index 0000000..b915d2c
--- /dev/null
+++ b/is_lmi/Icegeoid_ISN93.tif
diff --git a/is_lmi/README.is_lmi b/is_lmi/README.is_lmi
new file mode 100644
index 0000000..68dfb98
--- /dev/null
+++ b/is_lmi/README.is_lmi
@@ -0,0 +1,97 @@
+# README.is_lmi
+
+The files in this section result from the conversion of datasets originating
+from [National Land Survey of Iceland](https://atlas.lmi.is/LmiData/index.php?id=626468364600)
+
+## Included grids
+
+### Iceland: ISN93 -> ISN2016
+
+2D geodetic transformation from Iceland Datum ISN93 to Iceland Datum ISN2016
+
+*Source*: [National Land Survey of Iceland](https://atlas.lmi.is/LmiData/index.php?id=626468364600)  
+*Format*: NTv2  
+*License*: [Creative Commons Attribution 4.0](https://creativecommons.org/licenses/by/4.0/)  
+*Credit*: (c)2017-2019 Grid created by Guðmundur Valsson.
+
+The transformation is designed so that the ISN2016 points should fit within 1-2 cm accuracy, but in the area around the earthquake area in the South and North of Vatnajökull the accuracy is worse.
+
+Recommended to transform from ISN93 2D geodetic (EPSG:4659) to ISN2016 2D geodetic (EPSG:8086), also this transformation has code in EPSG Dataset (EPSG:9232)
+
+* ISN93_ISN2016.tif
+
+### Iceland: ISN2004 -> ISN2016
+
+2D geodetic transformation from Iceland Datum ISN2004 to Iceland Datum ISN2016
+
+*Source*: [National Land Survey of Iceland](https://atlas.lmi.is/LmiData/index.php?id=626468364600)  
+*Format*: NTv2  
+*License*: [Creative Commons Attribution 4.0](https://creativecommons.org/licenses/by/4.0/)  
+*Credit*: (c)2017-2019 Grid created by Guðmundur Valsson.
+
+The transformation is designed so that the ISN2016 points should fit within 1-2 cm accuracy, but in the area around the earthquake area in the South and North of Vatnajökull the accuracy is worse.
+
+Recommended to transform from ISN2004 2D geodetic (EPSG:5324) to ISN2016 2D geodetic (EPSG:8086), also this transformation has code in EPSG Dataset (EPSG:9233)
+
+* ISN2004_ISN2016.tif
+
+### Iceland: ISN2004 ellipsoidal height  -> ISH2004 vertical datum
+
+*Source*: [National Land Survey of Iceland](https://atlas.lmi.is/LmiData/index.php?id=626468364600)  
+*Format*: GTX  
+*License*: [Creative Commons Attribution 4.0](https://creativecommons.org/licenses/by/4.0/)  
+*Credit*: (c)2017-2019 Grid created by Guðmundur Valsson.
+
+Geoid created in cooperation between the National Land Survey of Iceland and DTU Space in  December 2010. 
+A gravimetric geoid was fitted to National Height System ISH2004 (EPSG:8089) with over 300 GNSS levelling points
+with ellipsoidal height in ISN2004 (EPSG:5323). The std. of the fit was 1.7 cm.
+
+Grid to convert from ISN2004 ellipsoidal height (EPSG:5323) to national vertical datum ISH2004 (EPSG:8089)
+The horizontal grid coordinates are referenced to ISN2004
+
+See the `ISL` init-file for examples of use.
+
+* Icegeoid_ISN2004.tif 
+
+### Iceland: ISN93 ellipsoidal height  -> ISH2004 vertical datum
+
+*Source*: [National Land Survey of Iceland](https://atlas.lmi.is/LmiData/index.php?id=626468364600)  
+*Format*: GTX  
+*License*: [Creative Commons Attribution 4.0](https://creativecommons.org/licenses/by/4.0/)  
+*Credit*: (c)2017-2019 Grid created by Guðmundur Valsson.
+
+Icegeoid_ISN2004.tif plus corrections for vertical difference between ISN93(EPSG:4945) and ISN2004 (EPSG:5323).
+
+Grid to convert from ISN93 ellipsoidal height (EPSG:4945) to national vertical datum ISH2004 (EPSG:8089)
+The horizontal grid coordinates are referenced to ISN93
+
+See the `ISL` init-file for examples of use.
+
+* Icegeoid_ISN93.tif
+
+### Iceland: ISN2016 ellipsoidal height  -> ISH2004 vertical datum
+
+*Source*: [National Land Survey of Iceland](https://atlas.lmi.is/LmiData/index.php?id=626468364600)  
+*Format*: GTX  
+*License*: [Creative Commons Attribution 4.0](https://creativecommons.org/licenses/by/4.0/)  
+*Credit*: (c)2017-2019 Grid created by Guðmundur Valsson.
+
+Icegeoid_ISN2004.tif plus corrections for vertical difference between ISN2016(EPSG:8085) and ISN2004 (EPSG:5323).
+
+Grid to convert from ISN2016 ellipsoidal height (EPSG:8085) to national vertical datum ISH2004 (EPSG:8089)
+The horizontal grid coordinates are referenced to ISN2016
+
+See the `ISL` init-file for examples of use.
+
+* Icegeoid_ISN2016.tif
+
+## Included init-files
+
+### Iceland: System definitions
+
+*Source*: [National Land Survey of Iceland](https://atlas.lmi.is/LmiData/index.php?id=626468364600)  
+*License*: [Creative Commons Attribution 4.0](https://creativecommons.org/licenses/by/4.0/)  
+
+Init file with various transformations and conversions of relevance in Iceland.
+
+* ISL
diff --git a/regenerate_index_html.py b/regenerate_index_html.py
index ebc3aa7..340f5a0 100644
--- a/regenerate_index_html.py
+++ b/regenerate_index_html.py
@@ -13,6 +13,8 @@ for dirname in glob.glob('*'):
 total_size = 0
 set_files = set()
 for dirname in sorted(dirnames):
+    if '_' not in dirname:
+        continue
     filenames = []
     readme_filename = None
     for f in glob.glob(dirname + '/*'):