Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1015 1897 (N63W10 L=126)   409     51        0 ( 995/1870-1036/1918)  49.4 W06-W13 N59-N66
      2 1199 1552 (N30W16 L=132)  1511    349        1 (1060/1455-1324/1647) -43.0 W07-W22 N22-N37
      3 1447 1584 (N31W35 L=152)   759    182        1 (1373/1535-1529/1632) -31.0 W30-W40 N27-N34
      4  846 1579 (N34E09 L=108)   340     49        2 ( 817/1556- 879/1599) -21.8 E11-E06 N32-N36
      5  554 1476 (N29E30 L= 86)   337    114        0 ( 497/1447- 605/1500)  14.0 E34-E26 N27-N30
      6  438 1362 (N23E38 L= 78)  1089    205        2 ( 400/1305- 507/1437) -66.0 E41-E32 N18-N27
      7  911 1382 (N20E05 L=112)   231     83        0 ( 879/1362- 936/1412)   0.0 E07-E03 N19-N22
      8  176 1349 (N24E68 L= 48)  1641    108        3 ( 141/1298- 203/1400)  71.6 E74-E64 N21-N27
      9  574 1358 (N21E28 L= 89)   220     51        1 ( 548/1335- 595/1377) -45.0 E29-E26 N19-N22
     10  874 1325 (N17E08 L=109)   288    102        3 ( 820/1305- 914/1356)   9.5 E10-E04 N16-N19
     11 1838 1151 (N02W59 L=176)   644    114        2 (1777/1130-1890/1167)  -8.7 W52-W65 N00-N02
     12  958  792 (S15E06 L=111)   320    164        0 ( 934/ 719-1009/ 843)  60.3 E07-E02 S19-S12
     13  792  823 (S13E16 L=101)   200     38        1 ( 769/ 805- 816/ 839)   0.0 E16-E13 S13-S11
     14 1813  811 (S19W57 L=174)   377     56        2 (1782/ 785-1841/ 839)  21.8 W54-W60 S20-S16
     15  851  630 (S25E14 L=103)   596    168        0 ( 830/ 566- 900/ 708) -78.7 E15-E11 S29-S19
     16 1629  533 (S36W46 L=162)  1947    334        1 (1514/ 444-1743/ 660)  45.0 W36-W53 S42-S27
     17 1397  472 (S39W25 L=142)   711    118        1 (1356/ 413-1427/ 525) -58.0 W21-W29 S43-S35
     18 1296  252 (S58W22 L=139)   516     88        1 (1245/ 236-1347/ 268) -12.1 W15-W29 S59-S55
  

    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