Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  588 1628 (N19E45 L=182)   678    152        1 ( 556/1573- 628/1697) -30.3 E51-E39 N15-N23
      2  318 1514 (N08E63 L=164)  1757    128        7 ( 285/1416- 349/1619) -54.5 E72-E56 N01-N15
      3  503 1463 (N07E44 L=183)   254     52        1 ( 485/1435- 527/1485) -21.8 E45-E40 N05-N08
      4  369 1270 (S07E48 L=179)   290     71        0 ( 330/1248- 405/1291)   0.0 E51-E44 S07-S05
      5  316 1214 (S11E51 L=176)   348     55        0 ( 291/1194- 346/1231) -11.3 E53-E47 S11-S09
      6  890 1161 (S03E11 L=216)   230     49        2 ( 862/1138- 920/1179)   0.0 E12-E08 S03-S01
      7  609 1059 (S14E25 L=202)   317    105        1 ( 578/1013- 630/1114) -80.5 E27-E24 S17-S11
      8  429  993 (S22E37 L=190)   805    261        2 ( 341/ 916- 522/1079) -20.6 E45-E29 S24-S17
      9  221  896 (S31E58 L=169)   454     63        1 ( 200/ 867- 250/ 920) -56.3 E60-E54 S32-S29
     10  204  792 (S38E63 L=164)   526     33        0 ( 174/ 778- 234/ 808)   5.2 E68-E57 S38-S36
     11 1316  424 (S33W39 L=266)   305     69        1 (1280/ 411-1359/ 441)  18.4 W36-W42 S34-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