Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  730 1813 (N42E42 L= 92)  1410    124        2 ( 679/1764- 778/1853) -14.0 E51-E35 N37-N45
      2 1151 1664 (N37E02 L=132)   243     39        0 (1134/1636-1163/1689) -76.0 E03-E01 N34-N38
      3  480 1643 (N26E52 L= 83)   611     97        0 ( 444/1607- 512/1673) -20.0 E57-E46 N23-N27
      4 1521 1632 (N42W28 L=163)   713     74        2 (1487/1601-1559/1661)  20.6 W24-W32 N39-N45
      5 1663 1303 (N22W38 L=173)  1288    405        0 (1555/1163-1768/1454) -36.4 W29-W48 N14-N29
      6 1135 1333 (N14W02 L=137)   199     56        0 (1109/1314-1165/1345)  33.7 W00-W03 N12-N14
      7  592 1060 (S10E27 L=108)   267     50        1 ( 566/1037- 629/1079)  26.6 E28-E24 S11-S08
      8  691 1047 (S09E20 L=114)   219     68        0 ( 662/1025- 728/1073) -11.3 E22-E17 S10-S07
      9  119 1000 (S17E72 L= 62)  5113    251        9 (  76/ 890- 191/1071) -79.7 E85-E61 S23-S12
     10 1540 1009 (N01W32 L=166)   461    162        0 (1470/ 987-1607/1036)  11.3 W26-W36 N00-N03
     11 1848  951 (N05W59 L=193)   516     83        1 (1824/ 907-1867/ 986) -33.7 W55-W62 N02-N06
     12  368  899 (S22E43 L= 92)   254    100        1 ( 315/ 871- 404/ 921) -14.0 E47-E39 S23-S20
     13  857  875 (S17E08 L=127)   252     53        0 ( 821/ 857- 887/ 890)  26.6 E10-E06 S18-S16
     14  167  740 (S33E71 L= 63)   635     72        0 ( 153/ 699- 179/ 777)  59.0 E73-E68 S35-S30
     15 1499  685 (S18W37 L=171)   359     83        1 (1460/ 656-1530/ 720)  53.1 W34-W38 S20-S15
     16 1541  579 (S23W44 L=178)   547     85        2 (1500/ 553-1586/ 607)  -6.3 W39-W48 S24-S22
     17  584  428 (S51E28 L=106)   702     83        0 ( 559/ 377- 607/ 468) -74.1 E30-E26 S55-S48
     18 1508  444 (S32W49 L=183)   349     78        1 (1482/ 413-1532/ 467) -63.4 W45-W50 S34-S30
     19 1120  363 (S46W22 L=157)   547    118        1 (1098/ 300-1145/ 424) -52.1 W18-W26 S51-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