Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1207 1832 (N52W52 L= 60)   460     47        1 (1179/1809-1232/1852)  -8.1 W48-W57 N50-N54
      2  304 1284 (N37E51 L=317)   815    180        0 ( 255/1228- 364/1323) -21.8 E58-E43 N32-N39
      3 1050 1290 (N22W08 L= 15)   259    128        0 ( 991/1255-1095/1313) -50.2 W04-W09 N18-N24
      4  935 1131 (N15E03 L=  4)   998    257        2 ( 895/1038- 994/1254)  85.9 E04-E01 N08-N22
      5  427 1124 (N25E38 L=330)   300     77        0 ( 398/1104- 464/1152) -26.6 E40-E34 N22-N26
      6  154  929 (N16E70 L=297)  1142    112        2 ( 103/ 905- 191/ 967)  11.3 E83-E63 N14-N18
      7 1645  891 (S16W36 L= 44)   389    105        1 (1622/ 839-1667/ 938)  51.3 W33-W38 S18-S13
      8  305  837 (N09E53 L=315)   871    167        0 ( 262/ 807- 368/ 892) -11.3 E56-E46 N07-N12
      9 1232  606 (S22W03 L= 11)   314     71        0 (1201/ 571-1261/ 643)  31.0 W00-W05 S24-S20
     10  850  344 (S32E30 L=337)   559     89        1 ( 828/ 301- 873/ 388)  40.6 E33-E26 S34-S28
     11  678  293 (S32E47 L=320)  1109    150        2 ( 643/ 223- 709/ 354)  54.0 E52-E43 S38-S27
  

    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