Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1640 1636 (N30W56 L=203)   637    189        0 (1615/1593-1678/1663)  24.0 W51-W60 N27-N32
      2 1361 1648 (N35W33 L=180)   310     53        0 (1334/1631-1396/1662)   0.0 W31-W36 N33-N35
      3  307 1467 (N35E50 L= 97)  2041    438        2 ( 211/1437- 404/1509) -17.7 E62-E40 N32-N40
      4 1815 1426 (N15W64 L=211)   540    112        0 (1793/1397-1844/1445)  15.9 W61-W69 N13-N16
      5  267 1243 (N21E51 L= 96)   527    164        1 ( 234/1209- 307/1274) -29.7 E54-E46 N18-N22
      6  556 1260 (N18E27 L=120)   251     83        0 ( 517/1245- 595/1273) -21.8 E29-E24 N17-N19
      7  193 1156 (N16E58 L= 88)   409     93        1 ( 175/1115- 206/1191) -59.0 E60-E56 N13-N18
      8  244 1138 (N15E53 L= 94)   612     81        2 ( 222/1104- 270/1182)  59.0 E55-E50 N12-N17
      9  137 1133 (N16E65 L= 81)   895    164        4 ( 109/1100- 158/1174)  51.3 E69-E62 N13-N18
     10  648 1125 (N09E22 L=125)   466     97        0 ( 600/1079- 686/1168) -45.0 E24-E19 N06-N12
     11  838 1019 (N01E11 L=136)   382    104        1 ( 797/ 962- 872/1076)  45.0 E14-E08 S02-N04
     12  163 1034 (N09E63 L= 84)   446     76        1 ( 147/1007- 183/1057)  21.8 E64-E59 N07-N10
     13 1861  998 (S11W60 L=206)   361     59        0 (1848/ 964-1878/1026) -90.0 W58-W61 S12-S09
     14 1766  978 (S11W50 L=196)   291     69        1 (1741/ 953-1784/1005) -56.3 W47-W51 S12-S09
     15 1308  696 (S24W14 L=161)   680    229        0 (1187/ 662-1419/ 740)   0.0 W06-W22 S25-S21
     16  438  495 (S26E49 L= 97)  1310    355        2 ( 406/ 377- 489/ 625)  45.0 E56-E41 S34-S18
     17 1668  556 (S37W44 L=191)   323     71        1 (1654/ 526-1685/ 594)  90.0 W42-W45 S39-S34
     18 1073  406 (S41E05 L=142)   468     87        1 (1027/ 386-1112/ 439)   0.0 E08-E01 S42-S38
  

    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