Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  809 1854 (N47E46 L=340)   938    152        2 ( 776/1815- 867/1882)   0.0 E51-E39 N43-N50
      2  594 1453 (N14E37 L=349)  1520    403        3 ( 505/1342- 721/1577) -21.4 E48-E25 N09-N19
      3 1560 1319 (N28W29 L= 55)   610    197        0 (1492/1244-1615/1422) -36.0 W22-W33 N24-N32
      4 1215 1285 (N17W06 L= 32)   368     78        1 (1178/1252-1252/1321) -31.0 W03-W08 N15-N18
      5  589 1147 (S04E29 L=358)   585    157        1 ( 518/1097- 664/1186)  45.0 E32-E24 S08-N00
      6 1235 1130 (N08W10 L= 36)   266     22        0 (1215/1109-1259/1152) -33.7 W08-W11 N07-N09
      7  847  962 (S09E09 L= 17)   669    238        1 ( 736/ 932- 965/1000)  11.3 E16-E01 S11-S07
      8  206  797 (S33E57 L=330)   870    217        0 ( 169/ 750- 249/ 844)  35.0 E61-E51 S36-S29
      9 1295  536 (S25W29 L= 56)   258     53        0 (1263/ 523-1326/ 548)   0.0 W26-W31 S25-S23
     10  876  513 (S37W03 L= 29)   468    156        0 ( 802/ 491- 958/ 537)  29.1 E02-W08 S39-S33
  

    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