"""
from sqlalchemy import and_
from db.common import Session
from db.n2k import *
if __name__ == '__main__':
tgt = r"D:\work\ms.monina\wp4\florkart_2012_occurrences_eunis_characteristic_species_n2k_areas_germany_all.txt"
tgt = r"z:\bla.txt"
output = list()
session = Session()
# retrieving all habitat types
habitat_types = session.query(HabitatType).all()
habitat_types = session.query(HabitatType).filter(HabitatType.type_id == '2310').all()
for ht in habitat_types[:]:
habitat_type = ht.type_id
# retrieving all sites with the current habitat type
sites_with_habitat_type = session.query(SitesHabitatTypes).filter(
SitesHabitatTypes.type_id == habitat_type).all()
# bailing out if no sites were found
if not len(sites_with_habitat_type):
continue
print "Working on habitat type %s (%s)..." % (habitat_type, ht.shortname)
if len(tokens) > 1:
ht = tokens[0]
ht_map_ids[ht] = list(set(tokens[1].split(", ")))
# caching centroids
m_centroids = cache_centroids(md_src)
q_centroids = cache_centroids(qd_src)
# preparing other stuff
session = Session()
output = list()
wkt_output = list()
for ht_id in sorted(ht_map_ids.keys()):
#for ht_id in ['2310']:
ht_name = session.query(HabitatType).filter(
HabitatType.type_id == ht_id).one().shortname
print "Working on %s (%s)..." % (ht_id, ht_name)
maxent_ht_name = " ".join((ht_id, ht_name)).lower().replace(" ", "_")
all_map_ids = ht_map_ids[ht_id]
# retrieving all mtb tiles where current habitat type is present
m_ids = [map_id for map_id in all_map_ids if map_id.endswith('0')]
# retrieving all quadrant tiles where current habitat type is present
q_ids = [map_id for map_id in all_map_ids if not map_id.endswith('0')]
# generalizing quadrant ids to corresponding mtb ids, i.e. 25101 -> 25100
mq_ids = set([map_id[:-1] + '0' for map_id in q_ids])
# calculating intersection between mtb ids and mtb/quadrant ids
mq_intersect = [map_id for map_id in m_ids if map_id in mq_ids]
print "\t+ Typical species present in %d map tiles overall (%d quadrant tiles, %d mtb tile(s))" % (
if len(tokens) > 1:
ht = tokens[0]
ht_map_ids[ht] = list(set(tokens[1].split(", ")))
# caching centroids
m_centroids = cache_centroids(md_src)
q_centroids = cache_centroids(qd_src)
# preparing other stuff
session = Session()
output = list()
wkt_output = list()
for ht_id in sorted(ht_map_ids.keys()):
#for ht_id in ['2310']:
ht_name = session.query(HabitatType).filter(HabitatType.type_id == ht_id).one().shortname
print "Working on %s (%s)..." % (ht_id, ht_name)
maxent_ht_name = " ".join((ht_id, ht_name)).lower().replace(" ", "_")
all_map_ids = ht_map_ids[ht_id]
# retrieving all mtb tiles where current habitat type is present
m_ids = [map_id for map_id in all_map_ids if map_id.endswith('0')]
# retrieving all quadrant tiles where current habitat type is present
q_ids = [map_id for map_id in all_map_ids if not map_id.endswith('0')]
# generalizing quadrant ids to corresponding mtb ids, i.e. 25101 -> 25100
mq_ids = set([map_id[:-1] + '0' for map_id in q_ids])
# calculating intersection between mtb ids and mtb/quadrant ids
mq_intersect = [map_id for map_id in m_ids if map_id in mq_ids]
print "\t+ Typical species present in %d map tiles overall (%d quadrant tiles, %d mtb tile(s))" % (len(all_map_ids), len(q_ids), len(m_ids))
"""
from sqlalchemy import and_
from db.common import Session
from db.n2k import *
if __name__ == '__main__':
tgt = r"D:\work\ms.monina\wp4\florkart_2012_occurrences_eunis_characteristic_species_n2k_areas_germany_all.txt"
tgt = r"z:\bla.txt"
output = list()
session = Session()
# retrieving all habitat types
habitat_types = session.query(HabitatType).all()
habitat_types = session.query(HabitatType).filter(
HabitatType.type_id == '2310').all()
for ht in habitat_types[:]:
habitat_type = ht.type_id
# retrieving all sites with the current habitat type
sites_with_habitat_type = session.query(SitesHabitatTypes).filter(
SitesHabitatTypes.type_id == habitat_type).all()
# bailing out if no sites were found
if not len(sites_with_habitat_type):
continue
print "Working on habitat type %s (%s)..." % (habitat_type,
