Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1116 1616 (N32E02 L= 46)   276    105        0 (1050/1598-1153/1629)  18.4 E06-W00 N30-N33
      2  311 1470 (N14E57 L=351)  4528    371       13 ( 220/1316- 446/1586)  50.9 E65-E47 N04-N22
      3  463 1514 (N18E45 L=  4)   766    119        2 ( 413/1472- 509/1565)  45.0 E47-E41 N15-N21
      4  519 1469 (N15E39 L= 10)   344     59        2 ( 483/1442- 545/1492)  26.6 E41-E36 N13-N17
      5 1843 1078 (N10W56 L=104)   764    128        1 (1824/1013-1861/1141)  76.0 W53-W57 N05-N13
      6  151 1079 (S10E65 L=343)   741     89        3 ( 130/1033- 167/1125) -78.7 E68-E62 S12-S07
      7  222 1004 (S14E56 L=352)  1094    156        5 ( 188/ 948- 266/1089)  53.1 E60-E51 S17-S09
      8  975  849 (S17E01 L= 48)   654    226        0 ( 926/ 756-1002/ 954) -70.0 E04-W01 S22-S11
      9 1522  762 (S14W35 L= 84)  1115    442        0 (1441/ 619-1599/ 928) -87.0 W30-W39 S22-S03
     10  280  852 (S23E51 L=358)  2933    539        5 (  94/ 773- 550/ 901)   7.1 E84-E28 S27-S18
     11 1214  855 (S13W13 L= 62)   287    118        0 (1155/ 842-1277/ 874)   7.1 W09-W17 S13-S12
     12  574  709 (S30E27 L= 22)  2241    453        6 ( 365/ 626- 727/ 849)  29.5 E46-E16 S36-S19
     13 1286  751 (S18W20 L= 68)   534    272        0 (1190/ 663-1371/ 847)  58.4 W15-W23 S24-S11
     14 1772  746 (S10W56 L=104)   300     59        1 (1738/ 731-1801/ 765)  26.6 W52-W58 S11-S08
     15 1120  407 (S43W18 L= 67)  2505    439        7 ( 914/ 366-1334/ 457)  10.3 W01-W34 S45-S38
     16  802  390 (S49E08 L= 41)   261     44        0 ( 779/ 369- 826/ 410)  45.0 E10-E05 S51-S47
     17 1594  371 (S36W61 L=110)   788     80        0 (1556/ 350-1634/ 395)  11.3 W55-W66 S37-S33
  

    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