Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1407 1359 (N19W21 L= 13)   392    113        2 (1351/1336-1468/1384)  -7.1 W17-W25 N17-N19
      2 1652 1288 (N17W39 L= 30)   844    187        1 (1619/1178-1678/1359) -79.7 W36-W41 N10-N21
      3 1106 1284 (N11W03 L=354)  1049    280        0 ( 985/1188-1242/1357) -23.6 E05-W11 N06-N15
      4 1809 1137 (N10W53 L= 44)   499    101        2 (1782/1106-1842/1174)  63.4 W50-W56 N08-N12
      5 1708 1113 (N07W44 L= 35)   249     72        0 (1673/1099-1737/1126)  11.3 W41-W46 N06-N08
      6  156 1001 (S11E64 L=287)   671    157        1 ( 127/ 958- 192/1041)  29.7 E67-E59 S13-S08
      7  919  953 (S10E06 L=346)   300    122        0 ( 852/ 932- 979/ 982)  -7.1 E09-E01 S11-S09
      8  596  726 (S26E26 L=325)   991    343        0 ( 480/ 606- 716/ 849)  41.6 E36-E17 S34-S17
      9  417  595 (S36E42 L=309)   328     75        1 ( 386/ 579- 450/ 614)   9.5 E44-E38 S36-S34
     10 1545  512 (S30W44 L= 35)   781    224        1 (1494/ 455-1610/ 578) -23.2 W37-W51 S33-S26
  

    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