def _get_bundle_from_path(path):
'''
This method tries to instance every bundle using the given directory. When
a bundle is able to identify the files present in the directory the instance
of that bundle is returned to the caller.
'''
return get_bundle_from_path(path, '../../../mapper', BUNDLE_PACKAGE)
def _get_bundle_from_path(path):
"""
This method tries to instance every bundle using the given directory. When
a bundle is able to identify the files present in the directory the instance
of that bundle is returned to the caller.
"""
return get_bundle_from_path(path, "../../../mapper", BUNDLE_PACKAGE)
def test_rasterize_vector_and_write_to_disk(self):
from madmex.mapper.data.raster import new_options_for_create_raster_from_reference, create_raster_tiff_from_reference
from madmex.core.controller.commands import get_bundle_from_path
from madmex.mapper.data import raster
BUNDLE_PACKAGE = 'madmex.mapper.bundle'
landmask_path = '/Users/erickpalacios/Documents/CONABIO/MADMEXdata/eodata/footprints/country_mexico/'
bundle = get_bundle_from_path(landmask_path, '../../../mapper',
BUNDLE_PACKAGE)
extents_dictionary = {
u'x_range': 7521.0,
u'y_range': 7741.0,
u'properties': {
u'projection':
'PROJCS["UTM Zone 15, Northern Hemisphere",GEOGCS["Unknown datum based upon the WGS 84 ellipsoid",DATUM["Not specified (based on WGS 84 spheroid)",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",-93],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAMETER["false_northing",0],UNIT["Meter",1]]',
u'geotransform': (523185.0, 30.0, 0.0, 2033715.0, 0.0, -30.0)
},
u'x_offset': 'array([ 10., 30., 0., 50.])',
u'y_offset': 'array([-0., -0., -0., -0.])'
}
options_to_create = new_options_for_create_raster_from_reference(
extents_dictionary, raster.DATA_SHAPE,
(int(extents_dictionary['x_range']),
int(extents_dictionary['y_range']), 1), {})
image = '/Users/erickpalacios/Documents/CONABIO/Tareas/Redisenio_MADMEX/clasificacion_landsat/landsat8/classification/rasterize2.tif'
#The next line create a tiff with empty array
data = create_raster_tiff_from_reference(extents_dictionary, image,
None, options_to_create)
bundle.rasterize(data, [1], [1])
new_options_for_create_raster_from_reference(extents_dictionary,
raster.DATASET, data,
options_to_create)
create_raster_tiff_from_reference(extents_dictionary, image,
data.ReadAsArray()) #created in disk
def _get_bundle_from_path(path):
'''
This method tries to instance every bundle using the given directory. When
a bundle is able to identify the files present in the directory the instance
of that bundle is returned to the caller.
'''
return get_bundle_from_path(path, '../../../mapper', PREPROCESSING_PACKAGE)
def test_rasterize_vector_in_memory_and_then_write_to_disk(self):
from madmex.mapper.data.raster import new_options_for_create_raster_from_reference, create_raster_tiff_from_reference
from madmex.core.controller.commands import get_bundle_from_path
from madmex.mapper.data import raster
BUNDLE_PACKAGE = 'madmex.mapper.bundle'
landmask_path = '/Users/erickpalacios/Documents/CONABIO/MADMEXdata/eodata/footprints/country_mexico/'
bundle = get_bundle_from_path(landmask_path, '../../../mapper', BUNDLE_PACKAGE)
extents_dictionary = {u'x_range': 7521.0, u'y_range': 7741.0, u'properties': {u'projection': 'PROJCS["UTM Zone 15, Northern Hemisphere",GEOGCS["Unknown datum based upon the WGS 84 ellipsoid",DATUM["Not specified (based on WGS 84 spheroid)",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",-93],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAMETER["false_northing",0],UNIT["Meter",1]]', u'geotransform': (523185.0, 30.0, 0.0, 2033715.0, 0.0, -30.0)}, u'x_offset': 'array([ 10., 30., 0., 50.])', u'y_offset': 'array([-0., -0., -0., -0.])'}
options_to_create = new_options_for_create_raster_from_reference(extents_dictionary, raster.DATA_SHAPE, (int(extents_dictionary['x_range']), int(extents_dictionary['y_range']), 1), {})
data = create_raster_tiff_from_reference(extents_dictionary, '', None, options_to_create) #created in memory
#new_options_for_create_raster_from_reference(extents_dictionary, raster.DATASET, data, options_to_create)
bundle.rasterize(data, [1], [1])
image = '/Users/erickpalacios/Documents/CONABIO/Tareas/Redisenio_MADMEX/clasificacion_landsat/landsat8/classification/rasterize3.tif'
create_raster_tiff_from_reference(extents_dictionary, image, data.ReadAsArray()) #created in disk
