Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  822 1594 (N35E12 L=324)   294     68        1 ( 802/1566- 837/1629)  76.0 E13-E10 N33-N37
      2 1607 1508 (N26W44 L= 20)   506     74        2 (1585/1478-1639/1549) -78.7 W42-W46 N23-N28
      3 1786 1357 (N15W57 L= 33)  1125    191        5 (1726/1282-1845/1418) -48.4 W52-W63 N10-N19
      4 1551 1378 (N17W37 L= 13)   269     72        2 (1519/1360-1584/1397) -26.6 W34-W39 N16-N18
      5 1627 1289 (N11W41 L= 17)   303     63        1 (1607/1263-1655/1321) -76.0 W39-W43 N09-N13
      6 1243 1149 (N04W14 L=350)  1565    380        4 (1111/ 967-1316/1269) -62.1 W06-W17 S06-N11
      7  819 1106 (N03E12 L=324)   519    206        1 ( 783/1049- 874/1162) -56.3 E14-E08 N00-N06
      8  334  916 (S05E47 L=289)   342     49        0 ( 312/ 893- 359/ 942) -56.3 E48-E44 S06-S03
      9 1774  820 (S17W52 L= 28)   467     62        0 (1745/ 781-1803/ 857)  51.3 W49-W55 S19-S14
     10  259  757 (S14E57 L=279)  1657    296        3 ( 218/ 662- 318/ 825)  32.7 E63-E49 S19-S10
     11 1505  438 (S42W37 L= 13)  1870    366        2 (1344/ 382-1675/ 518)  13.0 W22-W49 S46-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