Guest User

Untitled

a guest
Jun 18th, 2018
80
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 2.56 KB | None | 0 0
  1. class Parcel(models.Model):
  2.  
  3. objectid = models.BigIntegerField()
  4. apn = models.CharField(max_length=80)
  5. lot_size = models.BigIntegerField()
  6. neibrhc = models.CharField(max_length=80)
  7. street_nbr = models.CharField(max_length=80)
  8. street_nam = models.CharField(max_length=80)
  9. city = models.CharField(max_length=80)
  10. zip = models.BigIntegerField()
  11. tra = models.CharField(max_length=80)
  12. subdivisio = models.CharField(max_length=80)
  13. lot = models.CharField(max_length=80)
  14. landuse = models.CharField(max_length=80)
  15. lu_general = models.CharField(max_length=80)
  16. lu_specif = models.CharField(max_length=80)
  17. lu_detail = models.CharField(max_length=80)
  18. lu_use = models.CharField(max_length=80)
  19. lu_sec_use = models.CharField(max_length=80)
  20. area = models.FloatField()
  21. perimeter = models.FloatField()
  22. shape_are = models.FloatField()
  23. shape_len = models.FloatField()
  24. geom = models.MultiPolygonField(srid=3492)
  25.  
  26. def __str__(self):
  27. return 'City: %s' % self.city
  28.  
  29.  
  30. class Primary(models.Model):
  31. id = models.CharField(max_length=12,
  32. primary_key=True)
  33. address = models.CharField(max_length=30)
  34. city = models.CharField(max_length=30)
  35. zip = models.CharField(max_length=10)
  36. state = models.CharField(max_length=2)
  37. unit = models.CharField(max_length=13)
  38. primpgon = models.BigIntegerField()
  39. numpgons = models.BigIntegerField()
  40. x = models.FloatField()
  41. y = models.FloatField()
  42. censusbloc = models.CharField(max_length=15)
  43. objectid = models.BigIntegerField()
  44. geom = models.MultiPolygonField(srid=4326)
  45.  
  46.  
  47.  
  48.  
  49. from django.contrib.gis.geos import Polygon
  50. parcel = Parcel.objects.filter(objectid="22520").values()
  51. geojson = parcel[0]['geom'].geojson
  52. format_l = ast.literal_eval(geojson)
  53. coordinates_data = format_l['coordinates'][0][0]
  54. poly = Polygon(coordinates_data, srid=3492)
  55. poly_parcel = poly.transform(4326, clone=True)
  56. polygon = poly_parcel.coords
  57. print(polygon)
  58.  
  59.  
  60. (((-121.49440799279022, 38.55725848782723),
  61. (-121.49438844924772, 38.557303481514126),
  62. (-121.4943760310021, 38.5573320694682),
  63. (-121.49436263531841, 38.557362909896675),
  64. (-121.49402385986245, 38.557269114460084),
  65. (-121.49406987333441, 38.55716268909225),
  66. (-121.49440799279022, 38.55725848782723)),)
  67.  
  68. Primary.objects.filter(geom__contains=polygon).values()
  69. Primary.objects.filter(geom__contained=polygon).values()
Add Comment
Please, Sign In to add comment