Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1510 1688 (N32W56 L=330)   505     66        1 (1477/1674-1543/1705) -36.9 W53-W58 N30-N33
      2 1679 1419 (N13W54 L=329)   410     90        1 (1642/1388-1717/1444) -38.7 W51-W57 N10-N14
      3 1726 1285 (N04W53 L=327)   402     52        1 (1699/1263-1754/1303) -11.3 W50-W55 N03-N05
      4  432 1102 (N23E37 L=237)  2359    603        3 ( 242/ 995- 654/1222) -32.8 E53-E22 N13-N33
      5  180  970 (N18E63 L=212)   541     78        2 ( 156/ 948- 208/1001) -36.9 E65-E59 N17-N20
      6  358  727 (N02E51 L=224)   296     89        1 ( 340/ 681- 374/ 771)  76.0 E51-E49 N00-N04
      7  573  329 (S27E51 L=224)   413     59        1 ( 556/ 307- 592/ 353)  53.1 E53-E48 S29-S25
      8  649  279 (S32E48 L=226)   637    116        1 ( 620/ 251- 672/ 316)  35.5 E51-E44 S34-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