Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1283 1643 (N45W25 L=345)  1000    196        2 (1196/1603-1375/1684) -19.7 W17-W31 N41-N48
      2 1754 1203 (N12W54 L= 14)  1390    213        2 (1677/1138-1818/1288) -42.3 W47-W59 N07-N17
      3 1587 1156 (N10W38 L=359)   476    107        0 (1568/1105-1603/1206) -80.5 W36-W39 N07-N13
      4  877 1137 (N11E09 L=312)   572    139        2 ( 843/1084- 925/1196) -60.3 E11-E06 N07-N14
      5 1637  999 (N00W41 L=  2)   319     72        0 (1598/ 982-1677/1010)   0.0 W37-W44 N00-N01
      6 1562  763 (S14W36 L=356)   474    143        1 (1542/ 692-1580/ 810) -82.9 W34-W37 S18-S10
      7 1637  793 (S12W42 L=  2)   274     34        2 (1615/ 777-1662/ 810)   0.0 W39-W43 S13-S11
      8 1448  740 (S15W27 L=348)   303    108        1 (1423/ 712-1491/ 773) -31.0 W25-W30 S16-S12
      9  445  586 (S24E44 L=277)   427     78        0 ( 424/ 548- 478/ 615)  39.8 E47-E40 S26-S21
  

    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