radproc.arcgis.create_idarray

radproc.arcgis.create_idarray(projectionFile, idRasterGermany, clipFeature, idRaster, extendedNationalGrid=True)

Creates an ID Array for a study area.

Creates a new ID raster for Germany, clips it to study area and converts the raster to a one-dimensional numpy array.

Parameters:
projectionFile : string
Path to a file containing stereographic projection definition. File type may be Feature Class, Shapefile, prj file or grid.
idRasGermany : string
Path and name for the output ID raster of Germany to be created.
clipFeature : string
Path to the clip feature defining the extent of the study area. File type may be Shapefile or Feature Class. If clipFeature == None, the nationwide ID raster will not be clipped.
idRas : string
Path and name for the output ID raster of the study area to be created.
extendedNationalGrid : bool (optional, default: True)
True: extended 900 x 1100 national RADOLAN grid, False: 900x900 national grid.
Returns:
idArr : one-dimensional numpy array
containing ID values of dtype int32
Note:

Note

The RADOLAN data are provided in a custom stereographic projection defined by the DWD. As there is no transformation method available yet, it is not possible to directly perform any geoprocessing tasks with RADOLAN and geodata with other spatial references. Nevertheless, ArcGIS is able to perform a correct on-the-fly transformation to display the data together. The clip function uses this as a work-around solution to “push” the clip feature into the RADOLAN projection. Hence, the function works with geodata in different projections, but the locations of the cells might be slightly inaccurate.

To use the custom RADOLAN projection as projectionFile for the output ID raster, you can specify the prj-file provided in radproc.sampledata:

>>> from radproc.sampledata import get_projection_file_path
>>> projectionFile = get_projection_file_path()