Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  774 1284 (N19E17 L=220)   483     65        1 ( 739/1257- 813/1320) -11.3 E18-E13 N17-N21
      2 1511 1255 (N17W33 L=269)   539    196        1 (1457/1176-1541/1315)  71.6 W29-W35 N12-N21
      3 1339 1279 (N19W21 L=257)   300     60        1 (1298/1259-1373/1299)   9.5 W17-W23 N18-N20
      4 1840 1209 (N13W63 L=299)   433    118        0 (1822/1170-1862/1240) -68.2 W60-W65 N10-N15
      5 1385 1153 (N11W23 L=260)   432    123        2 (1332/1114-1445/1198) -23.2 W19-W27 N08-N14
      6  484 1170 (N12E36 L=200)   257     95        0 ( 443/1152- 525/1186) -15.9 E39-E32 N10-N12
      7  516 1047 (N04E33 L=203)   230     81        0 ( 505/1002- 534/1090)  78.7 E33-E31 N01-N06
      8 1578 1033 (N03W36 L=273)   239     83        0 (1549/1006-1613/1063)  31.0 W34-W39 N01-N05
      9  752  978 (N00E17 L=220)   311     76        0 ( 737/ 931- 769/1021) -90.0 E17-E15 S02-N03
     10  336  960 (S02E48 L=189)  1437    272        4 ( 259/ 914- 409/1011)   0.0 E55-E41 S04-N01
     11 1402  850 (S07W24 L=261)   259     82        0 (1383/ 818-1431/ 881)  53.1 W22-W26 S09-S05
     12 1730  781 (S13W50 L=287)  1078    153        1 (1696/ 711-1758/ 865)  84.3 W47-W52 S17-S07
     13  624  671 (S19E27 L=210)   461    101        1 ( 586/ 618- 654/ 716) -80.5 E29-E24 S22-S16
     14 1374  575 (S26W25 L=261)   268     87        0 (1347/ 527-1393/ 614)  71.6 W23-W25 S28-S22
     15 1078  480 (S32W04 L=241)   798    310        0 (1022/ 379-1131/ 611)  68.7 W00-W07 S40-S22
     16 1607  578 (S26W43 L=280)   458     74        1 (1573/ 556-1639/ 602)  11.3 W40-W45 S27-S24
  

    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