Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  564 1786 (N32E66 L= 75)  1606    113        1 ( 538/1722- 593/1842) -27.8 E78-E58 N26-N36
      2 1503 1699 (N56W18 L=159)  1437    115        3 (1453/1654-1536/1747) -66.0 W12-W23 N50-N62
      3 1213 1657 (N40E07 L=134)   548    117        1 (1153/1623-1256/1676)   0.0 E11-E02 N39-N41
      4  908 1470 (N19E19 L=122)  1273    284        1 ( 850/1363- 972/1591) -34.5 E27-E11 N13-N25
      5 1619 1571 (N49W31 L=172)   439     43        0 (1597/1549-1647/1587)  33.7 W28-W35 N46-N51
      6  360 1532 (N10E62 L= 79)  1331    169        3 ( 301/1481- 421/1581)  53.1 E66-E56 N06-N14
      7  405 1487 (N08E54 L= 87)   926    146        2 ( 360/1447- 438/1547)  81.9 E57-E51 N05-N12
      8 1337 1526 (N35W06 L=147)   359     65        1 (1307/1504-1369/1546) -11.3 W03-W08 N32-N36
      9 1056 1275 (N11E04 L=137)   190     18        1 (1040/1254-1076/1294) -45.0 E05-E03 N09-N12
     10 1664 1164 (N21W37 L=178)   473     89        0 (1614/1144-1708/1185)   6.3 W32-W41 N20-N22
     11  398  831 (S30E38 L=103)   360     72        0 ( 371/ 814- 434/ 856)   0.0 E40-E34 S31-S29
     12  716  807 (S24E13 L=128)   344    119        0 ( 672/ 791- 776/ 834)  15.9 E16-E09 S25-S21
     13  824  769 (S23E05 L=136)   318    152        0 ( 767/ 743- 876/ 793)  33.7 E08-E01 S26-S21
     14  662  748 (S29E15 L=126)   331     58        1 ( 639/ 717- 694/ 778)  90.0 E16-E13 S31-S26
     15 1646  583 (S12W54 L=195)   439     39        0 (1629/ 558-1667/ 609)  90.0 W52-W56 S13-S09
     16 1427  571 (S19W38 L=179)   435     55        0 (1405/ 542-1451/ 604) -36.9 W36-W40 S20-S16
     17  878  529 (S37W05 L=146)   405     60        1 ( 843/ 512- 917/ 554)  36.9 W02-W07 S38-S35
     18 1485  426 (S25W51 L=192)   523     91        0 (1442/ 409-1530/ 445)   6.3 W46-W55 S26-S23
     19  961  285 (S52W25 L=166)   516     39        1 ( 935/ 259- 986/ 319) -29.7 W21-W28 S53-S49
  

    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