Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  859 1809 (N49E37 L=254)   573     77        1 ( 830/1780- 902/1825)  36.9 E40-E33 N46-N51
      2  651 1718 (N37E47 L=244)   851    126        2 ( 623/1661- 692/1781) -25.2 E57-E39 N33-N41
      3  967 1685 (N41E19 L=272)   362     70        0 ( 940/1651- 988/1709) -33.7 E22-E16 N38-N42
      4  447 1639 (N28E61 L=231)   487     51        1 ( 423/1620- 470/1666) -20.6 E64-E56 N26-N29
      5 1555 1560 (N44W33 L=325)   572     80        1 (1513/1536-1592/1578)   0.0 W28-W38 N42-N46
      6 1393 1443 (N32W18 L=310)   279     87        0 (1356/1418-1424/1477) -18.4 W14-W20 N30-N33
      7 1073 1367 (N21E03 L=289)   354    104        0 (1036/1337-1122/1403) -21.8 E05-W00 N19-N22
      8  275 1216 (S01E56 L=236)   933    155        2 ( 229/1168- 327/1276) -17.1 E63-E49 S03-N01
      9 1742  622 (S12W61 L=353)  1203    232        1 (1652/ 551-1816/ 679)  42.5 W55-W68 S18-S07
     10  929  565 (S31W02 L=294)   262     60        0 ( 892/ 552- 962/ 580)  26.6 W00-W04 S31-S29
     11  977  368 (S44W11 L=303)   541    139        1 ( 941/ 312-1020/ 409)  90.0 W09-W13 S49-S40
     12 1543  341 (S34W62 L=354)  1131    126        2 (1498/ 318-1586/ 377)  10.3 W56-W69 S36-S31
     13 1133  241 (S51W33 L=325)   361     49        0 (1107/ 217-1159/ 262) -12.5 W28-W37 S53-S49
  

    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