Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  566 1696 (N29E48 L=334)  1459    222        4 ( 461/1657- 636/1738)   3.4 E59-E41 N25-N33
      2  377 1619 (N22E61 L=320)   561     77        0 ( 345/1586- 408/1647)  59.0 E63-E59 N19-N24
      3  476 1609 (N22E50 L=331)   511    101        0 ( 434/1564- 515/1644)  63.4 E52-E47 N18-N24
      4 1297 1347 (N17W11 L= 33)   319     57        1 (1268/1321-1322/1369) -18.4 W09-W13 N15-N18
      5  944  752 (S24E00 L= 21)   771    195        1 ( 859/ 731-1047/ 771)  16.7 E05-W05 S25-S21
      6 1712  640 (S15W56 L= 77)   407     20        0 (1684/ 622-1738/ 654)   9.5 W52-W58 S16-S14
      7 1543  474 (S29W50 L= 72)   537     68        0 (1512/ 440-1575/ 508)  90.0 W48-W51 S31-S25
      8  651  350 (S57E21 L=  1)   541     46        0 ( 629/ 324- 678/ 381) -68.2 E24-E17 S59-S54
  

    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