Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  783 1685 (N33E28 L=308)   676    214        0 ( 703/1603- 845/1737) -26.6 E36-E20 N27-N36
      2  588 1661 (N29E42 L=294)   422     86        2 ( 572/1619- 606/1710) -60.3 E44-E39 N26-N33
      3  661 1465 (N17E29 L=307)   349     78        1 ( 637/1417- 676/1503) -71.6 E31-E27 N13-N19
      4 1401 1435 (N24W19 L=355)   329    120        0 (1374/1386-1421/1490)  82.9 W17-W20 N20-N28
      5 1267 1330 (N16W11 L=347)   403    124        0 (1227/1290-1318/1375) -23.2 W07-W14 N13-N17
      6 1739 1290 (N21W46 L= 22)   306     99        0 (1709/1254-1779/1328) -33.7 W43-W49 N19-N23
      7 1632 1296 (N19W36 L= 13)   445     63        1 (1594/1270-1658/1321) -11.3 W33-W38 N18-N21
      8  597  901 (S18E25 L=311)   490    128        2 ( 562/ 829- 637/ 951)  69.4 E27-E22 S22-S14
      9  204  852 (S24E61 L=276)   409     66        0 ( 175/ 820- 228/ 885) -33.7 E64-E58 S25-S21
     10  342  844 (S24E46 L=291)   510    150        1 ( 289/ 820- 400/ 883)   0.0 E50-E40 S24-S21
     11  838  762 (S24E08 L=328)   600    193        2 ( 783/ 684- 874/ 852)  74.7 E11-E06 S29-S18
     12 1715  648 (S17W55 L= 32)   323     59        1 (1685/ 629-1738/ 670)  56.3 W52-W57 S18-S15
  

    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