Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1203 1766 (N47W03 L=307)   389     62        1 (1165/1748-1245/1779)  15.9 W00-W07 N45-N48
      2 1334 1699 (N43W15 L=318)   421     67        1 (1306/1668-1367/1725)  31.0 W12-W17 N40-N45
      3  630 1578 (N25E34 L=269)   362    101        0 ( 600/1521- 656/1610) -39.8 E37-E30 N21-N26
      4  479  704 (S30E35 L=268)   241     47        0 ( 450/ 689- 507/ 722)  21.8 E37-E32 S31-S29
      5  546  642 (S34E30 L=273)   307     90        0 ( 498/ 618- 587/ 667)  20.6 E34-E26 S35-S32
      6 1430  580 (S27W34 L=338)   324     75        2 (1397/ 554-1467/ 603)  36.9 W32-W36 S28-S24
      7  740  459 (S45E15 L=289)   262     39        2 ( 722/ 426- 757/ 482)  90.0 E16-E13 S47-S42
      8 1264  441 (S38W27 L=331)   376     85        0 (1242/ 409-1296/ 467)  90.0 W25-W30 S40-S36
  

    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