Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1621 1507 (N38W44 L=131)  1397    255        2 (1532/1403-1695/1596) -50.2 W37-W48 N30-N44
      2  931 1508 (N32E11 L= 76)   215     49        0 ( 910/1491- 954/1525) -18.4 E12-E09 N30-N32
      3 1076 1407 (N26W00 L= 87)   276     47        2 (1056/1386-1107/1429) -18.4 E01-W02 N24-N27
      4  186 1319 (N12E70 L= 17)   765     76        2 ( 167/1291- 202/1341) -45.0 E74-E67 N10-N13
      5  818  812 (S13E11 L= 76)   302     94        0 ( 803/ 760- 834/ 857) -71.6 E12-E10 S16-S10
      6 1624  719 (S13W44 L=131)   307     62        0 (1603/ 692-1656/ 738)  45.0 W42-W46 S14-S11
      7 1704  674 (S15W52 L=139)   486     88        2 (1680/ 640-1740/ 704)  45.0 W49-W55 S17-S13
      8  429  639 (S28E40 L= 47)   305    100        0 ( 396/ 620- 471/ 670)  -9.5 E42-E36 S29-S26
      9 1728  586 (S21W58 L=145)   645    127        4 (1691/ 558-1760/ 615)  26.6 W54-W61 S22-S18
     10 1251  473 (S32W22 L=109)   238     47        1 (1229/ 455-1277/ 492) -14.0 W19-W23 S33-S30
  

    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