Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1162 1621 (N33E02 L=341)   375     91        0 (1124/1603-1212/1635)  36.9 E04-W00 N31-N34
      2 1040 1499 (N22E08 L=336)   279     80        1 (1010/1472-1073/1524) -21.8 E09-E04 N20-N23
      3  978 1423 (N16E10 L=334)   220     85        1 ( 953/1405-1017/1456) -14.0 E11-E07 N15-N17
      4 1209 1397 (N19W05 L=348)   276     70        0 (1176/1382-1242/1421)  36.9 W02-W06 N17-N20
      5 1799 1102 (N14W51 L= 34)   471     81        1 (1780/1061-1819/1147) -59.0 W49-W53 N11-N16
      6  275  962 (S20E51 L=292)  1002    255        1 ( 218/ 895- 338/1048)  45.0 E57-E45 S24-S14
      7 1162  678 (S24W15 L=358)   262     91        0 (1125/ 647-1202/ 710)  59.0 W13-W17 S26-S21
      8 1588  622 (S18W46 L= 29)   657    162        1 (1510/ 583-1650/ 663)  36.9 W40-W49 S22-S14
      9  960  588 (S34W04 L=347)   346    143        0 ( 894/ 572-1034/ 609)  26.6 W00-W08 S35-S30
     10  707  588 (S38E14 L=329)   228     54        0 ( 684/ 565- 735/ 608)  45.0 E15-E12 S40-S36
     11  239  573 (S46E70 L=274)   703     68        0 ( 218/ 549- 268/ 602) -26.6 E74-E65 S47-S43
     12  933  279 (S57W13 L=356)   667    109        1 ( 871/ 255- 977/ 309)   4.4 W06-W19 S58-S54
  

    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