Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  934 1791 (N38E35 L= 64)   486     67        1 ( 897/1773- 972/1812)   8.1 E38-E31 N36-N39
      2  820 1693 (N28E35 L= 64)   384    122        0 ( 787/1637- 845/1736) -36.9 E39-E31 N24-N30
      3 1104 1572 (N28E12 L= 87)   717    226        1 (1029/1490-1150/1659) -26.6 E20-E06 N23-N32
      4  748 1494 (N13E30 L= 70)   308    118        0 ( 730/1439- 772/1545) -74.1 E31-E28 N09-N16
      5  927 1343 (N09E14 L= 85)   258     94        1 ( 887/1314- 968/1374)  59.0 E16-E12 N06-N11
      6 1433 1342 (N24W15 L=115)   345     77        1 (1394/1312-1458/1367)  36.9 W12-W16 N21-N25
      7 1086 1302 (N11E04 L= 95)   648    263        0 ( 986/1205-1143/1358) -18.4 E11-W01 N07-N13
      8 1638 1093 (N16W35 L=134)   735    175        1 (1576/1036-1699/1151) -15.3 W29-W40 N14-N18
      9 1515 1100 (N12W26 L=126)   242    119        0 (1498/1052-1542/1148)  90.0 W24-W26 N09-N15
     10 1860  716 (N04W71 L=170)  1028    126        1 (1822/ 675-1902/ 773)  48.8 W66-W75 N00-N08
     11 1036  620 (S28W13 L=113)   584    146        0 ( 976/ 594-1109/ 645)  35.5 W09-W17 S30-S24
     12 1534  446 (S21W55 L=154)  1062    268        1 (1456/ 354-1580/ 548) -70.3 W49-W59 S28-S14
     13 1064  449 (S37W23 L=122)   371    106        0 (1026/ 408-1108/ 487)  81.9 W21-W23 S40-S33
     14 1286  455 (S29W37 L=136)   370    128        0 (1229/ 437-1342/ 469)  35.5 W33-W40 S32-S27
     15 1453  441 (S24W49 L=148)   465     69        1 (1416/ 420-1488/ 468)  36.9 W46-W51 S26-S22
     16 1383  332 (S33W53 L=152)   579     83        0 (1341/ 312-1425/ 353)  23.2 W48-W56 S34-S30
  

    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