Пример #1
0
	def run(self):
		for year in self._years:
			for feature in self._feature_collection.get_features():
				path = feature.get('PATH')
				row = feature.get('ROW')
				geometry = ee.Geometry.MultiPolygon(feature.get('GEOMETRY'))

				image_collection = self._collection.filter(ee.Filter.eq('WRS_PATH', int(path))).filter(ee.Filter.eq('WRS_ROW', int(row)))

				for period in self._periods:
					images_by_period = ImageCollection(image_collection).filter_by_period(year, feature.get(period), self._offset)

					if self._apply_brdf:
						images_by_period = images_by_period.apply_brdf()

					if self._clip_geometry:
						images_by_period = images_by_period.clip_geometry()

					if self._apply_mask:
						images_by_period = images_by_period.apply_qamask()

					if self._bands:
						images_by_period = images_by_period.apply_bands(self._bands)

					if self._reducers:
						image_reduced = Image(ImageCollection(images_by_period).apply_reducers(self._reducers))
						image_reduced = image_reduced.rename(image_reduced.bandNames().map(
							lambda band: ee.String(period).cat('_').cat(band))
						)
						final_name = "{0}_{1}_{2}".format(settings.COLLECTION_PREFIX, "{0}{1}".format(path, row), str(year))
						final_image = image_reduced.clip(geometry).select(settings.GENERATION_VARIABLES + settings.GENERATION_EXTRA_VARIABLES).set('year', year).set('system:footprint', geometry)

						self.add_image_in_batch(final_name, {"image": final_image,
															 "year": int(year),
															 "path": int(path),
															 "row": int(row),
															 "geometry": geometry})
					else:
						images_by_period = images_by_period.map(
							lambda image: Image(image).rename(Image(image).bandNames().map(
								lambda band: ee.String(period).cat('_').cat(band))
							)
						)
						images_by_period = images_by_period.toList(settings.MAX_IMAGES)
						for i in xrange(images_by_period.size().getInfo()):
							image = Image(images_by_period.get(i))
							date = image.date().format('yyyyMMdd').getInfo()
							final_name = "{0}_{1}_{2}".format(settings.COLLECTION_PREFIX, "{0}{1}".format(path, row), str(date))
							final_image = Image(image).clip(geometry).select(settings.GENERATION_VARIABLES + settings.GENERATION_EXTRA_VARIABLES).set('year', year).set('system:footprint', geometry)

							self.add_image_in_batch(final_name, {"image": final_image,
																 "year": int(year),
																 "path": int(path),
																 "row": int(row),
																 "geometry": geometry})
Пример #2
0
    def run(self):
        features = self._feature_collection.get_features()
        for year in self._years:
            for feature in features:
                path = feature.get('PATH')
                row = feature.get('ROW')
                geometry = ee.Geometry.MultiPolygon(feature.get('GEOMETRY'))

                neighbors = self.get_neighbors(features, int(path), int(row))
                neighbors_paths = map(lambda x: x[0], neighbors)
                neighbors_rows = map(lambda x: x[1], neighbors)
                neighbors_geometry = ee.Geometry.MultiPolygon(map(lambda x: x[2], neighbors))

                image_collection = self._collection \
                    .filter(ee.Filter.inList('WRS_PATH', ee.List(neighbors_paths))) \
                    .filter(ee.Filter.inList('WRS_ROW', ee.List(neighbors_rows)))

                final_image = Image()
                for period in self._periods:
                    period_dates = feature.get(period)
                    images_by_period = ImageCollection(image_collection).filter_by_period(year, period_dates, self._offset)

                    if self._clip_geometry:
                        images_by_period = images_by_period.clip_geometry()

                    if self._apply_mask:
                        images_by_period = images_by_period.apply_qamask()

                    if self._bands:
                        images_by_period = images_by_period.apply_bands(self._bands)

                    image_reduced = Image(images_by_period.apply_reducers(self._reducers))

                    final_image = final_image.addBands(
                        image_reduced.rename(image_reduced.bandNames().map(lambda band: ee.String(period).cat('_').cat(band)))
                    )

                if settings.GENERATION_EXTRA_BANDS:
                    final_image = final_image.addBands(Image(final_image).get_bands(settings.GENERATION_EXTRA_BANDS))

                final_name = "{0}_{1}_{2}".format(settings.COLLECTION_PREFIX, "{0}{1}".format(path, row), str(year))
                final_image = final_image.select(settings.GENERATION_VARIABLES + settings.GENERATION_EXTRA_VARIABLES).set('year', year)

                self.add_image_in_batch(final_name,
                                        {"image": final_image, "year": int(year), "path": int(path), "row": int(row), "geometry": geometry,
                                         'neighbors': neighbors_geometry})
Пример #3
0
    def run(self):
        for year in self._years:
            for feature in self._feature_collection.get_features():
                path = feature.get('PATH')
                row = feature.get('ROW')
                geometry = ee.Geometry.MultiPolygon(feature.get('geometry'))

                image_collection = self._collection.map(lambda image: image.clip(geometry))

                images = []
                for period in self._periods:
                    period_dates = feature.get(period)

                    images_by_period = ImageCollection(image_collection).filter_by_period(year, period_dates, self._offset)

                    images_by_period = images_by_period.apply_smoothing()

                    if self._clip_geometry:
                        images_by_period = images_by_period.clip_geometry()

                    if self._apply_mask:
                        images_by_period = images_by_period.apply_qamask()

                    if self._bands:
                        images_by_period = images_by_period.apply_bands(self._bands)

                    image_reduced = Image(images_by_period.apply_reducers(self._reducers))
                    image_reduced = image_reduced.rename(image_reduced.bandNames().map(
                        lambda band: ee.String(period).cat('_').cat(band))
                    )
                    images.append(image_reduced)

                final_image = Image.cat(images)

                if settings.GENERATION_EXTRA_BANDS:
                    final_image = final_image.addBands(Image(final_image).get_bands(settings.GENERATION_EXTRA_BANDS))

                final_name = "{0}_{1}_{2}".format(settings.COLLECTION_PREFIX, "{0}{1}".format(path, row), str(year))
                final_image = final_image.clip(geometry).select(settings.GENERATION_VARIABLES + settings.GENERATION_EXTRA_VARIABLES).set('year',
                                                                                                                                         year).set(
                    'system:footprint', geometry)

                self.add_image_in_batch(final_name,
                                        {"image": final_image, "year": int(year), "path": int(path), "row": int(row), "geometry": geometry})