Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1388 1719 (N46W20 L=259)  1608    209        3 (1295/1662-1482/1755)  -6.7 W11-W28 N42-N49
      2  289 1521 (N19E64 L=175)  2618    156        2 ( 254/1438- 326/1601) -53.1 E71-E59 N13-N25
      3 1223 1411 (N20W08 L=247)  1232    343        1 (1090/1300-1341/1506) -30.5 E01-W16 N14-N25
      4 1767 1387 (N28W51 L=290)  1104    168        3 (1718/1319-1811/1457) -53.1 W47-W55 N23-N31
      5 1858 1269 (N22W60 L=299)  1187    210        0 (1841/1169-1875/1369)  72.9 W56-W65 N15-N28
      6 1644 1312 (N20W38 L=277)   720    127        0 (1605/1274-1693/1365) -45.0 W34-W42 N18-N23
      7 1762 1307 (N22W49 L=288)   561     84        1 (1714/1285-1796/1324)  12.5 W44-W53 N20-N23
      8  137 1034 (S12E69 L=170)   587     68        2 ( 120/1004- 154/1072)  76.0 E71-E65 S13-S09
      9  333 1019 (S12E46 L=192)   412     98        1 ( 303/ 974- 370/1062)  45.0 E48-E43 S14-S09
     10  766  891 (S16E14 L=224)   206    110        0 ( 719/ 868- 814/ 910)  26.6 E17-E11 S18-S14
     11  223  843 (S23E59 L=180)   459     85        0 ( 183/ 825- 260/ 863)   0.0 E64-E54 S24-S21
     12 1387  596 (S26W31 L=269)  2109    488        3 (1211/ 464-1566/ 748)  52.3 W21-W39 S37-S14
     13  491  696 (S31E34 L=205)   543    101        0 ( 454/ 653- 529/ 737)  40.6 E37-E30 S33-S27
     14 1704  600 (S21W56 L=295)   538    118        0 (1675/ 553-1735/ 643)  80.5 W54-W59 S23-S17
     15  392  524 (S43E49 L=190)   336     38        1 ( 373/ 495- 410/ 547)  38.7 E52-E46 S45-S41
     16 1327  351 (S44W38 L=277)   355     90        0 (1307/ 324-1363/ 383) -20.6 W34-W42 S45-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