def get_ndvi(latlng_pair): latlng = [] latlng.append(list(latlng_pair)) l = GeoCoordinatesToWorldwideReferenceSystem() image_map = l.getPathRowToLatLngListMapFromLatLngList(latlng) for path in image_map.keys(): value = image_map[path] for row in value.keys(): directory_path = construct_directory_path(path, row) latlng_pairs = value[row] band_file = get_band_filename(directory_path) xy_pairs = get_xy_from_latlng(latlng_pairs, directory_path+band_file) xy_pairs_bands = get_bands(xy_pairs, directory_path) print band_file print xy_pairs_bands metadata = ParseMetaData(band_file.split('_')[0], directory_path) metadata.get_meta_filename() metadata.process_metadata() m_b4 = metadata.get_radiance_mult_band_4() m_b5 = metadata.get_radiance_mult_band_5() a_b4 = metadata.get_radiance_add_band_4() a_b5 = metadata.get_radiance_add_band_5() ndvi = CalcNDVI(m_b4, m_b5, a_b4, a_b5) band4And5Pairs = map(lambda allBands : [allBands[3], allBands[4]], xy_pairs_bands) return ndvi.computeNDVIList(band4And5Pairs)
def get_ndvi(latlng_pair): latlng = [] latlng.append(list(latlng_pair)) l = GeoCoordinatesToWorldwideReferenceSystem() image_map = l.getPathRowToLatLngListMapFromLatLngList(latlng) for path in image_map.keys(): value = image_map[path] for row in value.keys(): directory_path = construct_directory_path(path, row) latlng_pairs = value[row] band_file = get_band_filename(directory_path) xy_pairs = get_xy_from_latlng(latlng_pairs, directory_path + band_file) xy_pairs_bands = get_bands(xy_pairs, directory_path) print band_file print xy_pairs_bands metadata = ParseMetaData(band_file.split('_')[0], directory_path) metadata.get_meta_filename() metadata.process_metadata() m_b4 = metadata.get_radiance_mult_band_4() m_b5 = metadata.get_radiance_mult_band_5() a_b4 = metadata.get_radiance_add_band_4() a_b5 = metadata.get_radiance_add_band_5() ndvi = CalcNDVI(m_b4, m_b5, a_b4, a_b5) band4And5Pairs = map(lambda allBands: [allBands[3], allBands[4]], xy_pairs_bands) return ndvi.computeNDVIList(band4And5Pairs)
gt, ct = gh.open_and_transform(band_file) pixel_pairs = gh.coordarray_to_pixelarray(latlng_pairs, ct, gt) return pixel_pairs def get_band_filename(directory_path): files = os.listdir(directory_path) return files[0] def get_bands(xy_pairs, directory_path): k = stackGeoTIFF(directory_path) k.load_bands() return k.get_pixel_stacks(xy_pairs) latLngList = [[27.486772, -81.459955],[27.397890, -81.296839],[27.425886, -81.173685],[27.364792, -81.398872],[27.361307, -81.413538]] l = GeoCoordinatesToWorldwideReferenceSystem() image_map = l.getPathRowToLatLngListMapFromLatLngList(latLngList) print image_map xy_pairs = [] for path in image_map.keys(): value = image_map[path] for row in value.keys(): directory_path = construct_directory_path(path, row) latlng_pairs = value[row] band_file = get_band_filename(directory_path) xy_pairs = get_xy_from_latlng(latlng_pairs, directory_path+band_file) xy_pairs_bands = get_bands(xy_pairs, directory_path) print band_file print xy_pairs_bands metadata = ParseMetaData(band_file.split('_')[0], directory_path)
def get_bands(xy_pairs, directory_path): k = stackGeoTIFF(directory_path) k.load_bands() return k.get_pixel_stacks(xy_pairs) latLngList = [ [27.486772, -81.459955], [27.397890, -81.296839], [27.425886, -81.173685], [27.364792, -81.398872], [27.361307, -81.413538], ] l = GeoCoordinatesToWorldwideReferenceSystem() image_map = l.getPathRowToLatLngListMapFromLatLngList(latLngList) print image_map xy_pairs = [] for path in image_map.keys(): value = image_map[path] for row in value.keys(): directory_path = construct_directory_path(path, row) latlng_pairs = value[row] band_file = get_band_filename(directory_path) xy_pairs = get_xy_from_latlng(latlng_pairs, directory_path + band_file) xy_pairs_bands = get_bands(xy_pairs, directory_path) print band_file print xy_pairs_bands metadata = ParseMetaData(band_file.split("_")[0], directory_path)