Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1520 1565 (N30W40 L=133)  1994    381        3 (1361/1460-1651/1673) -39.8 W30-W48 N22-N38
      2 1697 1598 (N32W62 L=154)  2132    169        3 (1628/1542-1763/1639) -35.0 W54-W69 N27-N34
      3  794 1480 (N27E13 L= 80)   316     89        0 ( 767/1438- 827/1510)  59.0 E14-E10 N24-N29
      4 1066 1453 (N24W05 L= 98)   344    145        0 (1002/1403-1101/1489) -39.8 W00-W06 N20-N26
      5  368 1372 (N23E45 L= 47)   811     85        1 ( 300/1323- 420/1425)  31.0 E51-E41 N20-N26
      6 1241 1358 (N17W15 L=108)   486    148        1 (1174/1334-1311/1386)   0.0 W11-W20 N15-N18
      7 1242  928 (S09W13 L=105)   723    194        0 (1170/ 865-1318/ 994)  33.7 W08-W17 S12-S05
      8 1204  688 (S23W10 L=102)   516    141        1 (1183/ 619-1221/ 762)  90.0 W08-W11 S28-S18
      9 1716  652 (S28W51 L=144)   498     66        1 (1689/ 623-1739/ 680) -53.1 W48-W54 S29-S25
     10 1639  486 (S39W50 L=142)   992     96        2 (1615/ 430-1669/ 540) -82.9 W45-W54 S43-S35
     11  709  297 (S49E36 L= 56)   574     60        1 ( 686/ 266- 747/ 324)   7.1 E40-E31 S51-S46
     12  559  279 (S49E56 L= 37)   540     69        1 ( 525/ 264- 602/ 300) -14.0 E61-E51 S50-S46
  

    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