Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  705 1662 (N32E39 L=197)   623    105        0 ( 651/1632- 741/1686)   8.1 E42-E34 N29-N34
      2 1020 1572 (N33E13 L=223)   920    195        1 ( 970/1501-1086/1632) -32.5 E17-E06 N29-N36
      3  671 1496 (N20E33 L=203)   387     85        0 ( 631/1480- 719/1512)  26.6 E36-E30 N18-N22
      4 1581 1132 (N16W34 L=270)   584     71        2 (1557/1096-1605/1172)  68.2 W31-W35 N13-N18
      5  168 1069 (S14E67 L=169)   939    127        0 ( 141/1015- 201/1111)  80.5 E71-E61 S17-S11
      6 1449  675 (S13W34 L=270)  1709    466        6 (1328/ 574-1568/ 751)  -2.9 W25-W45 S21-S09
      7 1252  659 (S18W22 L=258)   289     86        1 (1224/ 627-1290/ 691)  59.0 W20-W23 S20-S15
      8  724  661 (S29E12 L=224)   361    149        1 ( 655/ 621- 790/ 683)  37.9 E16-E07 S33-S26
      9  667  338 (S56E14 L=222)   576     68        2 ( 623/ 321- 702/ 362)  26.6 E20-E10 S59-S53
     10 1490  283 (S38W65 L=301)   628     70        0 (1461/ 265-1521/ 309)  -4.1 W58-W73 S38-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