Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  548 1608 (N15E48 L=281)   400     83        0 ( 532/1569- 564/1647) -59.0 E49-E45 N12-N17
      2 1287 1543 (N31W00 L=329)   368     92        0 (1241/1510-1336/1564)   9.5 E02-W04 N29-N32
      3 1694 1514 (N46W34 L=  3)   919    140        1 (1670/1462-1710/1558)  71.6 W31-W36 N41-N50
      4  499 1428 (N03E43 L=285)   502    117        0 ( 480/1393- 528/1483) -63.4 E45-E40 N00-N06
      5 1520 1431 (N32W19 L=348)   384     84        1 (1497/1398-1551/1475) -38.7 W16-W21 N29-N34
      6  697 1381 (N05E29 L=300)   357    109        0 ( 644/1344- 735/1406)   9.5 E31-E25 N03-N05
      7  249 1208 (S15E57 L=272)  2234    470        5 ( 136/1167- 371/1235)  12.3 E69-E45 S18-S11
      8 1819  986 (N16W52 L= 21)   644    110        1 (1798/ 923-1836/1036) -51.3 W49-W54 N12-N18
      9 1752  840 (N05W51 L= 19)   849    202        2 (1694/ 798-1809/ 889)   4.8 W45-W57 N03-N06
     10 1342  794 (S10W24 L=353)   265     94        0 (1322/ 746-1362/ 847) -71.6 W23-W25 S12-S06
     11 1884  788 (N08W68 L= 37)  1150    236        1 (1849/ 754-1935/ 826)  41.2 W62-W74 N05-N12
     12 1507  692 (S10W38 L=  7)   310     54        1 (1477/ 668-1532/ 714)  14.0 W36-W40 S11-S08
     13  819  630 (S35E01 L=328)   240     77        1 ( 776/ 617- 857/ 646)  36.9 E03-W01 S36-S32
     14  985  573 (S33W12 L=341)  1205    281        2 ( 831/ 536-1096/ 631)  40.6 W02-W17 S38-S26
     15 1588  593 (S13W48 L= 17)   314     60        0 (1550/ 575-1618/ 615)  45.0 W46-W50 S14-S10
     16 1430  418 (S27W47 L= 16)  1895    397        1 (1305/ 340-1529/ 492)   2.4 W34-W58 S28-S24
  

    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