Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  641 1721 (N58E13 L=247)  3232    504        6 ( 434/1617- 832/1822)   5.6 E40-W12 N54-N62
      2 1602 1602 (N22W57 L=317)  3834    452       16 (1511/1493-1704/1710) -82.9 W47-W67 N13-N29
      3 1668 1645 (N23W69 L=329)   538     53        1 (1649/1621-1685/1671)  36.9 W65-W73 N21-N24
      4 1176 1530 (N29W22 L=282)   522    205        1 (1103/1455-1239/1600) -76.0 W20-W23 N22-N34
      5 1552 1474 (N16W45 L=305)   357    132        0 (1521/1415-1571/1529)  76.0 W42-W47 N12-N20
      6  278 1386 (N41E52 L=208)   525     52        0 ( 258/1362- 299/1413)  36.9 E55-E49 N39-N42
      7  561  767 (S02E33 L=227)   361     84        0 ( 547/ 725- 575/ 814)  63.4 E34-E31 S04-N00
      8 1609  745 (S27W32 L=292)   253     50        2 (1586/ 719-1629/ 768) -45.0 W30-W33 S28-S25
      9 1329  534 (S33W08 L=268)  1804    527        2 (1114/ 431-1490/ 689)   5.9 E07-W22 S38-S27
     10 1490  427 (S47W20 L=280)   710    155        1 (1429/ 391-1545/ 488)  15.3 W14-W25 S48-S43
     11  391  410 (S21E65 L=195)  1855    128        3 ( 334/ 376- 460/ 451)  -2.6 E80-E58 S23-S18
     12 1068  154 (S58E32 L=228)   542     67        1 (1034/ 136-1098/ 176) -15.9 E35-E26 S60-S55
  

    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