Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1556 1700 (N30W61 L= 72)   553     70        1 (1533/1677-1577/1723)  21.8 W57-W63 N28-N31
      2 1843 1048 (S14W56 L= 68)  1409    170        1 (1793/ 997-1877/1100)  26.6 W50-W60 S16-S09
      3  208  935 (N17E60 L=312)   427     86        0 ( 186/ 902- 235/ 973) -68.2 E61-E57 N14-N19
      4 1201  834 (S08W06 L= 17)   408    119        0 (1144/ 794-1262/ 857)   0.0 W01-W09 S09-S06
      5  490  719 (N00E40 L=331)   769    187        1 ( 470/ 630- 510/ 809)  68.2 E42-E37 S05-N05
      6 1372  658 (S23W12 L= 24)   717    224        1 (1325/ 551-1396/ 760)  56.3 W07-W15 S29-S16
      7 1021  713 (S11E07 L=  4)   224     47        0 (1007/ 683-1037/ 746)  63.4 E08-E06 S12-S08
      8 1489  609 (S30W20 L= 31)   632    120        1 (1444/ 552-1538/ 670)  36.9 W15-W23 S32-S26
      9  994  434 (S28E17 L=354)   539    198        0 ( 913/ 375-1075/ 492)  15.9 E24-E10 S30-S25
  

    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