Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  455 1664 (N37E60 L=302)  2072    155        5 ( 423/1598- 493/1739) -15.7 E82-E50 N32-N41
      2  659 1609 (N35E35 L=327)  1760    246        6 ( 558/1565- 768/1643)  10.0 E43-E25 N31-N37
      3  557 1468 (N25E38 L=324)   385    125        1 ( 542/1417- 587/1530)  90.0 E39-E35 N21-N29
      4  625 1405 (N21E31 L=331)   551    201        0 ( 551/1357- 666/1486)  66.0 E35-E28 N17-N26
      5 1240 1394 (N26W11 L= 13)   361    152        0 (1203/1350-1310/1422)  20.6 W08-W16 N23-N28
      6 1023 1320 (N20E03 L=359)   229    100        1 ( 970/1307-1076/1329)   9.5 E06-W00 N18-N20
      7 1457 1193 (N15W26 L= 29)   215     90        0 (1411/1184-1503/1202)   0.0 W23-W29 N14-N16
      8 1767 1097 (N12W52 L= 54)   383     72        1 (1751/1060-1782/1136) -68.2 W50-W53 N09-N14
      9 1447  995 (N03W27 L= 29)   470    156        1 (1390/ 958-1518/1022)  20.6 W23-W31 N00-N04
     10 1756  703 (S12W57 L= 59)  3050    442        3 (1627/ 586-1844/ 897)  57.0 W48-W63 S20-N00
     11  437  592 (S32E41 L=321)  1137    223        1 ( 367/ 509- 522/ 646)  32.3 E51-E31 S38-S26
  

    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