Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  294 1308 (N17E55 L=260)   958    200        3 ( 253/1246- 324/1370) -54.5 E59-E51 N13-N20
      2  324 1234 (N13E51 L=265)  1257    318        4 ( 237/1186- 418/1297)  25.0 E59-E43 N09-N17
      3  417 1166 (N09E42 L=274)   261     71        2 ( 395/1133- 441/1187)  53.1 E43-E39 N06-N10
      4  144  938 (S07E70 L=245)   683    105        1 ( 128/ 886- 158/ 984) -63.4 E73-E67 S10-S04
      5 1516  815 (S09W33 L=348)  1463    339        4 (1410/ 687-1611/ 959)  59.5 W26-W39 S17-N00
      6  274  925 (S07E53 L=262)   550    157        2 ( 230/ 891- 310/ 948)  20.6 E58-E49 S08-S05
      7  279  754 (S18E54 L=261)  1716    356        5 ( 228/ 609- 338/ 867) -83.3 E59-E49 S27-S10
      8  513  660 (S22E34 L=281)   866    271        1 ( 455/ 593- 550/ 796)  57.3 E39-E30 S27-S13
      9  376  716 (S20E45 L=270)   581    111        2 ( 340/ 668- 409/ 767) -56.3 E47-E42 S22-S16
  

    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