Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  776 1767 (N49E37 L=293)   495     66        0 ( 740/1747- 818/1784)  26.6 E40-E33 N46-N50
      2 1797 1434 (N36W66 L= 37)  2372    220        2 (1746/1342-1833/1525)  81.3 W59-W73 N29-N42
      3  329 1440 (N19E59 L=272)  1732    204        6 ( 270/1360- 401/1504)  55.0 E64-E52 N13-N23
      4 1434 1360 (N27W24 L=355)  1011    190        1 (1408/1270-1466/1448)  77.0 W22-W27 N20-N33
      5  278 1320 (N10E59 L=271)  1671    257        3 ( 202/1253- 390/1384)  27.9 E66-E49 N05-N15
      6 1262 1199 (N14W13 L=344)  1184    249        3 (1158/1166-1377/1227)  12.1 W06-W20 N11-N17
      7  296 1138 (S00E52 L=278)   431     52        0 ( 266/1122- 336/1154)   0.0 E55-E48 S01-N00
      8  962  724 (S18E01 L=330)  1025    394        1 ( 899/ 570-1026/ 870) -68.7 E05-W04 S27-S09
      9  326  824 (S19E48 L=283)   301     39        0 ( 299/ 810- 353/ 838)  21.8 E50-E45 S20-S18
     10  360  743 (S24E46 L=285)   433     79        2 ( 331/ 713- 387/ 773)  63.4 E48-E42 S26-S22
     11 1624  530 (S24W53 L= 23)  2072    337        1 (1517/ 438-1730/ 665)  62.1 W45-W58 S30-S13
     12 1006  314 (S48W10 L=341)  1788    323        1 ( 956/ 190-1055/ 444) -87.8 W05-W14 S62-S36
  

    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