Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1212 1653 (N38E08 L= 14)  2251    574        7 (1092/1457-1312/1831) -35.0 E26-W04 N27-N50
      2  419 1686 (N19E66 L=316)  2715    270        2 ( 360/1635- 471/1747)  49.4 E69-E61 N14-N24
      3  616 1585 (N17E42 L=340)   619    208        1 ( 593/1519- 649/1675) -48.4 E46-E37 N13-N22
      4  511 1590 (N15E50 L=332)   452     72        1 ( 474/1571- 542/1617)  53.1 E53-E47 N12-N16
      5  373 1510 (N07E58 L=324)   654     96        2 ( 324/1454- 409/1546)  56.3 E60-E55 N02-N09
      6  470 1490 (N08E49 L=334)   456     94        1 ( 422/1449- 497/1520)  59.0 E50-E45 N04-N09
      7 1672 1276 (N29W34 L= 56)   625    173        0 (1605/1200-1703/1358) -35.0 W27-W37 N24-N31
      8  772 1040 (S11E15 L=  8)   237     88        0 ( 726/1006- 816/1076)  56.3 E16-E12 S13-S07
      9 1841  913 (N13W58 L= 80)  2575    381        5 (1802/ 757-1887/1049) -39.5 W51-W69 N06-N20
     10  165  926 (S32E60 L=322)   535    134        0 ( 145/ 884- 198/ 968)  50.2 E61-E56 S34-S28
     11 1533  932 (N04W32 L= 54)   227     18        0 (1513/ 912-1550/ 952)  71.6 W30-W33 N02-N05
     12 1278  775 (S12W21 L= 43)   560    207        0 (1202/ 705-1383/ 821)  61.4 W18-W26 S17-S06
     13  190  796 (S39E57 L=325)   779    175        2 ( 149/ 780- 224/ 820)  20.6 E61-E52 S41-S37
     14 1005  795 (S18W05 L= 27)   245     95        1 ( 944/ 784-1050/ 813)  26.6 W01-W07 S19-S16
     15 1026  651 (S26W11 L= 33)   703    176        0 ( 930/ 632-1113/ 676)  31.0 W05-W15 S28-S21
     16 1201  498 (S29W28 L= 50)   381     69        0 (1159/ 473-1237/ 525)   0.0 W24-W31 S30-S28
     17  920  463 (S41W12 L= 34)   263     52        0 ( 899/ 442- 946/ 489)   0.0 W10-W14 S42-S38
     18 1014  457 (S38W19 L= 41)   232     29        0 ( 993/ 441-1042/ 473)  45.0 W16-W20 S39-S36
  

    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