Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1114 1724 (N52W08 L=211)   321     44        1 (1087/1699-1142/1742)   0.0 W05-W10 N49-N53
      2  543 1519 (N34E40 L=163)  1813    406        2 ( 395/1445- 719/1604) -13.0 E61-E22 N29-N39
      3  221 1250 (N15E64 L=139)   608    102        1 ( 206/1197- 242/1294) -80.5 E65-E62 N11-N17
      4  726  688 (S18E19 L=184)   495    207        0 ( 633/ 665- 819/ 720)  -8.7 E25-E12 S19-S16
      5  391  584 (S27E49 L=154)  2314    314        4 ( 314/ 430- 460/ 690)  28.9 E77-E39 S40-S19
      6 1702  646 (S21W52 L=255)   408     32        2 (1683/ 627-1725/ 668) -36.9 W49-W54 S22-S19
      7  784  611 (S23E16 L=187)   396    154        0 ( 740/ 554- 845/ 647)  41.2 E19-E11 S27-S20
      8 1162  495 (S31W11 L=213)   246     81        0 (1122/ 477-1201/ 513)  26.6 W07-W13 S32-S29
      9 1562  417 (S38W48 L=251)   718    103        0 (1535/ 376-1587/ 454) -49.4 W43-W51 S41-S34
     10 1443  362 (S42W39 L=241)   980     96        1 (1404/ 336-1494/ 396)  14.0 W34-W44 S44-S39
     11 1212  308 (S47W18 L=221)  1211    199        2 (1131/ 285-1303/ 346)   7.1 W10-W27 S48-S43
  

    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