Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  298 1304 (N15E51 L=144)  1987    330        3 ( 261/1168- 366/1447) -71.6 E57-E45 N06-N24
      2 1717 1415 (N22W51 L=247)   369     34        0 (1695/1391-1741/1439) -56.3 W49-W53 N20-N23
      3  670 1314 (N15E22 L=173)   623    193        3 ( 598/1290- 768/1339)  -4.8 E27-E15 N13-N16
      4  217 1208 (N10E59 L=137)   644    121        1 ( 182/1149- 245/1260)  50.2 E61-E56 N06-N12
      5  332 1177 (N07E47 L=149)   667    154        2 ( 288/1102- 378/1226)  49.4 E50-E42 N02-N10
      6  891  835 (S14E08 L=187)   938    316        0 ( 825/ 701- 927/ 979)  70.6 E13-E05 S22-S05
      7 1367  820 (S15W22 L=217)   265     84        1 (1339/ 770-1387/ 848) -51.3 W19-W23 S18-S13
      8 1446  792 (S17W27 L=223)   281     99        0 (1397/ 775-1492/ 816)   0.0 W23-W30 S18-S15
      9  890  471 (S38E11 L=185)   756    211        1 ( 782/ 450- 989/ 504)   6.7 E19-E02 S39-S35
     10 1686  443 (S39W62 L=258)   555     70        1 (1651/ 423-1722/ 463)  -4.1 W57-W71 S40-S37
     11  901  310 (S51E12 L=183)   290     44        2 ( 874/ 290- 928/ 337)   9.5 E15-E09 S52-S48
  

    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