Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1292 1479 (N28W28 L=183)   525    103        1 (1251/1445-1344/1512)   7.1 W25-W33 N26-N28
      2  443 1332 (N37E36 L=119)  1031    276        0 ( 382/1241- 509/1427) -51.3 E42-E30 N29-N44
      3  212 1219 (N32E57 L= 98)   545    112        2 ( 187/1199- 245/1242)   0.0 E60-E53 N30-N33
      4  230 1079 (N23E56 L= 99)  1104    230        1 ( 210/ 953- 268/1163)  81.3 E58-E51 N15-N28
      5  626 1143 (N21E23 L=132)   240     69        1 ( 595/1125- 661/1162) -36.9 E24-E20 N20-N23
      6 1698 1118 (S03W45 L=200)   990    143        2 (1642/1080-1753/1152) -18.4 W40-W49 S05-N00
      7 1563 1074 (S02W34 L=189)   391    128        0 (1531/ 993-1598/1131)  45.0 W30-W37 S06-N00
      8  887  907 (N03E10 L=145)  1101    217        2 ( 822/ 813- 959/1002) -76.0 E12-E06 S03-N09
      9 1567  844 (S16W30 L=185)   396     72        1 (1554/ 799-1581/ 889)  80.5 W28-W31 S18-S12
     10  873  577 (S17E19 L=136)   899    198        4 ( 798/ 490- 915/ 628) -78.7 E22-E16 S22-S12
     11  711  282 (S35E43 L=112)  1043    198        2 ( 664/ 236- 754/ 329)  18.4 E49-E37 S38-S31
  

    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