Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1399 1716 (N46W21 L=259)  1695    211        3 (1301/1657-1496/1753) -10.0 W12-W30 N42-N49
      2  293 1519 (N19E63 L=175)  2660    169        3 ( 258/1434- 328/1598) -54.0 E70-E58 N13-N24
      3 1230 1411 (N20W08 L=247)  1278    332        1 (1097/1306-1343/1508) -27.9 E01-W16 N15-N25
      4 1771 1390 (N28W52 L=290)  1014    177        2 (1722/1317-1815/1454) -48.8 W47-W56 N23-N31
      5 1860 1241 (N20W60 L=298)  1263    251        0 (1840/1124-1877/1359)  78.7 W56-W65 N12-N27
      6 1653 1306 (N20W39 L=277)   714    133        1 (1613/1268-1701/1355) -23.2 W35-W43 N18-N22
      7 1767 1304 (N22W50 L=288)   571     77        2 (1720/1281-1803/1325)  14.0 W44-W53 N20-N23
      8  142 1229 (N00E71 L=167)   804     61        2 ( 120/1200- 169/1261)  18.4 E74-E66 S01-N02
      9  145 1029 (S12E67 L=171)   728     80        3 ( 123/ 997- 177/1074)  76.0 E70-E62 S13-S09
     10  336 1014 (S12E46 L=192)   447    101        1 ( 307/ 965- 371/1064)  56.3 E48-E43 S15-S09
     11  775  890 (S16E14 L=224)   195     96        0 ( 736/ 870- 821/ 905)  18.4 E16-E10 S17-S15
     12  223  843 (S23E59 L=180)   459     78        1 ( 186/ 829- 260/ 863)   0.0 E63-E54 S23-S21
     13 1390  592 (S27W31 L=269)  2156    498        3 (1211/ 462-1573/ 748)  52.3 W22-W40 S37-S14
     14  502  701 (S30E33 L=205)   553    120        0 ( 462/ 654- 546/ 745)  36.9 E37-E29 S33-S27
     15 1688  675 (S17W51 L=290)   383     84        1 (1654/ 646-1725/ 710) -14.0 W47-W55 S18-S15
     16 1701  592 (S21W57 L=295)   643    114        1 (1656/ 551-1738/ 640)  66.8 W54-W59 S24-S17
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y