Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1049 1808 (N44E29 L=219)  1278    183        2 ( 988/1733-1118/1858) -86.2 E32-E25 N36-N51
      2 1382 1497 (N32W08 L=256)   896    158        1 (1355/1412-1403/1559) -56.3 W04-W10 N27-N37
      3 1615 1392 (N34W27 L=276)  1482    371        0 (1521/1243-1698/1509) -30.1 W17-W36 N27-N40
      4 1173 1401 (N19E02 L=246)   313     92        1 (1149/1352-1204/1443)  90.0 E03-E01 N15-N22
      5 1745 1039 (N17W45 L=293)   893    200        2 (1715/ 955-1783/1127) -41.6 W40-W50 N12-N21
      6 1011  796 (S19W06 L=254)   409    186        1 ( 941/ 768-1081/ 837)  12.5 W00-W09 S21-S15
      7  984  437 (S41W19 L=267)   610    154        2 ( 920/ 397-1017/ 488)  83.7 W15-W20 S45-S36
      8  842  403 (S48W11 L=259)   288     62        0 ( 817/ 384- 877/ 430)   0.0 W07-W14 S49-S46
      9 1220  296 (S41W45 L=293)   484     88        0 (1179/ 270-1261/ 321)   6.3 W40-W49 S41-S39
     10 1465  287 (S32W65 L=313)  1027    106        2 (1403/ 260-1514/ 307)  10.6 W58-W74 S33-S29
  

    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