Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  919 1737 (N51W19 L=142)   653    125        1 ( 882/1680- 953/1791)  30.3 W13-W25 N47-N54
      2 1295 1712 (N36W43 L=166)   582    138        1 (1236/1675-1360/1738)  -4.8 W37-W49 N34-N37
      3 1245 1595 (N31W32 L=155)   291     38        1 (1217/1574-1275/1610) -33.7 W30-W33 N29-N32
      4  238 1343 (N43E56 L= 67)   381     67        0 ( 212/1312- 270/1366)  14.0 E59-E51 N41-N44
      5  798 1351 (N29E05 L=118)   221     67        0 ( 764/1338- 828/1364) -45.0 E06-E02 N27-N30
      6 1723 1268 (S01W50 L=173)  1056    176        2 (1649/1218-1797/1302) -32.5 W44-W55 S05-N02
      7  332 1075 (N24E44 L= 79)  2197    418        4 ( 183/ 987- 508/1209) -34.2 E60-E31 N15-N34
      8  618  983 (N13E24 L= 99)   207    120        0 ( 603/ 938- 644/1044) -90.0 E25-E23 N09-N16
      9 1804  942 (S21W49 L=172)   507     94        1 (1775/ 896-1824/ 994)  51.3 W45-W51 S23-S18
     10 1912  936 (S26W63 L=186)   634     45        1 (1883/ 918-1942/ 952) -10.3 W58-W69 S26-S23
     11  259  750 (N07E58 L= 65)  2989    345        9 ( 192/ 628- 345/ 850)  81.3 E68-E50 N00-N13
     12 1262  730 (S17W06 L=129)   659    314        0 (1158/ 637-1381/ 826)  13.2 E02-W15 S20-S14
     13 1602  664 (S32W27 L=150)  1617    376        2 (1498/ 555-1738/ 807)  19.2 W17-W40 S36-S26
     14  897  759 (S06E14 L=110)   980    294        1 ( 830/ 661-1019/ 802)   4.1 E20-E06 S10-S04
     15  478  594 (S06E45 L= 78)   371     76        1 ( 454/ 566- 513/ 623)  18.4 E47-E41 S07-S05
     16  586  510 (S14E41 L= 82)   829    184        2 ( 532/ 439- 633/ 592)  30.3 E48-E35 S16-S08
     17  837  168 (S44E49 L= 74)   819    109        0 ( 789/ 142- 880/ 208)   7.6 E56-E41 S45-S41
  

    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