Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  522 1702 (N45E51 L=142)  1855    360        4 ( 479/1638- 582/1766) -31.0 E60-E43 N40-N50
      2  973 1658 (N45E08 L=185)   296     37        1 ( 948/1637- 999/1674)   0.0 E10-E05 N43-N46
      3  479 1577 (N36E47 L=146)   348     31        0 ( 456/1555- 504/1599)  33.7 E48-E45 N33-N37
      4  239 1292 (N15E61 L=132)  5444    733       12 ( 163/1162- 320/1482)  79.2 E71-E54 N06-N27
      5 1837 1348 (N25W64 L=257)   605    132        1 (1822/1318-1854/1375)  45.0 W61-W66 N22-N26
      6  335 1239 (N12E50 L=143)   419     59        0 ( 317/1206- 354/1268)  90.0 E51-E48 N10-N14
      7  728 1051 (N03E19 L=174)   229     59        0 ( 699/1017- 756/1079) -36.9 E20-E16 N01-N04
      8  228  846 (S13E59 L=134)  2552    318        7 ( 179/ 735- 277/ 931) -76.0 E65-E53 S19-S07
      9 1731  827 (S08W51 L=243)   482     87        1 (1705/ 786-1748/ 865) -59.0 W47-W52 S10-S05
     10 1051  455 (S35W05 L=197)   758    128        1 ( 982/ 434-1117/ 475)  -6.3 W00-W09 S36-S33
     11  689  423 (S39E24 L=169)  1011    223        3 ( 598/ 390- 792/ 455)  13.2 E32-E15 S42-S36
  

    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