Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1381 1389 (N21W31 L=196)   569    158        0 (1302/1367-1450/1411) -24.0 W26-W35 N18-N23
      2 1614  996 (S07W37 L=202)   633    180        0 (1574/ 922-1658/1077) -74.7 W34-W39 S12-S01
      3 1715 1032 (S07W45 L=210)   280     56        0 (1704/1001-1727/1067)  76.0 W44-W46 S09-S05
      4 1825  932 (S16W55 L=220)   990    130        1 (1778/ 879-1853/ 987)  71.6 W49-W58 S19-S13
      5 1484  775 (S17W24 L=189)   553    174        1 (1445/ 703-1502/ 865)  60.9 W20-W26 S21-S11
      6 1086  825 (S06W00 L=165)   239    137        1 (1034/ 777-1135/ 864)  26.6 E03-W03 S08-S04
      7 1261  698 (S17W08 L=173)   430    182        1 (1189/ 665-1319/ 759)  18.4 W03-W12 S18-S14
      8 1686  583 (S36W42 L=207)  1738    328        1 (1637/ 490-1737/ 720)  68.2 W38-W52 S43-S27
      9  365  555 (S12E57 L=108)   969    235        1 ( 336/ 450- 408/ 634)  69.0 E62-E53 S19-S06
     10 1204  432 (S35E00 L=165)  2966    524        9 (1043/ 276-1330/ 600)  35.1 E16-W11 S44-S24
     11  595  551 (S15E37 L=127)   436    151        1 ( 522/ 533- 667/ 562) -11.3 E42-E32 S16-S13
     12  996  261 (S44E22 L=143)   514     51        1 ( 959/ 234-1025/ 292) -90.0 E23-E19 S47-S40
  

    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