Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  368 1594 (N41E57 L=332)   528     72        1 ( 341/1573- 397/1609)  26.6 E59-E54 N39-N42
      2 1684 1550 (N35W57 L= 86)   909    116        1 (1645/1507-1710/1579) -78.7 W53-W59 N31-N36
      3  541 1480 (N34E36 L=354)  1001    226        2 ( 459/1429- 628/1521) -11.3 E43-E28 N30-N36
      4  378 1283 (N20E46 L=344)  2944    674        5 ( 235/1080- 483/1434) -39.2 E64-E35 N07-N29
      5  961 1385 (N27E03 L= 26)   435    152        0 ( 891/1359-1035/1400) -12.5 E08-W01 N24-N27
      6 1169 1183 (N13W09 L= 39)   253     66        0 (1128/1159-1206/1208) -21.8 W06-W11 N11-N15
      7 1565  989 (N00W35 L= 65)   373    107        1 (1524/ 957-1608/1013)   9.5 W32-W38 S01-N01
      8  500  911 (S03E34 L=355)   403     64        1 ( 469/ 887- 538/ 931)  11.3 E36-E31 S03-S01
      9  799  830 (S08E14 L= 15)   432    128        0 ( 733/ 787- 852/ 863) -35.5 E18-E11 S10-S05
     10  380  677 (S18E47 L=343)  1144    339        2 ( 312/ 582- 434/ 779)  45.0 E56-E40 S24-S10
     11 1442  580 (S26W28 L= 58)   587    126        1 (1379/ 566-1504/ 595)   0.0 W23-W33 S26-S24
     12  586  519 (S28E33 L=356)   955    154        2 ( 536/ 481- 647/ 556)  21.8 E37-E27 S31-S25
  

    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