Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  631 1804 (N31E58 L=230)  1257     99        1 ( 609/1756- 656/1852) -29.7 E68-E54 N26-N34
      2  701 1774 (N30E49 L=239)   426     46        1 ( 675/1755- 731/1798)  11.3 E51-E45 N28-N31
      3  453 1697 (N19E64 L=224)  2601    240        3 ( 375/1628- 510/1787) -74.1 E70-E60 N13-N27
      4  995 1714 (N34E24 L=264)   670    100        1 ( 968/1666-1026/1765) -71.6 E26-E22 N29-N38
      5  427 1591 (N12E58 L=230)   400     59        1 ( 395/1573- 455/1609)  45.0 E60-E55 N10-N13
      6 1419 1541 (N35W10 L=298)   685    114        1 (1380/1497-1448/1593)  77.5 W07-W11 N31-N40
      7  277 1413 (S02E63 L=225)   732     66        1 ( 248/1383- 311/1438)   7.1 E66-E58 S03-N00
      8 1381 1327 (N20W12 L=300)   362    124        0 (1336/1269-1416/1383)  77.5 W10-W13 N15-N24
      9  166 1230 (S14E70 L=218)  3864    230        5 ( 112/1143- 238/1303)  34.7 E81-E62 S19-S09
     10 1848  844 (N08W61 L=349)   578     69        1 (1814/ 820-1876/ 875)  33.7 W57-W64 N06-N10
     11 1067  613 (S29W15 L=303)   920    276        1 ( 941/ 592-1182/ 660)  37.6 W07-W20 S32-S22
     12 1623  527 (S16W56 L=344)   334     75        0 (1588/ 507-1657/ 558)  68.2 W53-W57 S17-S12
  

    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