Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1490 1601 (N27W44 L=238)  1185    205        2 (1417/1560-1557/1639)   8.1 W37-W51 N25-N28
      2 1341 1558 (N27W32 L=226)   306     81        0 (1324/1520-1365/1588)  45.0 W29-W34 N25-N29
      3  835 1558 (N37E03 L=191)   771    129        2 ( 761/1524- 895/1583) -24.0 E08-W01 N34-N39
      4 1588 1356 (N11W42 L=236)   308     49        0 (1557/1340-1619/1374) -21.8 W39-W44 N09-N12
      5  574 1278 (N23E24 L=170)   807    189        0 ( 494/1229- 654/1332) -38.7 E30-E19 N19-N28
      6 1272 1265 (N11W19 L=213)   292     60        1 (1243/1240-1310/1294) -45.0 W16-W20 N09-N12
      7  931 1245 (N15E02 L=192)   610    177        0 ( 869/1183- 974/1286) -51.3 E06-W00 N11-N17
      8  213 1223 (N26E55 L=139)   538    105        0 ( 185/1193- 253/1242)  18.4 E58-E50 N24-N27
      9 1759 1215 (S00W52 L=246)   335     62        0 (1733/1198-1787/1233)   0.0 W49-W55 S01-N00
     10  772  422 (S31E30 L=165)   534     99        0 ( 724/ 394- 815/ 455)  14.0 E33-E25 S32-S29
     11  961  213 (S52E28 L=166)   435     47        1 ( 933/ 195- 984/ 235) -63.4 E30-E24 S54-S50
  

    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