Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  842 1813 (N49E40 L=255)   503     68        0 ( 814/1782- 879/1832)  78.7 E42-E36 N45-N51
      2  648 1703 (N36E46 L=248)   907    177        1 ( 609/1641- 710/1771) -20.9 E57-E36 N32-N40
      3  952 1686 (N41E21 L=274)   365     86        0 ( 923/1658- 990/1716) -20.6 E24-E16 N39-N42
      4  433 1644 (N28E63 L=231)   482     32        2 ( 414/1627- 457/1664)  -6.3 E67-E58 N26-N29
      5  849 1573 (N30E23 L=271)   249    102        0 ( 820/1535- 897/1597) -20.6 E26-E18 N28-N31
      6 1548 1560 (N44W32 L=327)   325     34        0 (1530/1537-1565/1581)  53.1 W30-W34 N42-N46
      7 1044 1374 (N21E05 L=289)   295     92        1 (1006/1350-1089/1405)  -9.5 E07-E01 N19-N22
      8  279 1196 (S03E55 L=239)   435     39        0 ( 259/1178- 309/1218)   0.0 E57-E51 S03-S01
      9 1755  658 (S10W61 L=355)   801    105        1 (1714/ 628-1805/ 683)  29.1 W55-W67 S12-S07
     10  899  575 (S30W00 L=294)   332     83        1 ( 858/ 556- 943/ 593)  26.6 E02-W02 S32-S29
     11  954  383 (S43W09 L=303)   457    119        0 ( 917/ 333- 995/ 422)  77.5 W07-W11 S48-S39
     12 1549  348 (S33W62 L=356)   620     78        0 (1514/ 331-1588/ 369)   4.8 W56-W68 S34-S32
  

    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