Difference between revisions of "MGI / Balkans coordinate systems"

m (Validation of CRS parameters)
m (Other coordinate systems of interest)
Line 290: Line 290:
 
To generate the grid cell id in the [http://www.ec-gis.org European Reference Grid System] ([http://spatialreference.org/ref/epsg/3035/ EPSG:3035]), we can run:
 
To generate the grid cell id in the [http://www.ec-gis.org European Reference Grid System] ([http://spatialreference.org/ref/epsg/3035/ EPSG:3035]), we can run:
  
<geshi lang=R lines=0>
+
<pre>
 
> library(maptools); library(rgdal)
 
> library(maptools); library(rgdal)
  
Line 324: Line 324:
 
> out <- paste(getwd(), "grid_croatia_etrs_laea_1k", sep="/")
 
> out <- paste(getwd(), "grid_croatia_etrs_laea_1k", sep="/")
 
> writeSpatialShape(HR1kmgrid.spoly, out)
 
> writeSpatialShape(HR1kmgrid.spoly, out)
</geshi>
+
</pre>
  
 
This will create a [http://spatial-analyst.net/DATA/grid_croatia_etrs_laea_1k.zip '''shape file'''] with grid nodes placed exactly on the reference grid. Each grid cell has an unique name, which is abbreviation of the XY coordinates of the origin of the grid. E.g. <tt>4572_2625</tt> corresponds to the grid node <tt>X=4572500</tt>, <tt>Y=2625500</tt>.
 
This will create a [http://spatial-analyst.net/DATA/grid_croatia_etrs_laea_1k.zip '''shape file'''] with grid nodes placed exactly on the reference grid. Each grid cell has an unique name, which is abbreviation of the XY coordinates of the origin of the grid. E.g. <tt>4572_2625</tt> corresponds to the grid node <tt>X=4572500</tt>, <tt>Y=2625500</tt>.

Revision as of 23:01, 18 November 2018