Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  977 1788 (N59W16 L= 47)  1338    307        2 ( 842/1742-1103/1830)  -3.6 W00-W32 N56-N60
      2 1400 1519 (N30W36 L= 67)   604    238        0 (1299/1444-1481/1570) -54.0 W30-W38 N23-N34
      3  530 1473 (N41E32 L=359)   499    164        1 ( 470/1406- 572/1518) -33.7 E39-E27 N36-N44
      4  366 1341 (N34E46 L=345)  2879    477        2 ( 284/1138- 465/1517)  65.4 E52-E39 N21-N45
      5 1223 1457 (N29W21 L= 52)   241     79        0 (1184/1438-1257/1481) -53.1 W18-W22 N27-N31
      6 1838 1337 (N09W68 L= 99)   518     68        2 (1817/1314-1861/1363) -53.1 W66-W72 N07-N11
      7  193 1226 (N27E65 L=326)  2365    228        3 ( 135/1130- 252/1313) -21.4 E85-E55 N21-N32
      8  566 1196 (N23E28 L=  4)   295    135        0 ( 511/1169- 614/1247) -35.5 E31-E24 N20-N26
      9  813 1173 (N18E11 L= 20)   395    205        0 ( 747/1133- 908/1225) -41.6 E14-E05 N14-N22
     10  449 1012 (N12E38 L=354)   458    151        2 ( 406/ 956- 469/1067) -74.1 E40-E35 N09-N16
     11  370  870 (N05E46 L=346)   304     83        0 ( 351/ 836- 395/ 914) -68.2 E46-E43 N02-N07
     12 1580  826 (S15W33 L= 64)   989    197        1 (1496/ 774-1658/ 875)   8.7 W26-W39 S16-S12
     13 1441  797 (S14W22 L= 54)   216     41        0 (1422/ 777-1461/ 818) -56.3 W21-W23 S15-S12
     14  709  598 (S15E27 L=  5)  1144    394        0 ( 602/ 526- 791/ 718)  31.0 E34-E19 S19-S08
     15  431  628 (S11E46 L=345)  1024    168        1 ( 410/ 535- 470/ 692)  63.4 E50-E41 S16-S06
     16  550  622 (S12E37 L=354)   282     61        0 ( 521/ 608- 584/ 639) -26.6 E39-E34 S13-S10
     17 1410  346 (S47W20 L= 51)   877    117        2 (1364/ 316-1464/ 377)   0.0 W15-W25 S49-S43
     18 1203  251 (S51E01 L= 31)  1646    166        5 (1140/ 203-1264/ 310)  22.6 E07-W05 S56-S45
     19 1195  167 (S62E04 L= 27)   654     42        1 (1167/ 150-1223/ 184) -15.9 E07-W00 S65-S58
  

    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