Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  637 1587 (N17E39 L=185)   263     63        0 ( 617/1564- 666/1617) -21.8 E41-E36 N15-N18
      2  347 1556 (N11E62 L=162)   772     69        1 ( 321/1519- 373/1590) -26.6 E66-E58 N08-N12
      3  535 1451 (N07E40 L=184)   306     81        0 ( 514/1414- 568/1484) -23.2 E43-E36 N05-N08
      4  403 1260 (S07E44 L=180)   325     85        0 ( 359/1238- 442/1280)   0.0 E48-E41 S07-S05
      5  324 1201 (S11E50 L=175)   470    111        0 ( 286/1164- 372/1223)  23.2 E52-E45 S13-S09
      6  932 1150 (S02E08 L=216)   183     37        0 ( 906/1130- 953/1165)   0.0 E09-E06 S03-S01
      7  391 1046 (S19E41 L=183)   287     29        1 ( 371/1018- 412/1068) -71.6 E42-E39 S20-S17
      8  289  983 (S25E50 L=175)   338     76        0 ( 270/ 959- 318/1019) -31.0 E51-E46 S25-S22
      9  487  953 (S23E32 L=192)   323    100        0 ( 464/ 913- 522/ 996) -38.7 E34-E28 S24-S20
     10  232  879 (S32E56 L=168)   348     41        1 ( 214/ 850- 254/ 896) -71.6 E58-E53 S33-S30
     11  216  786 (S38E60 L=164)   609     81        0 ( 178/ 771- 249/ 811)   9.5 E66-E54 S38-S35
  

    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