def idl(**kwargs): for key, value in kwargs.items(): if key == "infile" and value is not None: infile = value try: if infile.endswith(".geojson"): aoi_resp = multipoly(infile) try: for things in aoi_resp["features"]: ovall.append(things["geometry"]["coordinates"]) except Exception: for things in json.loads(aoi_resp)["features"]: ovall.append(things["geometry"]["coordinates"]) # print(list_depth(ovall)) aoi_geom = ovall if (list_depth(ovall) == 2 and json.loads(aoi_resp) ["features"][0]["geometry"]["type"] == "MultiPolygon"): temp["coordinates"] = aoi_geom[0] elif (list_depth(ovall) == 2 and json.loads(aoi_resp)["features"][0]["geometry"] ["type"] != "MultiPolygon"): aoi_geom = ovall[0][0] temp["type"] = "Polygon" temp["coordinates"] = aoi_geom elif list_depth(ovall) == 1: aoi_geom = ovall[0] temp["type"] = "Polygon" temp["coordinates"] = aoi_geom elif list_depth(ovall) == 0: aoi_geom = ovall temp["type"] = "Polygon" temp["coordinates"] = aoi_geom else: print( "Please check GeoJSON: Could not parse coordinates" ) elif infile.endswith(".json"): with open(infile) as aoi: aoi_resp = json.load(aoi) aoi_geom = aoi_resp["config"][0]["config"][ "coordinates"] elif infile.endswith(".kml"): getcoord = kml2coord(infile) aoi_geom = getcoord["features"][0]["geometry"][ "coordinates"] except Exception as e: print("Could not parse geometry") print(e) if key == "item" and value is not None: try: item = value except Exception as e: sys.exit(e) if key == "start" and value is not None: try: start = time2utc(value) st = filters.date_range("acquired", gte=start) except Exception as e: sys.exit(e) if key == "end" and value is not None: try: end = time2utc(value) ed = filters.date_range("acquired", lte=end) except Exception as e: sys.exit(e) if key == "asset" and value is not None: try: asset = value except Exception as e: sys.exit(e) if key == "cmin": if value == None: try: cmin = 0 except Exception as e: print(e) if value is not None: try: cmin = float(value) except Exception as e: print(e) if key == "cmax": if value == None: try: cmax = 1 except Exception as e: print(e) elif value is not None: try: cmax = float(value) except Exception as e: print(e) if key == "num": if value is not None: num = value elif value == None: num = 1000000 if key == "outfile" and value is not None: outfile = value try: open(outfile, "w") except Exception as e: sys.exit(e) if key == "ovp": if value is not None: ovp = int(value) elif value == None: ovp = 0.01 if key == "filters" and value is not None: for items in value: ftype = items.split(":")[0] if ftype == "string": try: fname = items.split(":")[1] fval = items.split(":")[2] # stbase={'config': [], 'field_name': [], 'type': 'StringInFilter'} stbase["config"] = fval.split(",") # fval stbase["field_name"] = fname except Exception as e: print(e) elif ftype == "range": fname = items.split(":")[1] fgt = items.split(":")[2] flt = items.split(":")[3] # rbase={'config': {'gte': [], 'lte': []},'field_name': [], 'type': 'RangeFilter'} rbase["config"]["gte"] = int(fgt) rbase["config"]["lte"] = int(flt) rbase["field_name"] = fname print("Running search for a maximum of: " + str(num) + " assets") l = 0 [head, tail] = os.path.split(outfile) sgeom = filters.geom_filter(temp) aoi_shape = shape(temp) if not aoi_shape.is_valid: aoi_shape = aoi_shape.buffer(0) # print('Your Input Geometry is invalid & may have issues:A valid Polygon may not possess anyoverlapping exterior or interior rings.'+'\n') date_filter = filters.date_range("acquired", gte=start, lte=end) cloud_filter = filters.range_filter("cloud_cover", gte=cmin, lte=cmax) asset_filter = filters.permission_filter("assets." + str(asset.split(",")[0]) + ":download") # print(rbase) # print(stbase) if len(rbase["field_name"]) != 0 and len(stbase["field_name"]) != 0: and_filter = filters.and_filter(date_filter, cloud_filter, asset_filter, sgeom, stbase, rbase) elif len(rbase["field_name"]) == 0 and len(stbase["field_name"]) != 0: and_filter = filters.and_filter(date_filter, cloud_filter, asset_filter, sgeom, stbase) elif len(rbase["field_name"]) != 0 and len(stbase["field_name"]) == 0: and_filter = filters.and_filter(date_filter, cloud_filter, asset_filter, sgeom, rbase) elif len(rbase["field_name"]) == 0 and len(stbase["field_name"]) == 0: and_filter = filters.and_filter(date_filter, cloud_filter, asset_filter, sgeom) item_types = [item] req = filters.build_search_request(and_filter, item_types) res = client.quick_search(req) for things in res.items_iter( 1000000): # A large number as max number to check against try: all_assets = [ assets.split(":")[0].replace("assets.", "") for assets in things["_permissions"] ] if things["properties"]["quality_category"] == "standard" and all( elem in all_assets for elem in asset.split(",")): itemid = things["id"] footprint = things["geometry"] s = shape(footprint) if item.startswith("SkySat"): epsgcode = "3857" else: epsgcode = things["properties"]["epsg_code"] if aoi_shape.area > s.area: intersect = (s).intersection(aoi_shape) elif s.area >= aoi_shape.area: intersect = (aoi_shape).intersection(s) proj_transform = pyproj.Transformer.from_proj( pyproj.Proj(4326), pyproj.Proj(epsgcode), always_xy=True ).transform # always_xy determines correct coord order print( "Processing " + str(len(ar) + 1) + " items with total area " + str("{:,}".format(round(sum(far)))) + " sqkm", end="\r", ) if (transform(proj_transform, (aoi_shape)).area > transform( proj_transform, s).area): if (transform(proj_transform, intersect).area / transform(proj_transform, s).area * 100) >= ovp: ar.append( transform(proj_transform, intersect).area / 1000000) far.append(transform(proj_transform, s).area / 1000000) with open(outfile, "a") as csvfile: writer = csv.writer(csvfile, delimiter=",", lineterminator="\n") writer.writerow([itemid]) elif (transform(proj_transform, s).area >= transform( proj_transform, aoi_shape).area): if (transform(proj_transform, intersect).area / transform( proj_transform, aoi_shape).area * 100) >= ovp: ar.append( transform(proj_transform, intersect).area / 1000000) far.append(transform(proj_transform, s).area / 1000000) with open(outfile, "a") as csvfile: writer = csv.writer(csvfile, delimiter=",", lineterminator="\n") writer.writerow([itemid]) if int(len(ar)) == int(num): break except Exception as e: pass except (KeyboardInterrupt, SystemExit) as e: print("\n" + "Program escaped by User") sys.exit() num_lines = sum( 1 for line in open(os.path.join(head, tail.split(".")[0] + ".csv"))) print("Total number of item ids written to " + str( os.path.join(head, tail.split(".")[0] + ".csv") + " ===> " + str(num_lines))) print("Total estimated cost to quota: " + str("{:,}".format(round(sum(far)))) + " sqkm") print("Total estimated cost to quota if clipped: " + str("{:,}".format(round(sum(ar)))) + " sqkm")
def to_permissions(asset_types): return filters.permission_filter( *['assets.%s:download' % a for a in asset_types])
def filters(self): populated_filters = [] start_date = None end_date = None start_datetime = None end_datetime = None if not self.startDateEdit.dateTime().isNull(): start_datetime = self.startDateEdit.dateTime() start_date = start_datetime.toString(Qt.ISODate) if not self.endDateEdit.dateTime().isNull(): end_datetime = self.endDateEdit.dateTime() end_date = end_datetime.toString(Qt.ISODate) if start_datetime and end_datetime: if start_datetime < end_datetime: date_filter = date_range('acquired', gte=start_date, lte=end_date) populated_filters.append(date_filter) else: self._show_message('Start date later than end date.', level=Qgis.Warning, duration=10) elif start_date: start_date_filter = date_range('acquired', gte=start_date) populated_filters.append(start_date_filter) elif end_date: end_date_filter = date_range('acquired', lte=end_date) populated_filters.append(end_date_filter) # TODO: double check actual domain/range of sliders sliders = self.frameRangeSliders.findChildren( PlanetExplorerRangeSlider) for slider in sliders: slide_filter = None range_low, range_high = slider.range() if slider.filter_key == 'cloud_cover': range_low /= 100.0 range_high /= 100.0 if range_low != slider.min and range_high != slider.max: slide_filter = range_filter(slider.filter_key, gte=range_low, lte=range_high) elif range_low != slider.min: slide_filter = range_filter(slider.filter_key, gte=range_low) elif range_high != slider.max: slide_filter = range_filter(slider.filter_key, lte=range_high) if slide_filter: populated_filters.append(slide_filter) s_ids = self.leStringIDs.text() if s_ids: ids_actual = [] s_ids.replace(" ", "") for s_id in s_ids.split(','): for text_chunk in s_id.split(':'): for pattern in self.id_regex: if pattern.match(text_chunk): ids_actual.append(text_chunk) if ids_actual: s_ids_list = ['id'] s_ids_list.extend(ids_actual) string_ids_filter = string_filter(*s_ids_list) populated_filters.append(string_ids_filter) else: self._show_message('No valid ID present', level=Qgis.Warning, duration=10) if self.chkBxCanDownload.isChecked(): dl_permission_filter = permission_filter('assets:download') populated_filters.append(dl_permission_filter) # Ground_control can be 'true', 'false, or a numeric value # Safest to check for not 'false' if self.chkBxGroundControl.isChecked(): gc_filter = not_filter(string_filter('ground_control', 'false')) populated_filters.append(gc_filter) return populated_filters
def idc(idlist, item, asset, geometry): try: if geometry.endswith('.geojson'): with open(geometry) as aoi: aoi_resp = json.load(aoi) if list_depth(aoi_resp['features'][0]['geometry'] ['coordinates']) == 0: aoi_geom = aoi_resp['features'][0]['geometry'][ 'coordinates'] elif list_depth(aoi_resp['features'][0]['geometry'] ['coordinates']) == 1: aoi_geom = aoi_resp['features'][0]['geometry'][ 'coordinates'][0] else: print('Please check GeoJSON: Could not parse coordinates') elif geometry.endswith('.json'): with open(geometry) as aoi: aoi_resp = json.load(aoi) aoi_geom = aoi_resp['config'][0]['config']['coordinates'] elif geometry.endswith('.kml'): getcoord = kml2coord(geometry) aoi_geom = getcoord except Exception as e: print('Could not parse geometry') print(e) l = [] stbase = {'config': [], 'field_name': 'id', 'type': 'StringInFilter'} with open(idlist) as f: try: reader = csv.reader(f) for row in reader: item_id = row[0] if str(item_id.isalpha()) == 'False': l.append(item_id) except Exception as e: print(e) stbase['config'] = l temp['coordinates'] = aoi_geom sgeom = filters.geom_filter(temp) aoi_shape = shape(temp) asset_filter = filters.permission_filter('assets.' + str(asset) + ':download') and_filter = filters.and_filter(asset_filter, stbase) item_types = [item] req = filters.build_search_request(and_filter, item_types) res = client.quick_search(req) num_lines = sum( 1 for line in open(os.path.join(idlist.split('.')[0] + '.csv'))) n = 1 print('Checking permission for assets in idlist...') for things in res.items_iter( num_lines): # A large number as max number to check against itemid = things['id'] #print('Processing '+str(n)+' of '+str(num_lines)) n = n + 1 footprint = things["geometry"] s = shape(footprint) if item.startswith('SkySat'): epsgcode = '3857' else: epsgcode = things['properties']['epsg_code'] if aoi_shape.area > s.area: intersect = (s).intersection(aoi_shape) elif s.area >= aoi_shape.area: intersect = (aoi_shape).intersection(s) proj = partial(pyproj.transform, pyproj.Proj(init='epsg:4326'), pyproj.Proj(init='epsg:' + str(epsgcode))) if transform(proj, aoi_shape).area > transform(proj, s).area: ar.append(transform(proj, intersect).area / 1000000) far.append(transform(proj, s).area / 1000000) elif transform(proj, s).area > transform(proj, aoi_shape).area: ar.append(transform(proj, intersect).area / 1000000) far.append(transform(proj, s).area / 1000000) print('\n' + 'Total estimated cost to quota: ' + str("{:,}".format(round(sum(far)))) + ' sqkm') print('Total estimated cost to quota if clipped: ' + str("{:,}".format(round(sum(ar)))) + ' sqkm')
def idl(**kwargs): for key, value in kwargs.items(): if key == 'infile' and value is not None: infile = value try: if infile.endswith('.geojson'): with open(infile) as aoi: aoi_resp = json.load(aoi) for things in aoi_resp['features']: ovall.append(things['geometry']['coordinates']) #print(list_depth(ovall)) if len(ovall) > 1: aoi_geom = ovall else: if list_depth(ovall) == 0: aoi_geom = ovall elif list_depth(ovall) == 1: aoi_geom = ovall[0] elif list_depth(ovall) == 2: aoi_geom = ovall[0][0] else: print( 'Please check GeoJSON: Could not parse coordinates' ) aoi_geom == aoi_geom #print(aoi_geom) elif infile.endswith('.json'): with open(infile) as aoi: aoi_resp = json.load(aoi) aoi_geom = aoi_resp['config'][0]['config'][ 'coordinates'] elif infile.endswith('.kml'): getcoord = kml2coord(infile) aoi_geom = getcoord except Exception as e: print('Could not parse geometry') print(e) if key == 'item' and value is not None: try: item = value except Exception as e: sys.exit(e) if key == 'start' and value is not None: try: start = value st = filters.date_range('acquired', gte=start) except Excpetion as e: sys.exit(e) if key == 'end' and value is not None: end = value ed = filters.date_range('acquired', lte=end) if key == 'asset' and value is not None: try: asset = value except Exception as e: sys.exit(e) if key == 'cmin': if value == None: try: cmin = 0 except Exception as e: print(e) if value is not None: try: cmin = float(value) except Exception as e: print(e) if key == 'cmax': if value == None: try: cmax = 1 except Exception as e: print(e) elif value is not None: try: cmax = float(value) except Exception as e: print(e) if key == 'num': if value is not None: num = value elif value == None: num = 1000000 if key == 'outfile' and value is not None: outfile = value try: open(outfile, 'w') except Exception as e: sys.exit(e) if key == 'ovp': if value is not None: ovp = int(value) elif value == None: ovp = 1 if key == 'filters' and value is not None: for items in value: ftype = items.split(':')[0] if ftype == 'string': try: fname = items.split(':')[1] fval = items.split(':')[2] #stbase={'config': [], 'field_name': [], 'type': 'StringInFilter'} stbase['config'] = fval.split(',') #fval stbase['field_name'] = fname except Exception as e: print(e) elif ftype == 'range': fname = items.split(':')[1] fgt = items.split(':')[2] flt = items.split(':')[3] #rbase={'config': {'gte': [], 'lte': []},'field_name': [], 'type': 'RangeFilter'} rbase['config']['gte'] = int(fgt) rbase['config']['lte'] = int(flt) rbase['field_name'] = fname print('Running search for a maximum of: ' + str(num) + ' assets') l = 0 [head, tail] = os.path.split(outfile) if len(ovall) > 1: temp = {"coordinates": [], "type": "MultiPolygon"} temp['coordinates'] = aoi_geom else: temp = tempsingle temp['coordinates'] = aoi_geom sgeom = filters.geom_filter(temp) aoi_shape = shape(temp) date_filter = filters.date_range('acquired', gte=start, lte=end) cloud_filter = filters.range_filter('cloud_cover', gte=cmin, lte=cmax) asset_filter = filters.permission_filter('assets.' + str(asset) + ':download') # print(rbase) # print(stbase) if len(rbase['field_name']) != 0 and len(stbase['field_name']) != 0: and_filter = filters.and_filter(date_filter, cloud_filter, asset_filter, sgeom, stbase, rbase) elif len(rbase['field_name']) == 0 and len(stbase['field_name']) != 0: and_filter = filters.and_filter(date_filter, cloud_filter, asset_filter, sgeom, stbase) elif len(rbase['field_name']) != 0 and len(stbase['field_name']) == 0: and_filter = filters.and_filter(date_filter, cloud_filter, asset_filter, sgeom, rbase) elif len(rbase['field_name']) == 0 and len(stbase['field_name']) == 0: and_filter = filters.and_filter(date_filter, cloud_filter, asset_filter, sgeom) item_types = [item] req = filters.build_search_request(and_filter, item_types) res = client.quick_search(req) for things in res.items_iter( 1000000): # A large number as max number to check against itemid = things['id'] footprint = things["geometry"] s = shape(footprint) if item.startswith('SkySat'): epsgcode = '3857' else: epsgcode = things['properties']['epsg_code'] if aoi_shape.area > s.area: intersect = (s).intersection(aoi_shape) elif s.area >= aoi_shape.area: intersect = (aoi_shape).intersection(s) proj = partial(pyproj.transform, pyproj.Proj(init='epsg:4326'), pyproj.Proj(init='epsg:' + str(epsgcode))) if transform(proj, aoi_shape).area > transform(proj, s).area: if (transform(proj, intersect).area / transform(proj, s).area * 100) >= ovp: ar.append(transform(proj, intersect).area / 1000000) far.append(transform(proj, s).area / 1000000) with open(outfile, 'a') as csvfile: writer = csv.writer(csvfile, delimiter=',', lineterminator='\n') writer.writerow([itemid]) elif transform(proj, s).area > transform(proj, aoi_shape).area: if (transform(proj, intersect).area / transform(proj, aoi_shape).area * 100) >= ovp: ar.append(transform(proj, intersect).area / 1000000) far.append(transform(proj, s).area / 1000000) with open(outfile, 'a') as csvfile: writer = csv.writer(csvfile, delimiter=',', lineterminator='\n') writer.writerow([itemid]) if int(len(ar)) == int(num): break num_lines = sum( 1 for line in open(os.path.join(head, tail.split('.')[0] + '.csv'))) print('Total number of assets written to ' + str( os.path.join(head, tail.split('.')[0] + '.csv') + ' ===> ' + str(num_lines))) print('Total estimated cost to quota: ' + str("{:,}".format(round(sum(far)))) + ' sqkm') print('Total estimated cost to quota if clipped: ' + str("{:,}".format(round(sum(ar)))) + ' sqkm')
def filters(self): populated_filters = [] start_qdate = None end_qdate = None start_date = None end_date = None if not self.startDateEdit.dateTime().isNull(): start_qdate = self.startDateEdit.date() start_date = start_qdate.toString(Qt.ISODate) if not self.endDateEdit.dateTime().isNull(): end_qdate = self.endDateEdit.date().addDays(1) end_date = end_qdate.toString(Qt.ISODate) if start_qdate and end_qdate: if start_qdate < end_qdate: date_filter = date_range("acquired", gte=start_date, lte=end_date) populated_filters.append(date_filter) else: self._show_message("Start date later than end date.", level=Qgis.Warning, duration=10) elif start_date: start_date_filter = date_range("acquired", gte=start_date) populated_filters.append(start_date_filter) elif end_date: end_date_filter = date_range("acquired", lte=end_date) populated_filters.append(end_date_filter) # TODO: double check actual domain/range of sliders sliders = self.frameRangeSliders.findChildren( PlanetExplorerRangeSlider) for slider in sliders: slide_filter = None range_low, range_high = slider.range() if slider.filter_key == "cloud_cover": range_low /= 100.0 range_high /= 100.0 slider_max = 1.0 else: slider_max = slider.max if range_low != slider.min and range_high != slider_max: slide_filter = range_filter(slider.filter_key, gte=range_low, lte=range_high) elif range_low != slider.min: slide_filter = range_filter(slider.filter_key, gte=range_low) elif range_high != slider_max: slide_filter = range_filter(slider.filter_key, lte=range_high) if slide_filter: populated_filters.append(slide_filter) s_ids = self.leStringIDs.text() if s_ids: ids_actual = [] s_ids.replace(" ", "") for s_id in s_ids.split(","): for text_chunk in s_id.split(":"): for pattern in self.id_regex: if pattern.match(text_chunk): ids_actual.append(text_chunk) if ids_actual: s_ids_list = ["id"] s_ids_list.extend(ids_actual) string_ids_filter = string_filter(*s_ids_list) populated_filters.append(string_ids_filter) else: self._show_message("No valid ID present", level=Qgis.Warning, duration=10) instruments = [] for chk in [self.chkPs2, self.chkPs2Sd, self.chkPsbSd]: if chk.isChecked(): instruments.append(chk.text()) if instruments: instrument_filter = string_filter("instrument", *instruments) populated_filters.append(instrument_filter) server_filters = [] if self.chkBxCanDownload.isChecked(): dl_permission_filter = permission_filter("assets:download") server_filters.append(dl_permission_filter) # Ground_control can be 'true', 'false, or a numeric value # Safest to check for not 'false' if self.chkBxGroundControl.isChecked(): gc_filter = not_filter(string_filter("ground_control", "false")) server_filters.append(gc_filter) server_filters.extend([ f for f in populated_filters if f["field_name"] not in LOCAL_FILTERS ]) local_filters = [ f for f in populated_filters if f["field_name"] in LOCAL_FILTERS ] return server_filters, local_filters
def to_permissions(asset_types): return filters.permission_filter(*[ 'assets.%s:download' % a for a in asset_types ])
def idc(idlist, item, asset, geometry): if geometry is not None: try: if geometry.endswith(".geojson"): aoi_resp = multipoly(geometry) try: for things in aoi_resp["features"]: ovall.append(things["geometry"]["coordinates"]) except Exception: for things in json.loads(aoi_resp)["features"]: ovall.append(things["geometry"]["coordinates"]) # print(list_depth(ovall)) aoi_geom = ovall if (list_depth(ovall) == 2 and json.loads(aoi_resp)["features"][0]["geometry"]["type"] == "MultiPolygon"): temp["coordinates"] = aoi_geom[0] elif (list_depth(ovall) == 2 and json.loads(aoi_resp)["features"][0]["geometry"]["type"] != "MultiPolygon"): aoi_geom = ovall[0][0] temp["type"] = "Polygon" temp["coordinates"] = aoi_geom elif list_depth(ovall) == 1: aoi_geom = ovall[0] temp["type"] = "Polygon" temp["coordinates"] = aoi_geom elif list_depth(ovall) == 0: aoi_geom = ovall temp["type"] = "Polygon" temp["coordinates"] = aoi_geom else: print("Please check GeoJSON: Could not parse coordinates") except Exception as e: print("Could not parse geometry") print(e) l = [] stbase = {"config": [], "field_name": "id", "type": "StringInFilter"} with open(idlist) as f: try: reader = csv.reader(f) for row in reader: item_id = row[0] if str(item_id.isalpha()) == "False": l.append(item_id) except Exception as e: print(e) stbase["config"] = l # temp["coordinates"] = aoi_geom # sgeom = filters.geom_filter(temp) aoi_shape = shape(temp) asset_filter = filters.permission_filter("assets." + str(asset) + ":download") and_filter = filters.and_filter(asset_filter, stbase) item_types = [item] req = filters.build_search_request(and_filter, item_types) res = client.quick_search(req) num_lines = sum( 1 for line in open(os.path.join(idlist.split(".")[0] + ".csv"))) n = 1 print("Checking assets in idlist...") for things in res.items_iter( num_lines): # A large number as max number to check against itemid = things["id"] # print('Processing '+str(n)+' of '+str(num_lines)) n = n + 1 footprint = things["geometry"] s = shape(footprint) if item.startswith("SkySat"): epsgcode = "3857" else: epsgcode = things["properties"]["epsg_code"] if aoi_shape.area > s.area: intersect = (s).intersection(aoi_shape) elif s.area >= aoi_shape.area: intersect = (aoi_shape).intersection(s) proj = pyproj.Transformer.from_proj(pyproj.Proj(4326), pyproj.Proj(epsgcode), always_xy=True).transform if transform(proj, aoi_shape).area > transform(proj, s).area: ar.append(transform(proj, intersect).area / 1000000) far.append(transform(proj, s).area / 1000000) elif transform(proj, s).area > transform(proj, aoi_shape).area: ar.append(transform(proj, intersect).area / 1000000) far.append(transform(proj, s).area / 1000000) elif geometry is None: l = [] stbase = {"config": [], "field_name": "id", "type": "StringInFilter"} with open(idlist) as f: try: reader = csv.reader(f) for row in reader: item_id = row[0] if str(item_id.isalpha()) == "False": l.append(item_id) except Exception as e: print(e) stbase["config"] = l asset_filter = filters.permission_filter("assets." + str(asset) + ":download") and_filter = filters.and_filter(asset_filter, stbase) item_types = [item] req = filters.build_search_request(and_filter, item_types) res = client.quick_search(req) num_lines = sum( 1 for line in open(os.path.join(idlist.split(".")[0] + ".csv"))) n = 1 print("Checking assets in idlist...") for things in res.items_iter( num_lines): # A large number as max number to check against itemid = things["id"] # print('Processing '+str(n)+' of '+str(num_lines)) n = n + 1 footprint = things["geometry"] s = shape(footprint) if item.startswith("SkySat"): epsgcode = "3857" else: epsgcode = things["properties"]["epsg_code"] proj = pyproj.Transformer.from_proj(pyproj.Proj(4326), pyproj.Proj(epsgcode), always_xy=True).transform far.append(transform(proj, s).area / 1000000) print("\n" + "Total estimated cost to quota: " + str("{:,}".format(round(sum(far)))) + " sqkm") if not len(ar) == 0: print("Total estimated cost to quota if clipped: " + str("{:,}".format(round(sum(ar)))) + " sqkm")