Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  406 1578 (N31E59 L=179)  1488    149        2 ( 341/1526- 464/1620)  -5.0 E74-E51 N27-N33
      2 1681 1467 (N36W51 L=289)   630     78        1 (1647/1429-1710/1498) -90.0 W47-W53 N32-N37
      3 1692 1385 (N30W49 L=287)   419     39        1 (1668/1362-1709/1409) -45.0 W47-W51 N28-N31
      4 1658  802 (S07W45 L=283)  1002    199        2 (1619/ 730-1679/ 898) -63.4 W41-W47 S11-S01
      5  584  822 (S15E27 L=211)   818    263        1 ( 485/ 764- 713/ 857) -10.0 E34-E17 S17-S12
      6  873  670 (S22E07 L=231)   218    122        0 ( 854/ 621- 895/ 712) -78.7 E08-E04 S24-S18
      7 1174  582 (S25W14 L=252)   223     96        0 (1131/ 569-1220/ 601)   0.0 W11-W17 S25-S23
      8 1382  520 (S28W31 L=268)   310     45        1 (1356/ 500-1409/ 537)  33.7 W28-W32 S29-S26
      9  398  461 (S43E53 L=185)   634     60        0 ( 380/ 424- 418/ 493)  32.5 E58-E47 S46-S39
  

    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