Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  623 1774 (N41E57 L=322)   661     75        1 ( 589/1754- 660/1795)  31.0 E61-E53 N38-N42
      2  973 1555 (N32E14 L=  5)   793    296        1 ( 905/1432-1049/1665) -41.2 E22-E06 N24-N38
      3  692 1563 (N27E34 L=345)   352    118        1 ( 673/1515- 730/1607)  90.0 E35-E32 N23-N30
      4  339 1461 (N14E60 L=320)   754     96        1 ( 321/1419- 363/1509) -50.2 E63-E55 N11-N17
      5 1815 1410 (N39W65 L= 85)   627     70        0 (1798/1378-1826/1450)  28.6 W60-W72 N36-N42
      6 1292 1301 (N21W12 L= 32)   432     77        2 (1264/1266-1322/1338)  78.7 W10-W14 N18-N23
      7 1736 1311 (N30W49 L= 68)   453    114        0 (1682/1290-1787/1336)  -5.2 W43-W55 N28-N30
      8  575 1137 (S02E30 L=350)   321     63        1 ( 541/1117- 607/1154)  11.3 E32-E27 S02-N00
      9  305  961 (S17E49 L=330)   516    202        0 ( 279/ 875- 351/1025)  73.3 E51-E45 S21-S11
     10  409  932 (S17E39 L=340)   643    237        1 ( 357/ 849- 439/ 993) -58.0 E44-E36 S21-S12
     11 1692  901 (N04W47 L= 66)   461    102        3 (1660/ 853-1727/ 942)  71.6 W44-W49 N00-N06
     12 1175  855 (S08W12 L= 32)   930    266        2 (1077/ 783-1267/ 936)  53.1 W07-W16 S13-S01
     13  288  785 (S28E51 L=328)   625    155        0 ( 268/ 739- 338/ 849) -39.8 E53-E46 S30-S24
     14  767  609 (S31E10 L= 10)  1759    396        7 ( 624/ 558- 934/ 696)  35.0 E20-W00 S37-S22
     15 1644  570 (S17W54 L= 74)  1818    372        2 (1479/ 507-1797/ 616)  23.2 W41-W70 S23-S11
     16  406  562 (S42E42 L=337)  1712    239        5 ( 305/ 503- 514/ 605)  20.6 E62-E30 S48-S36
     17 1284  463 (S31W30 L= 50)   559    133        0 (1222/ 447-1351/ 483)  24.0 W25-W34 S32-S28
     18 1143  404 (S38W23 L= 42)   284     52        1 (1120/ 374-1167/ 433)  90.0 W21-W24 S40-S35
     19 1053  361 (S43W18 L= 38)   355     37        1 (1028/ 340-1084/ 378)  36.9 W15-W20 S44-S41
  

    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