Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  606 1759 (N67E34 L=196)   553     61        0 ( 580/1738- 639/1776)   0.0 E40-E26 N65-N69
      2  954 1611 (N45W08 L=238)   615    100        1 ( 905/1579-1000/1635) -26.6 W04-W12 N42-N48
      3 1742 1544 (N21W73 L=303)   893     95        1 (1722/1506-1767/1576)  33.7 W68-W81 N18-N22
      4 1608 1406 (N17W48 L=278)  1793    446        0 (1426/1344-1797/1457) -25.8 W35-W64 N09-N24
      5 1315 1318 (N19W24 L=254)   283    124        0 (1297/1267-1334/1368)  74.1 W21-W25 N15-N22
      6 1030 1317 (N24W06 L=236)   202     90        0 ( 992/1307-1081/1336) -18.4 W03-W09 N22-N25
      7  497 1181 (N24E33 L=197)  1545    339        1 ( 416/1035- 539/1298) -66.0 E40-E29 N14-N32
      8  655 1199 (N23E21 L=209)   302    101        0 ( 613/1174- 704/1233) -39.8 E23-E17 N20-N25
      9 1306 1152 (N09W19 L=249)   380    192        1 (1237/1117-1391/1173)  -5.7 W14-W24 N07-N10
     10 1610 1133 (N02W39 L=269)   266     71        2 (1593/1095-1625/1164)  56.3 W37-W41 N00-N03
     11  302  925 (N10E52 L=178)   457    109        2 ( 276/ 890- 348/ 958)  15.9 E54-E47 N08-N12
     12  394  683 (S05E49 L=181)   289     62        1 ( 374/ 660- 423/ 708) -45.0 E50-E46 S06-S03
     13  635  618 (S12E32 L=198)   360    111        0 ( 611/ 590- 679/ 662) -53.1 E33-E29 S14-S09
     14  517  498 (S19E46 L=184)   343     72        0 ( 493/ 465- 545/ 527)  23.2 E49-E42 S20-S16
     15  888  396 (S31E22 L=208)  1205    348        1 ( 769/ 305- 996/ 511)  24.6 E35-E11 S37-S24
     16  665  376 (S29E40 L=190)   321     84        1 ( 638/ 334- 693/ 411)  32.0 E44-E36 S32-S26
  

    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