Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1021 1674 (N31E20 L=139)  1372    339        0 ( 956/1559-1116/1730) -86.2 E24-E13 N22-N37
      2  505 1683 (N20E57 L=102)   628    107        2 ( 460/1654- 543/1716)  38.7 E60-E54 N17-N22
      3  413 1627 (N15E62 L= 97)  1365    120        4 ( 376/1560- 441/1684) -48.8 E68-E56 N10-N18
      4  555 1632 (N17E49 L=110)   846    185        2 ( 469/1598- 630/1663)  12.1 E56-E41 N14-N19
      5 1412 1626 (N41W08 L=167)   284     51        0 (1383/1603-1436/1649) -11.3 W05-W10 N39-N42
      6 1052 1073 (S04W00 L=159)   280     72        1 (1037/1031-1068/1116) -78.7 E01-W01 S06-S01
      7  633 1067 (S14E24 L=135)   215     80        1 ( 593/1048- 673/1090)  31.0 E26-E21 S15-S12
      8 1787  986 (N13W50 L=209)   342     49        0 (1755/ 967-1811/1004)   0.0 W46-W52 N11-N13
      9 1817  879 (N08W56 L=215)   739     64        2 (1791/ 848-1850/ 916) -45.0 W53-W59 N06-N09
     10  394  787 (S36E38 L=121)   720    165        0 ( 328/ 748- 448/ 854)  36.0 E45-E33 S39-S30
     11  999  825 (S19W04 L=163)   201     43        0 ( 973/ 805-1024/ 845) -14.0 W01-W05 S19-S18
     12 1821  734 (N01W63 L=222)   468     36        1 (1794/ 718-1841/ 755)  26.6 W60-W66 N00-N02
     13 1709  666 (S07W55 L=214)   846     97        1 (1670/ 630-1747/ 707)   0.0 W49-W60 S08-S04
     14 1373  663 (S17W31 L=190)   320    118        1 (1324/ 617-1402/ 688)   0.0 W27-W34 S19-S15
     15 1494  382 (S28W56 L=215)  1248    237        2 (1379/ 355-1597/ 415)  33.7 W46-W63 S32-S22
     16  742  276 (S64W08 L=167)  1522    150        3 ( 692/ 236- 806/ 338) -13.2 W00-W17 S67-S59
  

    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