Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  715 1792 (N40E51 L=137)  1727    185        4 ( 652/1726- 788/1854) -12.7 E69-E38 N36-N43
      2  775 1627 (N29E33 L=155)  1031    214        2 ( 662/1583- 856/1670)  48.0 E39-E28 N23-N33
      3 1604 1515 (N43W33 L=221)  1986    270        3 (1530/1399-1654/1603) -60.3 W26-W39 N34-N48
      4  444 1475 (N12E50 L=137)   719    114        0 ( 415/1424- 483/1522)  81.9 E51-E48 N08-N15
      5 1167 1395 (N23W00 L=188)  1011    328        4 (1041/1308-1242/1505) -23.6 E10-W06 N18-N27
      6 1452 1460 (N35W19 L=206)   423    161        1 (1389/1415-1521/1487)  -4.8 W13-W25 N33-N35
      7 1586 1269 (N25W32 L=219)  2212    516        4 (1511/1024-1675/1458) -54.7 W23-W41 N12-N36
      8 1062 1295 (N14E04 L=184)   355    130        1 (1013/1258-1118/1334)  -8.1 E07-W00 N12-N15
      9 1850 1173 (N26W60 L=247)  1335    128        2 (1819/1112-1872/1232) -76.0 W55-W62 N22-N30
     10  328 1070 (S14E48 L=140)  1429    315        5 ( 248/ 994- 402/1131)  23.2 E55-E41 S18-S09
     11  213  925 (S25E59 L=128)   580     82        0 ( 194/ 876- 228/ 971)  74.1 E62-E56 S28-S21
     12 1636  925 (N06W41 L=228)   291     43        1 (1615/ 904-1666/ 944)  11.3 W38-W43 N04-N06
     13  551  887 (S20E27 L=161)   411    120        2 ( 508/ 853- 589/ 939)  60.3 E29-E24 S23-S16
     14 1486  709 (S10W36 L=224)   317     55        0 (1463/ 684-1508/ 734) -14.0 W34-W38 S11-S08
     15 1462  644 (S14W37 L=224)   444    116        1 (1435/ 606-1492/ 691) -23.2 W33-W40 S16-S11
     16 1363  615 (S18W31 L=219)   509    131        0 (1297/ 596-1404/ 655)  45.0 W27-W34 S20-S15
     17  391  602 (S43E41 L=146)  1341    175        3 ( 306/ 572- 460/ 630)   9.0 E53-E33 S45-S39
     18 1542  497 (S21W50 L=238)   393    106        0 (1493/ 468-1587/ 519)  39.8 W47-W53 S24-S18
  

    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