Guest User

Untitled

a guest
Dec 14th, 2017
75
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.44 KB | None | 0 0
  1. import numpy as np
  2. data = np.genfromtxt('sn_bound_10deg.txt', skip_header = 7, skip_footer = 3)
  3. lat,lon = -25.439538,149.053418
  4. in_tile = False
  5. i = 0
  6. while(not in_tile):
  7. #print data[i,0],data[i,1],data[i, 4],data[i, 5],data[i, 2],data[i, 3]
  8. in_tile = lat >= data[i, 4] and lat <= data[i, 5] and lon >= data[i, 2] and lon <= data[i, 3]
  9. i += 1
  10.  
  11. vert = data[i-1, 0]
  12. horiz = data[i-1, 1]
  13. print('Vertical Tile:', vert, 'Horizontal Tile:', horiz)
Add Comment
Please, Sign In to add comment